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ABSTRACT 


The  Department  of  Defense  seeks  to  improve  space-based  imaging  capabilities  by 
building  larger  space  telescopes  designed  for  higher  orbits  in  order  to  provide  timely 
accurate  imagery  to  the  military  and  national  decision  makers.  Satellites  have  mass  and 
volume  constraints  due  to  the  limited  size  of  the  launch  vehicle.  Larger  aperture  imagery 
satellites  require  lightweight,  deployable,  segmented  primary  mirrors  to  meet  mass  and 
volume  constraints.  Lightweight  segmented  mirrors  exist  but  lack  the  structure  and  mass 
necessary  for  high-quality  optical  performance.  This  research  investigates  using  a 
deformable  mirror  with  a  large  lightweight  active  primary  mirror  to  improve  optical 
performance. 

Control  techniques  are  applied  to  simultaneously  manipulate  an  active  primary 
mirror  segment  and  deformable  mirror  in  a  woofer-tweeter  configuration.  Computer 
simulations  and  experiments  using  woofer-tweeter  control  techniques  show  reduction  of 
residual  wavefront  error  attributed  to  the  primary  mirror.  A  woofer-tweeter  gradient 
control  technique  was  shown  to  experimentally  reduce  wavefront  error  by  24%  using  the 
Naval  Postgraduate  School’s  3-m  diameter  segmented  mirror  telescope  test  bed.  The 
addition  of  a  deformable  mirror  to  an  active  space  telescope  in  a  woofer -tweeter 
configuration  is  a  technically  feasible  option  for  a  large-aperture  telescope  design,  and 
can  reduce  primary  mirror  requirements  potentially  reducing  cost. 
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I. 


INTRODUCTION 


The  National  Security  Space  Strategy  seeks  to  improve  U.S.  space  capabilities 
that  deliver  timely  and  accurate  space  services  to  both  the  military  and  national  decision 
makers  [1].  High-resolution,  space-based  imagery  is  one  of  these  capabilities.  Imagery 
collection  timeliness,  however,  is  constrained  by  orbital  mechanics  and  the  total  number 
of  satellites.  This  research  develops  technologies  to  improve  large-aperture,  high-altitude 
imagery  satellites  that  can  improve  space-based  imagery  timeliness  using  fewer  satellites 
while  maintaining  current  imagery  accuracy. 

Satellites  have  mass  and  volume  constraints  due  to  the  limited  size  of  the  launch 
vehicle.  Larger  aperture  imagery  satellites  require  lightweight  deployable  primary 
mirrors  to  meet  mass  and  volume  constraints.  Lightweight  primary  mirrors  exist  but  lack 
the  structure  and  mass  that  improve  the  mirror  optical  surface  performance,  which 
improves  the  image  quality.  Using  a  deformable  mirror  along  with  a  lightweight  primary 
mirror  can  allow  for  the  correction  of  surface  errors  common  with  large  lightweight 
primary  mirrors.  This  could  allow  the  development  of  high-altitude  imagery  satellites, 
providing  persistent  high-resolution  imagery  from  space.  The  addition  of  a  deformable 
mirror  also  provides  a  robust  correctable  design  that  can  reduce  the  primary  mirror 
optical  requirements,  reduce  optical  testing  during  manufacturing,  and  correct  for  on- 
orbit  external  disturbances. 

The  goal  of  this  research  is  to  investigate  the  feasibility  of  using  a  deformable 
mirror  with  current  state-of-the-art  large  aperture  active  mirrors  recently  developed  for 
segmented  space  telescopes  at  the  Naval  Postgraduate  School  Segmented  Mirror 
Telescope  Lab.  This  research  applies  and  describes  control  techniques  used  to 
simultaneously  control  a  lightweight  active  primary  mirror  segment  and  deformable 
mirror  to  reduce  the  residual  wavefront  error  attributed  to  the  primary  mirror.  The 
segmented  space  telescope  model  and  control  techniques  are  described  and  simulated 
results  are  presented.  The  control  techniques  are  experimentally  validated  using  the  Naval 
Postgraduate  School’s  3-m  diameter  segmented  mirror  telescope  (SMT)  test  bed. 
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A.  PROBLEM  SPACE 

Current  imagery  satellites  operate  in  low  Earth  orbit  (LEO)  and  use  large 
monolithic  primary  mirrors  to  achieve  high-resolution  images.  The  LEO  limits  the  time  a 
satellite  is  in  view  of  a  ground  target  and  the  constellation  size  limits  the  frequency  the 
target  is  revisited.  Raising  the  orbit  of  the  satellite  increases  the  time  the  satellite  is  in 
view  of  the  ground  target.  To  achieve  similar  ground  resolution  at  a  higher  altitude, 
however,  the  primary  mirror’s  diameter  must  increase.  The  primary  mirror  size  required 
for  a  medium  Earth  orbit  (MEO)  or  geostationary  Earth  orbit  (GEO)  high-resolution 
imagery  satellite  is  greater  than  the  diameter  of  available  launch  vehicle  fairings. 

The  launch  vehicle  volume  and  mass  constraints  impact  the  space  telescope 
primary  mirror  design.  The  volume  constraint  forces  a  trade  between  a  monolithic 
primary  mirror  and  a  segmented  mirror.  A  monolithic  primary  mirror  offers  a  simpler 
design,  but  a  deployable-segmented  mirror  allows  a  primary  mirror  with  a  diameter  larger 
than  the  launch  vehicle  fairing  diameter.  The  mass  constraint  impacts  primary  mirror 
stiffness  and  primary  mirror  optical  surface  figure.  Past  space  telescope  primary  mirrors 
depended  on  stiffness  achieved  through  structure  and  mass  to  maintain  their  optical 
surface  quality.  However,  a  constraint  on  mass  requires  larger  primary  mirrors  to  reduce 
density,  therefore  reducing  stiffness,  and  potentially  reducing  surface  figure  performance. 

Current  large  lightweight  mirror  technologies  do  not  meet  surface  error 
requirements  for  a  visible  telescope.  This  research  proposes  the  addition  of  a  deformable 
mirror  to  a  space  telescope  optical  assembly  in  a  woofer-tweeter  configuration  to 
compensate  for  these  surface  errors.  Researchers  considered  active  mirror  systems  for 
space  telescopes  since  the  early  developments  of  adaptive  optics.  Howell  [2]  provides  an 
overview  of  early  advances  in  large  space  telescope  optical  control.  This  study  uses  a 
deformable  mirror  and  active  primary  mirror  together  to  correct  residual  surface  errors  of 
the  active  replicated  primary  mirrors.  This  research  builds  upon  several  fields  of  study 
including  adaptive  optics,  telescope  and  mirror  design,  structures,  and  controls.  Several 
key  pieces  of  literature  are  discussed. 
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Babcock  [3]  introduced  the  modem  concept  of  adaptive  optics  in  1953,  and  an 
early  implementation  was  adding  tip/tilt  compensation  to  the  Wilson  Observatory  in  1956 
to  compensate  for  atmospheric  scintillation.  Adaptive  optics  technologies  continued  to 
develop  through  the  1990s  primarily  for  defense  purposes.  Hardy  [4],  Greenwood,  and 
Primmerman  [5]  provide  a  concise  overview  of  the  early  history  of  adaptive  optics.  In 
the  1980s  Nelson,  Mast,  and  Faber  [4]  recommended  building  a  ten-meter  telescope  using 
a  segmented  primary  mirror.  This  recommendation  resulted  in  the  development  of  the 
Keck  Observatory  included  the  necessary  wavefront  sensing  and  control  systems. 
Chanan  et  al.  provide  an  overview  of  the  techniques  used  to  actively  phase  the  primary 
mirror  segments  [5]-[8].  The  success  of  ground-based  segmented  telescopes  furthered 
research  in  segmented  space  telescopes. 

Lane  et  al.  [9]  provide  an  overview  of  past  efforts  to  develop  deployable  space 
telescope  technology.  This  includes  the  development  of  several  ground-based  technology 
demonstrations  and  experiments  from  the  late  1980s  to  the  early  2000s.  However,  the 
main  development  in  segmented  space  telescope  technology  is  the  6.5  meter  18  segment 
infrared  James  Webb  Space  Telescope  (JWST).  The  JWST  required  new  wavefront 
sensing  techniques  like  dispersed  fringe  sensing  presented  by  Shi  et  al.  [10]  and  control 
techniques  presented  by  Redding  et  al.  [11]. 

Researchers  continue  to  develop  technologies  to  support  visible  wavelength 
segmented  space  telescopes  with  a  focus  on  development  of  advanced  lightweight 
primary  mirror  segments.  Current  advanced  mirror  technologies  under  consideration 
include  actuated  hybrid  mirrors  (AHM)  described  by  Hickey  et  al.  [12],  [13],  and  carbon 
fiber  mirror  technologies  presented  by  Romeo  et  al.  [14]  and  Steeves  et  al.  [15].  Both 
technologies  are  lightweight  replicated  optics  and  can  include  embedded  actuators  to 
compensate  for  surface  figure  errors.  However,  current  AHM  mirrors  and  carbon  fiber 
mirrors  struggle  to  achieve  optical  surface  performance  of  traditional  monolithic  glass 
mirrors. 
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B. 


RESEARCH  RESOURCES 


To  address  the  design  challenges  of  a  large  aperture  deployable  space  telescope 
the  Department  of  Defense  (DOD)  developed  the  Segmented  Mirror  Demonstrator 
(SMD)  telescope.  The  SMD  telescope  developed  critical  segmented  mirror  telescope 
technologies  including  lightweight  active  mirror  segments,  deployment  mechanisms,  and 
wavefront  sensing  and  control  technologies  using  a  ground-based  demonstration.  After 
program  completion  the  DOD  transferred  the  SMD  telescope  to  the  Naval  Postgraduate 
School  and  renamed  the  telescope  the  Segmented  Mirror  Telescope  (SMT). 

The  SMT  is  a  3 -meter  deployable-segmented  telescope  with  six  hexagonal  active 
1 -meter  segments.  The  active  mirror  segments  have  many  advantages  over  traditional 
mirrors.  The  segments  are  lower  cost,  lower  mass,  and  actively  correct  surface  errors 
using  a  wavefront  sensing  and  control  system.  The  optical  telescope  assembly  includes 
the  control  system  as  part  of  the  optical  prescription.  The  SMT  has  992  actuator  inputs 
including  936  surface  parallel  face  sheet  actuators,  18  fine  control  actuator,  36  coarse 
control  actuators,  and  a  tip/tilt  fast  steering  mirror.  The  SMT  sensors  include  a  Shack- 
Hartmann  wavefront  sensor,  phase  diversity  sensor,  and  segment  edge  sensors.  Figure  1 
shows  the  SMT  located  at  the  Naval  Postgraduate  School. 


The  SMT  active  segments  have  residual  surface  errors  after  correction.  The 
segment  manufacturing  process  caused  residual  surface  error  consisting  of  high  spatial 
frequency  print-through.  The  fixture  points  on  the  back  of  the  segments  also  caused 
larger  low  spatial  frequency  errors.  Figure  2  shows  a  two-dimensional  plot  of  the  SMT 
segment  optical  path  difference  measured  in  waves  (632.8  nm)  after  applying  a 
correction.  The  high  spatial  frequency  print-through  error  caused  by  the  segment  ribs, 
low  spatial  frequency  error  caused  by  the  fine  control  actuator  attachment  points,  and 
mirror  surface  defects  are  apparent. 


SMT  Scgcmcnt  Residual  WFE 
RMS  WFE  =  0.313 
Pcak-to-Vallcy  WFE  =  2.4023 


Figure  2.  SMT  segment  residual  surface  error,  0.32468  waves  RMS  wavefront 
error  and  2.5464  waves  peak-to-valley  wavefront  error,  from  [16] 

For  this  research  the  SMT  test  bed  was  modified  to  include  a  deformable  mirror  in 
the  optical  path.  This  is  similar  to  a  woofer-tweeter  adaptive  optics  system  used  to 
compensate  for  atmospheric  turbulence  in  astronomical  telescopes  and  high-energy  lasers 
systems.  In  these  systems  the  woofer  deformable  mirror  compensates  for  low  spatial  and 
temporal  frequency  disturbances,  and  the  tweeter  deformable  mirror  compensates  for 
high  spatial  and  temporal  frequency  disturbances.  However,  a  space-based  woofer- 

tweeter  system  does  not  have  to  correct  atmospheric  high  temporal  frequency 
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disturbances,  but  rather  static  errors  and  low  frequency  disturbances  caused  by  a 
thermally  dynamic  environment  [17],  [18].  In  this  system  the  SMT  active  segment 
represents  the  woofer  and  the  additional  deformable  mirror  represents  the  tweeter. 

In  addition  to  active  segmented  mirror  technologies,  NPS  conducts  research  in 
lightweight  monolithic  mirrors  with  the  intent  of  reducing  mass  and  manufacturing  cost. 
Figure  3  shows  a  1 -meter  carbon  fiber  mirror  from  Composite  Mirror  Applications 
(CMA)  delivered  to  NPS  in  March  2012.  The  mirror  surface  shows  significant 
manufacturing  print  through  (quilting)  on  the  optical  surface  (see  Figure  4).  Despite  the 
poor  optical  quality,  the  relative  low  cost  and  short  manufacturing  time  is  an  attractive 
quality  if  the  optical  performance  can  be  improved  or  corrected.  Adding  a  deformable 
mirror  to  a  passive  carbon  fiber  primary  telescope  can  also  improve  optical  performance. 


Figure  3.  Carbon  fiber  reinforced  polymer  mirror,  1 -meter  without 

optical  coating 
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Figure  4.  Carbon  fiber  reinforced  polymer  mirror  surface  error,  RMS 
wavefront  error  of  0.472  waves,  peak-to-valley  wavefront  error  of  3.443  waves 

C.  RESEARCH  OBJECTIVES 

The  objective  of  this  research  is  to  advance  control  techniques  for  woofer-tweeter 
space  telescopes.  This  research  demonstrates  wavefront  error  improvement  on  the  SMT 
test  bed  using  new  control  techniques.  The  addition  of  active  optics  to  the  telescope  and 
new  control  techniques  may  enable  current  and  future  lightweight  primary  mirror 
designs.  The  specific  research  objectives  are  to: 

•  Develop  an  analytical  segmented  space  telescope  model  with  an  additional 
deformable  mirror  in  the  optical  telescope  assembly 

•  Develop  control  algorithms  to  simultaneously  control  an  active  primary 
mirror  and  deformable  mirror  to  reduce  wavefront  error 

•  Experimentally  validate  the  space  telescope  model  and  control  techniques 

•  Reduce  the  surface  error  of  the  SMT  segment  by  controlling  the  segment 
and  deformable  mirror  simultaneously 

D.  DISSERTATION  OUTLINE 

Chapter  II  provides  background  information  for  this  research  including  an 
overview  of  remote  sensing.  The  background  includes  an  overview  of  large  aperture 
mirror  technologies  and  limitations  of  the  current  technology.  The  second  chapter  also 
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includes  a  review  of  the  history  of  segmented  telescopes  and  recent  developments  in 
segmented  space  telescopes.  An  overview  of  adaptive  optics,  the  SMT  test  bed,  and 
previous  SMT  research  is  provided. 

Chapter  III  describes  the  analytical  space  telescope  modeling  approach,  and  the 
placement  of  the  deformable  mirror  in  the  SMT  optical  telescope  assembly.  In  Chapter  III 
the  analytical  model  of  the  SMT  and  the  additional  deformable  mirror  are  described,  as  is 
the  combination  of  these  models  to  create  a  dual  deformable  mirror  model. 

Chapter  IV  describes  using  woofer-tweeter  control  techniques  to  simultaneously 
control  the  SMT  active  primary  segment  and  the  deformable  mirror.  A  new  woofer- 
tweeter  gradient  control  technique  is  proposed  using  singular  value  weighting.  Chapter 
IV  also  introduces  a  technique  to  iteratively  update  the  influence  matrix  to  improve 
mirror  performance  when  there  is  uncertainty  in  the  influence  matrix  or  a  non-  linear 
response  in  the  SMT  AHM  mirror. 

Chapter  V  presents  computer-simulated  results  of  the  control  techniques 
described  in  Chapter  IV.  The  SMT  coarse  alignment,  coarse  phasing,  and  fine  phasing 
simulations  demonstrate  the  need  to  improve  the  SMT  segment  surface  error.  The 
woofer-tweeter  simulations  show  the  feasibility  of  adding  a  deformable  mirror  to  the 
optical  telescope  assembly  and  predict  the  performance  of  the  different  woofer-tweeter 
controllers.  The  simulated  results  show  the  iterative  constrained  minimization  controller 
and  global  gradient  controller  outperform  the  traditional  woofer-tweeter  control 
approaches.  The  simulation  results  also  show  that  controlling  the  woofer  and  tweeter 
simultaneously  outperforms  controlling  them  serially. 

Chapter  VI  describes  the  experimental  system  and  presents  the  experimental 
results.  The  results  are  compared  to  the  simulated  results  and  show  that  the  woofer- 
tweeter  outperformed  the  SMT  segment  performance  and  outperform  operating  the 
mirrors  serially.  The  SMT  AHM  actuator  response  and  SMT  mirror  surface  response  is 
experimentally  characterized  and  the  results  show  operating  regions  with  a  non-linear 
response.  Finally,  the  adaptive  influence  matrix  technique  showed  a  small  improvement 
in  residual  wavefront  error  when  using  an  uncertain  influence  matrix. 
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Chapter  VII  provides  a  summary  of  the  dissertation  contributions  with 
recommendations  for  future  work. 
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II.  BACKGROUND 


This  dissertation  develops  control  techniques  for  an  active  primary  mirror  space 
telescope  with  an  additional  deformable  mirror.  This  work  requires  an  understanding  of 
several  fields  of  study  including  optics,  structures,  and  controls.  This  background 
provides  an  overview  of  the  technologies  that  support  the  development  of  future  large 
aperture  space  telescopes.  The  chapter  begins  with  an  overview  of  space-based  remote 
sensing  characteristics  followed  by  a  review  of  space  telescope  mirror  technologies.  The 
chapter  then  provides  background  on  adaptive  optics  and  discusses  segmented  mirror 
telescopes.  The  chapter  finishes  with  a  background  on  segmented  mirror  telescope 
research  conducted  at  the  Naval  Postgraduate  School. 

A.  SPACE-BASED  EARTH  REMOTE  SENSING 

Spatial  resolution,  spectral  resolution,  and  temporal  resolution  are  performance 
characteristics  of  space-based  remote  sensing  systems.  Payload  and  satellite  designers 
trade  these  characteristics  to  optimize  the  performance  for  the  mission.  The  satellite 
payload  spatial  resolution  depends  on  angular  resolution  and  is  determined  by  the 
Rayleigh  criterion,  Equation  (2.1),  where  X  is  the  wavelength  and  D  is  the  aperture 
diameter.  This  spatial  resolution  translates  to  a  ground  resolution  using  the  small  angle 
approximation,  Equation  (2.2),  where  R  is  the  range  from  the  target  ground  resolution  is 
the  smallest  resolvable  spot  size  on  the  Earth. 

6  =  1.22—  (2.1) 

D 

Ground  Resolution  =  6  - R  (2.2) 

The  satellite  orbit  constrains  the  temporal  resolution.  Equation  (2.3)  determines 
the  orbital  period  of  a  satellite  in  circular  orbit  where  r  is  the  orbit  radius  and  /u  is  the 
standard  gravitational  parameter.  Therefore  a  satellite  at  a  higher  circular  orbit  has  a 
longer  view  of  a  target,  but  a  worse  ground  resolution  than  a  system  with  the  same  spatial 
resolution  at  a  lower  orbit.  Temporal  resolution  also  relates  to  the  revisit  rate  of  the 
satellite  over  a  specific  ground  target. 
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T  =  2n 


(2.3) 


v  M  y 

The  time  a  satellite  is  in  continuous  view  of  a  target  on  a  single  orbit  is  the 
satellite  dwell  time.  In  LEO  dwell  times  are  measured  in  minutes  and  depend  on  both  the 
orbit  and  the  off  nadir  pointing  capability  of  the  satellite.  As  the  satellite  orbit  approaches 
a  geosynchronous  orbit  the  dwell  time  approaches  infinity  meaning  that  there  is  persistent 
coverage  of  the  target  from  a  single  satellite. 

To  increase  dwell  time  and  maintain  the  same  ground  resolution  as  a  LEO 
imagery  satellite  the  angular  resolution  must  decrease  by  increasing  the  aperture 
diameter.  New  telescope  technologies  like  lightweight  large  aperture  replicated  mirrors 
and  segmented  mirror  telescopes  can  enable  future  large  aperture  remote  sensing 
telescopes  that  operate  in  high  altitude  orbits. 

B.  SPACE  TELESCOPE  MIRROR  TECHNOLOGIES 

This  section  provides  an  overview  of  space  telescope  mirrors  including  existing 
on-orbit  mirror  technologies  and  new  technologies  developed  for  large  aperture  space 
telescopes.  Advantages  and  disadvantages  of  each  technology  are  presented  as  well  as 
physical  properties.  This  survey  gives  the  reader  an  overview  of  current  mirror 
capabilities  and  a  historical  perspective. 

New  technologies  include  large  aperture  replicated  optics  and  large  lightweight 
active  mirrors.  In  most  cases  active  mirrors  are  also  replicated  optics.  Replicated  mirrors 
are  formed  using  a  mandrel  where  the  mirror  is  a  negative  spatial  image  of  the  mandrel. 
This  process  produces  identical  copies  where  the  surface  quality  closely  matches  the 
quality  of  the  mandrel.  Traditional  grinding  and  polishing  mirror-manufacturing 
processes  result  in  an  original  optical  element.  An  active  mirror  includes  actuators  that 
apply  controlled  forces  to  change  the  shape  of  the  mirror.  Active  mirrors  can  be 
replicated  or  manufactured  using  grinding  and  polishing  methods. 

Traditional  grinding  and  polishing  mirror-manufacturing  methods  can  take  years. 
For  example,  the  18  JWST  segments  took  8  years  to  complete  manufacturing  [19].  Large 
replicated  mirrors  are  desirable  because  they  offer  a  lower  cost,  high  volume,  and  shorter 
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manufacturing  process  over  traditional  grinding  and  polishing  methods  [20].  Segmented 
mirrors  require  multiple  copies  of  the  same  mirror  making  replicated  mirrors  desirable. 

Hickey,  Barbee,  Ealey,  and  Redding  make  the  argument  for  actively  controlled 
lightweight  mirrors  for  space  telescopes  [21].  They  assert  that  reduction  in  mass  and  the 
ability  to  correct  for  surface  figure  error,  testing  errors,  and  on  orbit  degrading  effects 
outweighs  increased  complexity.  This  technology  has  the  potential  to  reduce  program 
risk  of  large  primary  optics  by  relaxing  fabrication  requirements,  testing  requirements, 
integration  requirements  and  reducing  overall  schedule.  Reducing  mass  compared  to 
conventional  optics  can  also  reduce  system  cost  in  addition  to  cost  savings  associated 
with  reduced  risk.  These  savings  must  be  traded  against  the  complexity  of  adding  a 
wavefront  sensing  and  control  system  to  actively  correct  the  mirror  surface. 

1.  Ultra  Low  Expansion  Glass  Mirrors 

Ultra  Low  Expansion  (ULE)  glass  is  the  standard  for  primary  mirror  construction 
for  optical  telescopes.  The  Hubble  Space  Telescope  (HST)  primary  mirror  began 
development  in  the  late  1970s,  launched  in  1990,  and  remains  the  largest  primary  mirror 
in  space  at  2.4-meter  diameter.  The  Kepler  Space  Telescope  launched  in  2009  and  uses  a 
1.4-meter  ULE  glass  primary  mirror  similar  to  the  HST.  These  mirrors  have  great 
properties  for  space  telescope  missions,  as  ULE  glass  is  thermally  stable,  stiff,  and  can  be 
polished  to  optical  requirements  [22],  A  downside  to  ULE  glass  mirrors  is  mass  as  the 
area  density  (mass  /  optical  surface)  is  approximately  180  kg/m  [22],  The  current 
primary  mirrors  are  monolithic  and  scaling  those  mirrors  to  a  larger  diameter  will 
decrease  the  stiffness,  requiring  an  increase  in  support  structure  to  maintain  stiffness 
requirements  [22],  The  mass  and  stiffness  requirements  make  ULE  undesirable  for  a 
large  telescope. 

2.  Beryllium  Mirrors 

In  the  mid-1990s  NASA  began  developing  beryllium  mirror  technologies  for  the 

JWST  program.  Beryllium  mirrors  offer  several  advantages  over  ULE  glass,  including 

lighter  mass,  super  cryogenic  coefficient  of  thermal  expansion  (CTE),  and  thermal 

sensitivity.  The  thermal  properties  drove  the  decision  to  use  Beryllium  since  it 
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outperformed  ULE  over  the  telescope  temperature  operating  range  [23].  The  JWST 
beryllium  mirror  segments  have  an  areal  density  of  approximately  10  kg/m  .  The 
segments  are  also  active  and  have  center  actuators  that  compensate  for  mirror  curvature. 
In  addition  to  the  JWST  the  infrared  Spitzer  Space  Telescope  uses  beryllium  for  the  85- 
cm  primary  mirror  to  meet  required  thermal  properties. 

3.  Glass  Membrane  Facesheet  with  Active  Rigid  Support 

In  the  late  1990s  the  University  of  Arizona  developed  a  prototype  mirror  known 
as  Membrane  with  Active  Rigid  Support  (MARS)  to  demonstrate  advance  mirror 
concepts  for  the  JWST  program.  The  MARS  mirrors  were  non-replicated  active  mirrors. 
A  glass  membrane  was  ground  and  polished  from  a  glass  blank.  The  resulting  2mm  thick 
glass  membrane  was  attached  to  a  rigid  carbon  fiber  epoxy  support  structure.  The 
demonstrator  for  the  JWST  had  a  2m  diameter,  166  actuators,  and  an  areal  density  of 
13kg/m  [24],  [25].  Glass  membrane  mirror  designs  traded  the  total  number  of  actuators 
and  the  membrane  thickness.  The  actuators  corrected  surface  warping  of  the  membrane, 
but  the  discrete  actuator  locations  caused  rippling  of  the  membrane.  Adding  actuators  or 
increasing  the  thickness  of  the  membrane  could  have  fixed  this  problem. 

4.  Actuated  Hybrid  Mirrors 

Actuated  hybrid  mirrors  (AHM)  are  replicated  active  mirrors  that  are  an 
alternative  to  conventional  glass  primary  mirrors.  Developed  in  the  early  2000s  by  AOA 
Xinetics,  AHM  technology  combined  actuated  integrated  zonal  meniscus  mirrors  with  a 
nanolaminate  foils  [12],  [21],  [26],  [27].  The  integrated  zonal  meniscus  mirrors  combine 
an  optical  substrate,  a  mirror  surface,  and  embedded  actuators  in  the  support  structure 
parallel  to  the  mirrors  surface.  The  actuators  generate  bending  moments  to  control  the 
shape  of  the  mirror  surface.  The  mirror  surface  is  a  nanolaminate  foil  made  from  a 
magnetron  sputtering  process  developed  at  Lawrence  Livermore  National  Laboratory 
[27].  The  process  coats  a  polished  mandrel  surface  atom  by  atom  to  form  a  foil  that  is 
10-100  microns  thick.  The  magnetron  sputtering  process  first  deposits  a  carbon  layer  that 
serves  as  a  release  layer  from  the  mandrel.  The  process  then  applies  a  reflective  layer 
like  gold  to  form  the  optical  surface.  The  next  step  adds  alternating  crystalline  layers  and 
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amorphous  layers  of  zirconium  and  zirconium/copper  to  form  the  laminate  foil  [21].  The 
substrate  bonds  directly  to  the  foil  while  attached  to  the  mandrel. 

Figure  5  shows  the  elements  of  an  AHM.  The  substrate  is  a  rib-stiffened  silicon 
carbide  (SiC)  structure  cast  to  meet  the  required  optical  figure.  The  nanolaminate  face 
sheet  forms  on  a  precision  mandrel  and  bonds  to  the  surface  of  the  SiC  substrate.  Lead 
magnesium  niobate  (PMN)  electro  static  surface  parallel  actuators  bond  to  the  substrate 
ribs  with  an  epoxy.  Fine  control  actuators  and  backing  hardware  attaches  to  the  segment 
to  position  the  mirror.  The  result  is  a  replicable  mirror  segment  that  has  an  areal  density 
of  approximately  20  kg/m  .  AHM  mirrors  greater  than  1 -meter  diameter  achieved  a  RMS 
wavefront  error  of  42  nm  [21]. 


Nanolaminato 
0.2  kg/m2 


SiC  Substrate 
7.7  kg/m2 


Facesheet  Actuators 
0.5  kg/m2 


Fine  Control  Actuators 
0.6  kg/m2 


Ancillary  Hardware 
0.5  kg/m2 


Figure  5.  Actuated  Hybrid  Mirror  from  [21] 


5.  Carbon  Fiber  Mirrors 

Interest  in  replicated  Carbon  fiber  reinforced  polymer  (CFRP)  mirrors  began  in 
the  mid-1990s  because  the  material  properties  are  attractive  for  optical  systems.  CFRP  is 
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very  stiff,  has  a  low  CTE,  thermal  conductivity  similar  to  steel,  and  low  areal  density 
between  2-10  kg/m  depending  on  structural  configuration  [28].  The  fabrication  process 
uses  a  negative  shape  polished  mandrel  to  form  the  mirror  surface.  The  manufacturer 
layers  CFRP  prepreg  over  the  mandrel  and  after  curing  applies  the  structural  core  to  the 
back  of  the  face  sheet.  The  manufacturer  then  releases  the  mirror  from  the  mandrel  and 
applies  the  optical  coating  [14],  [29].  The  optical  performance  of  these  mirrors  depends 
on  the  surface  quality  of  the  mandrel  as  well  as  mitigating  fiber  and  core  print  through 
issues.  Figure  6  shows  the  1-meter  diameter  CFRP  mirror  at  the  Naval  Postgraduate 
School  with  an  areal  density  of  17.6  kg/m  including  the  substrate  structure.  The  residual 
wavefront  error  of  this  mirror  is  presented  later  in  Chapter  III. 


Figure  6.  1 -meter  carbon  fiber  reinforced  polymer  mirror 

CMA,  Naval  Research  Faboratory  (NRF),  and  Sandia  National  Faboratories 
(SNF)  developed  active  CFRP  mirrors  and  showed  that  astigmatism  is  inherent  to  large 
aperture  CFRP  mirrors  due  to  the  isotropic  nature  of  CFRP.  CMA,  NRF,  and  SNF 
mounted  actuators  to  the  back  of  the  CFRP  mirror  to  correct  low  order  aberrations  using 
closed  loop  feedback.  They  experimentally  demonstrated  an  actuation  method  to  change 
the  CFRP  mirror  radius  of  curvature  [30],  [31]. 

Steeves  et.  al  [15]  developed  another  CFRP  mirror  concept  called  carbon  shell 

mirrors  (CSMs).  The  mirror  substrate  is  a  thin  shell  made  from  CFRP.  The  reflective 
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layer  is  a  nano-laminate  bonded  to  the  front  of  the  CFRP  substrate.  An  active  layer  is 
bonded  to  the  back  of  the  substrate  to  provide  surface  parallel  actuation.  Figure  7  shows  a 
diagram  of  the  mirror.  CSMs  with  a  diameter  of  175  mm  have  been  prototyped  with  an 
areal  density  of  1.5-2. 7  kg/m  ,  but  do  not  include  a  support  structure.  The  current 
prototypes  do  not  meet  surface  requirements  for  visible  wavelength  applications,  but 
demonstrate  potential  for  large  correctable  space  telescope  mirrors. 


Reflective  NL  Facesheet 
CFRP  Substrate 
Ground  Electrode 
Active  Piezoelectric  Layer 
Patterned  Electrodes 
Electrode  Routing  Layer 


Figure  7.  Exploded  diagram  of  the  CSM  layers  from  [15] 


C.  ADAPTIVE  OPTICS 

Adaptive  optics  is  a  multidisciplinary  field  focused  on  improving  the  performance 
of  an  optical  signal.  The  SMT  design  and  this  research  are  founded  on  the  principles  of 
adaptive  optics.  This  section  provides  a  background  of  adaptive  optics  systems  including 
woofer-tweeter  adaptive  optics.  The  principle  components  of  an  adaptive  optics  system 
are  described  including  wavefront  sensors,  deformable  mirrors,  and  adaptive  optics 
control  techniques.  The  research  in  this  dissertation  builds  on  these  principles. 

1.  Background 

The  two  major  historical  applications  of  adaptive  optics  have  been  laser  beam 
control  and  astronomical  imaging.  Initial  work  in  adaptive  optics  began  when  Babcock 
[3]  introduced  the  concept  of  adaptive  optics  in  1953  to  compensate  for  atmospheric 
scintillation  when  making  astronomical  observations.  Later  in  1957,  Linnik  [32] 
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proposed  using  segmented  mirrors  to  compensate  for  atmospheric  turbulence  including 
the  use  of  a  reference  beacon. 

The  invention  of  the  laser  in  1960  led  to  the  use  of  adaptive  optics  in  laser  beam 
control  and  atmospheric  compensation  for  defense  purposes.  Hardy  [33]  and  Greenwood 
[34]  give  a  historical  overview  of  early  adaptive  optics  research  focused  on  defense 
applications  from  the  late  1960s  through  the  early  1990s.  One  use  of  adaptive  optics  was 
the  development  of  tactical  and  strategic  military  laser  systems.  Initial  adaptive  optics 
research  and  experiments  compensated  for  laser  propagation  through  the  atmosphere. 
Later  research  investigated  compensating  thermal  blooming  of  high-energy  lasers 
through  the  atmosphere  for  strategic  defense  applications. 

Image  compensation  adaptive  optic  applications  began  in  the  1980s  to  improve 
the  performance  of  large  ground  telescopes  used  by  the  DOD  for  space  surveillance. 
Adaptive  optics  compensated  for  atmospheric  turbulence  to  improve  telescope 
performance.  Hardy  [33]  summarizes  the  early  development  of  astronomical  telescope 
compensation  systems  including  the  use  of  laser  guide  stars.  These  techniques  were  later 
declassified  and  transferred  to  the  astronomical  telescope  community  in  the  early  1990s. 

Two  mirror  adaptive  optics  systems  were  introduced  in  the  1990s  to  compensate 
for  phase  and  amplitude  of  lasers  propagated  through  atmospheric  turbulence  and  thermal 
blooming  [35],  [36].  The  air  borne  laser  (ABL)  program  used  several  adaptive  optics 
systems  to  compensate  for  atmospheric  turbulence  for  the  high-energy  laser  weapon.  The 
ABL  program  introduced  a  woofer-tweeter  system  to  compensate  for  wavefront  errors 
caused  by  thermal  heating  of  the  optics  and  errors  from  the  laser  cavity  [37].  The 
deformable  mirror  called  a  woofer  compensates  for  low  spatial  and  temporal  frequency 
disturbances  while  the  tweeter  compensates  for  high  spatial  and  temporal  frequency 
disturbances. 

Hampton  et  al.  proposed  woofer-tweeter  systems  for  ground-based  astronomical 
telescope  atmospheric  correction  [38].  These  systems  addressed  the  adaptive  optic 
system  requirements  for  future  large  aperture  ground-based  telescopes.  Instead  of 
designing  a  single  deformable  mirror  to  meet  both  large  actuator  stroke  and  high  actuator 
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density  requirements  the  researchers  proposed  using  two  of  deformable  mirrors.  The 
woofer  deformable  mirror  has  large  stroke,  fewer  actuators,  and  correct  large  amplitude 
but  low  spatial  and  temporal  frequency  errors.  The  tweeter  deformable  mirror  has  high 
actuator  density  with  smaller  stroke  and  corrects  small  amplitude  high  spatial  and 
temporal  frequency  errors.  Gavel  and  Norton  demonstrated  woofer-tweeter  techniques 
using  the  Lick  3-meter  telescope  [39]. 

2.  Components 

The  primary  components  of  an  adaptive  optics  system  are  the  wavefront  sensor, 
wavefront  corrector,  and  control  system.  Configurations  of  these  components  depend  on 
the  application,  but  the  purpose  of  the  wavefront  sensor  is  to  provide  feedback,  the 
wavefront  corrector  compensates  the  wavefront,  and  the  control  computer  reconstructs 
the  wavefront  and/or  computes  the  control  signal  for  the  wavefront  compensator.  Figure 
8  shows  an  example  adaptive  optics  system. 


Telescope 


Figure  8.  Example  adaptive  optics  system 
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a.  Wavefront  Sensors 

Wavefront  sensors  in  adaptive  optics  systems  provide  feedback  to  the  wavefront 
corrector  and  operate  at  a  frequency  greater  than  the  disturbance  being  corrected.  Typical 
wavefront  sensors  in  adaptive  optics  systems  include  Shack-Hartmann,  curvature,  and 
phase  diversity.  A  Shack-Hartmann  wavefront  sensor  output  is  proportional  to  the 
wavefront  slope  measured  at  the  pupil  plane.  The  wavefront  can  be  reconstructed  from 
these  measurements.  A  curvature  sensor  and  phase  diversity  sensor  measure  irradiance  at 
two  locations  at  or  near  the  image  plane,  but  use  different  phase  retrieval  techniques. 

Figure  9  shows  a  one-dimensional  geometric  representation  of  a  Shack-Hartmann 
wavefront  sensor.  The  sensor  consists  of  an  array  of  lenses  in  front  of  a  camera.  Light 
passes  through  the  lenses  and  creates  small  focus  spots  on  the  camera.  The  spots  are 
proportional  to  the  local  wavefront  tilt  at  the  location  of  each  lens.  The  local  wavefront 
slopes  is  calculated  from  the  lateral  shifts,  Axy,  of  the  local  focal  point  on  the  sensor.  To 
estimate  the  phase  we  use  the  relationship  between  slope  and  phase  points  as  described 
by  Southwell  [40]. 


Figure  9.  Shack-Hartmann  wavefront  sensor,  after  [40] 
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A  curvature  sensor  uses  an  image  sensor  to  take  two  highly  out  of  focus  images  in 
two  planes  that  are  symmetric  about  the  focal  plane  [41].  The  intensity  of  the  image 
pixels  measures  the  irradiance  distribution  at  the  out  of  focus  planes.  Figure  10  shows  an 
example  of  a  curvature  sensor  described  by  Roddier  [42],  The  difference  in  irradiance 
distribution  between  the  two  planes  is  a  measure  of  the  local  wavefront  curvature.  Using 
the  sensor  data  and  solving  the  irradiance  transport  equation  with  the  proper  boundary 
conditions  estimates  the  phase  [43].  Figure  11  shows  an  example  of  the  curvature  sensor 
measurements,  difference  image  and  reconstructed  wavefront  for  a  wavefront  with 
astigmatism. 


Lens  1  Plane  1 


Plane  2 


Wavefront 


focal  length 


Figure  10.  Curvature  wavefront  sensor  example,  after  [42] 
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Figure  11.  Curvature  wavefront  sensor  example  of  wavefront  with  astigmatism. 


Phase  diversity  sensors  also  measure  intensity  to  retrieve  phase  information. 
However,  the  intensity  measurements  are  usually  asymmetric  with  a  measurement  at  the 
focal  plane  and  one  just  slightly  out  of  focus  on  the  order  of  a  single  wave.  The 
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wavefront  is  estimated  from  the  intensity  measurements  using  phase  retrieval  algorithms 
like  the  Gerchberg-Saxton  algorithm  [44]  or  modified  Gerchberg-Saxton  algorithm  [45]. 
These  algorithms  use  the  fact  that  the  Fourier  transform  relates  the  complex  amplitude 
field  at  the  pupil  and  the  image  plane.  The  image  intensity  measured  from  the  phase 
diversity  camera  provides  the  magnitude  of  the  complex  amplitude  but  not  the  phase. 
However,  by  iteratively  processing  multiple  defocused  images  one  can  correlate  the 
intensity  variations  between  the  images  and  compute  the  wavefront. 

The  wavefront  sensors  work  well  for  fine  phasing  and  surface  figure  control, 
however  the  sensors  cannot  determine  wavefront  phase  differences  greater  than  2k.  The 
27t  ambiguity  is  caused  by  alignment  piston  error  of  the  segments.  Coarse  segment 
phasing  requires  a  different  sensor  and  phase  retrieval  technique  that  can  retrieve  phase 
errors  greater  than  a  wavelength.  Dispersed  fringe  sensing  (DFS)  was  developed  to  phase 
the  JWST  segments  and  has  been  experimentally  tested  using  the  ground-based  Keck 
telescope.  DFS  can  produce  a  direct  measurement  of  the  magnitude  of  the  piston  error. 
The  technique  uses  a  grated  prism  placed  at  the  interface  of  two  segments  to  disperse 
broadband  light  from  a  point  source  onto  a  camera  to  create  a  fringe  image.  To  solve  for 
piston  a  least  squares  routine  uses  the  measured  intensity  of  the  fringe  images  to  solve  the 
parameters  of  the  fringe  equation  (2.4)  where  I  is  the  dispersed  fringe  intensity,  x  is  the 
dispersion  coordinate,  Io  is  the  mean  intensity,  y  is  the  fringe  visibility  between  0  and  1, 
and  <|>  is  the  phase  constant  [46]. 

/(x>y)  =  /o 

b.  Deformable  Mirror 

The  tip/tilt  mirror  is  one  of  the  simplest  wavefront  compensators  and  can  operate 
at  very  high  temporal  frequencies  but  only  compensates  for  tip/tilt  spatial  modes. 
Deformable  mirrors  can  correct  more  spatial  frequencies  by  adding  more  actuators,  but 
the  temporal  frequency  is  typically  reduced.  Figure  12  shows  a  plot  of  spatial  frequency 
versus  temporal  frequency  for  a  deformable  mirror.  The  spatial  frequencies  are 
represented  as  a  fraction  of  the  mirror  diameter.  The  deformable  mirror  correction 
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capabilities  are  also  dependent  on  actuator  stroke.  Stroke,  spatial  frequency,  and 
temporal  frequency  are  deformable  mirror  design  trades.  Generally  low  spatial  frequency 
correction  requires  large  amplitudes  while  higher  spatial  frequency  correction  requires 
small  amplitudes.  Therefore  most  adaptive  optic  systems  include  multiple  wavefront 
compensators  to  meet  wavefront  correction  requirements  in  terms  of  stroke,  spatial 
correction  and  temporal  correction. 


Figure  12.  Spatial  frequency  versus  temporal  frequency  of  wavefront 
compensating  mirrors,  after  [47] 


Figure  13  shows  examples  of  current  deformable  mirror  actuator  technologies  and 
geometries.  The  two  main  actuator  configurations  are  surface  normal  and  surface 
parallel.  Deformable  mirrors  with  surface  normal  actuators  create  a  local  influence  on 
the  mirror  surface  while  surface  parallel  actuators  create  both  a  local  influence  and  a 
global  influence.  In  both  cases  the  mirror  response  is  designed  to  be  linear,  and  the 
actuator  responses  obey  the  principle  of  linear  superposition.  This  means  that  the  entire 
mirror  response  to  a  set  of  actuator  commands  is  a  linear  combination  of  the  response  of 
the  individual  actuators  being  commanded. 
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Some  deformable  mirror  technologies  combine  different  types  of  actuators 
geometries.  For  example,  a  surface  parallel  mirror  may  have  global  control  actuators  like 
in  the  top  right  mirror  in  Figure  13  and  another  set  of  local  actuators  nested  within  the 
substrate  ribs.  This  configuration  increases  the  actuator  density  and  allows  for  the 
correction  of  more  spatial  frequencies,  but  the  actuators  experience  coupling  and  the 
assumption  of  linear  superposition  does  not  hold. 

Deformable  mirror  actuators  are  typically  piezoelectric  or  electrostrictive  and 
generate  a  displacement  producing  a  constant  strain  on  the  mirror  structure.  Lead 
zirconate  titanate  (PZT)  piezoelectric  actuators  produce  a  strain  linearly  proportional  to 
the  applied  electric  field.  Lead  magnesium  niobate  (PMN)  electrorestrictive  actuators 
produce  a  strain  proportional  to  the  square  of  the  applied  electric  field  [48]. 
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Figure  13.  Deformable  mirror  technologies,  from  [49] 


3.  Adaptive  Optics  Control 

Most  adaptive  optics  control  techniques  exploit  the  linear  behavior  of  deformable 
mirrors.  In  addition  to  the  assumption  of  linear  response  we  assume  the  deformable 
mirrors  dynamics  are  much  faster  than  the  bandwidth  of  the  feedback  sensors.  This 
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allows  the  deformable  mirror  control  system  to  use  a  static  controller  built  from  an 
influence  matrix.  An  influence  matrix  converts  deformable  mirror  actuator  inputs  to  a 
deformable  mirror  response,  usually  the  wavefront  sensor  output.  The  influence  matrix 
contains  column  vectors  called  influence  functions,  equation  (2.6).  To  form  the  influence 
function  the  mirror  is  set  to  a  biased  position  and  each  actuator  is  individually  actuated  in 
the  positive  and  negative  direction.  Equation  (2.5)  shows  the  positive  and  negative 
mirror  responses  are  averaged,  and  normalized  over  the  control  input  voltage,  v,  for  each 
actuator,  a,.  Equation  (2.7)  relates  the  input  control  vector  u  to  the  mirror  output  y  using 
the  influence  matrix  T.  The  reconstruction  matrix  relates  the  desired  output  to  a  set  of 
control  signals  by  taking  the  Moore-Penrose  pseudo-inverse  T,  equation  (2.8). 
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Most  adaptive  optics  systems  use  a  feedback  integral  controller.  In  this  control 
approach  the  control  signal  is  proportional  to  the  integral  of  the  error.  Equation  (2.10) 
implements  the  discrete  time  integral  control  using  the  reconstruction  matrix  R,  a  gain  g, 
and  error  signal  defined  in  equation  (2.9)  where  y  is  the  measured  wavefront  and  yo  is  the 
desired  wavefront. 

ek  =  y-y0  (2.9) 

uk+i  =  gRek+uk  (2-10) 

D.  SEGMENTED  TELESCOPES 

The  telescope  used  in  this  research  is  a  segmented  telescope.  The  SMT  was  one 
of  the  first  visible  large  aperture  segmented  telescopes  where  the  control  system  was 
included  in  the  optical  prescription.  This  section  provides  background  on  the  reason  for 
using  segmented  telescopes,  provides  a  review  of  segmented  telescope  research,  and 
describes  the  SMT. 
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1.  Background 

Manufacturers  and  researchers  can  fabricate  larger  and  lighter  monolithic  mirrors 
than  the  HST.  However,  scaling  existing  traditional  monolithic  mirrors  introduces  mirror 
stiffness  problems.  The  mirror  must  be  stiff  enough  to  allow  manufacturing,  maintain  the 
proper  shape,  and  survive  launch.  Increasing  the  mirror  diameter  without  increasing  the 
thickness  drastically  reduces  the  mirror’s  stiffness  requiring  the  addition  of  stiffening 
support  structure.  Large  aperture  telescopes  require  a  segmented  mirror  to  fit  in  the 
launch  vehicle  fairing  [22], 

The  transition  of  space-based  telescopes  from  monolithic  primary  mirrors  to 
segmented  mirrors  is  partially  based  on  the  success  of  ground  based  astronomical 
segmented  telescopes.  A  ground-based  segmented  telescope  was  proposed  in  the  late 
1970s  in  order  to  build  cost  effective  10-meter  diameter  telescopes.  The  result  was  the 
first  of  two  Keck  telescopes  completed  in  1993  using  36  segments.  The  Keck  telescope 
actively  aligns  the  segments  and  adjusts  the  segment  surface  figure  to  ensure  the 
segments  act  as  a  single  surface.  There  are  several  other  large  ground-based  segmented 
telescopes  that  use  active  optics  including  the  Hobby-Eberly  Telescope,  Gran  Telescop io 
Canaritas,  the  Southern  African  Large  Telescope  (SALT),  and  the  Large  Sky  Area  Multi- 
Object  Fibre  Spectroscopic  Telescope  (LAMOST).  Figure  14  shows  a  comparison  of 
telescope  primary  mirrors. 

Research  in  deployable  space  telescopes  has  been  ongoing  since  the  1970s. 
Research  focused  on  the  objective  of  phasing  the  telescope  segments.  This  required  work 
in  many  fields  including  structures,  wavefront  sensing,  and  controls.  Lane  et  al.  [9] 
provides  an  overview  of  past  research  efforts.  In  the  1980s  and  1990s  the  Lockheed 
Martin  adaptive  large  aperture  technology  (ALOT)  program  demonstrated  segmented 
mirror  alignment  and  wavefront  sensing  and  control  technologies.  Lockheed  Martin  also 
analyzed  structural  behavior  of  a  segmented  primary  mirror  on  the  advanced  structures 
and  controls  integrated  experiment  (ASCIE).  This  research  informed  initial  work  on  the 
JWST  in  the  mid-1990s. 
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Technology  demonstrations  and  risk  reduction  efforts  continued  into  the  2000s. 
Massachusetts  Institute  of  Technology  (MIT)  developed  a  sparse  aperture  interferometer 
test  bed  called  Adaptive  Reconnaissance  Golay-3  Optical  Satellite  (ARGOS).  ARGOS 
identified  the  multidisciplinary  design  challenges  of  integrating  building  interferometric 
array  and  integrating  in  a  spacecraft  simulator  [50].  Air  Force  Research  Laboratory 
(AFRL)  also  built  the  Deployable  Optical  Telescope  (DOT)  to  demonstrate  the 
deployment  and  calibration  of  a  three-segment  telescope. 

More  recent  technology  developments  include  the  SMT  and  NASA  UV/EO  space 
telescope  work  like  the  Advanced  Technology  Large  Aperture  Space  Telescope 
(ATLAST).  ATLAST  is  an  8  to  16-meter  ultraviolet  optical  infrared  space  observatory 
[51].  NASA  has  also  proposed  concepts  to  robotically  assemble  a  large  segmented 
telescope  on  orbit.  This  resulted  in  a  risk  reduction  program  for  ATLAST  called  Optical 
Testbed  and  Integration  of  ISS  experiment  (OpTIIX).  The  OpTIIX  program  intended  to 
demonstrate  on  orbit  telescope  assembly  technologies  and  end-to-end  wavefront  sensing 
and  control  of  an  optical  segmented  telescope  [52], 
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Figure  14.  Comparison  of  telescope  primary  mirrors,  from  [53] 


2.  SMT  Description 

The  3 -meter  6  segment  SMT  requires  the  segment  phasing  in  order  for  the  SMT 
to  operate  as  a  single  monolithic  mirror.  This  means  that  the  surface  of  each  mirror 
segment  must  align  within  close  proximity  of  the  design  prescription.  Generally  a  peak  to 
valley  wavefront  error  of  A/4  is  considered  phased.  The  segment  surface  RMS  wavefront 
error  should  also  have  a  RMS  wavefront  error  less  than  A/15.  Phasing  the  SMT  primary 
mirror  requires  the  SMT  actuators,  wavefront  sensors  and  a  control  system. 

The  3-meter  diameter  primary  mirror  has  six  hexagonal  1 -meter  diameter  AHM 
segments  with  156  PMN  surface  parallel  face-sheet  actuators  (FSAs)  each.  These 
actuators  correct  the  surface  of  the  segment.  Each  segment  is  aligned  using  six  coarse 
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control  actuators  (CCAs)  and  three  fine  control  actuators  (FCAs).  Each  FCA  is 
connected  to  a  bipod  and  each  leg  of  the  bipod  is  a  CCA.  The  CCA  has  a  stroke  length  of 
+/-  2.8mm  with  a  step  resolution  of  1.5 jum.  The  six  CCAs  per  a  segment  act  as  a  hexapod 
and  are  used  for  coarse  alignment  and  phasing.  The  FCA  actuators  are  PMN  actuators 
like  the  FSAs.  They  have  a  step  resolution  of  less  than  1  nm  and  a  total  stroke  of  +/-  5 jum. 
The  three  segment  FSAs  control  segment  piston,  tip,  and  tilt.  Figure  15  shows  the 
position  of  the  FSAs  and  CCAs  on  the  segment  back  structure. 


CCA2 


FCA1 


Figure  15.  SMT  FCA,  CCA,  and  GS  locations 


Figure  16  shows  the  SMT  optical  layout.  The  SMT  has  a  three-mirror  anastigmat 
optical  design.  The  primary,  secondary,  and  tertiary  mirrors  all  have  curvature  in  order  to 
remove  spherical  aberrations,  coma,  and  astigmatism.  The  primary  and  secondary  mirror 
work  similar  to  a  Cassegrain  design  and  the  tertiary  mirror  enables  a  large  system  focal 
length  of  60  meters  giving  a  f#  of  20.  The  magnification  of  the  tertiary  relay  mirror 
enables  the  long  system  focal  length. 
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Figure  16.  SMT  optical  layout,  from  [54] 


The  satellite  field  of  regard  (FOR)  is  the  target  area  that  the  telescope  can  see 
based  on  the  telescope  position,  while  the  field  of  view  (FOV)  is  the  area  that  the  sensor 
can  image  within  the  FOR.  The  SMT  FOR  is  0.5  degrees  circular  and  the  design 
provides  a  full  0.5  x  0.5  degree  rectangular  FOR.  The  sensor  FOV  is  2x2  arc  minutes  and 
can  only  image  a  small  portion  of  the  FOR.  The  SMT  includes  a  field  steering  mirror 
(FSM)  at  the  exit  pupil  to  steer  the  FOV  over  the  FOR.  The  FSM  also  removes 
mechanical  jitter  from  the  optical  assembly.  Figure  17  shows  an  example  of  the  SMT 
FOV  within  the  FOR. 
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Field  of  View 


Figure  17.  SMT  Field  of  Regard  and  Field  of  View 


The  SMT  wavefront  sensing  system  includes  edge  sensors,  jitter  sensor,  a  phase 
diversity  wavefront  sensor,  and  a  Shack-Hartmann  wavefront  sensor.  The  phase  capture 
range  of  the  sensors  varies.  The  Shack-Hartmann  and  phase  diversity  wavefront  sensors 
suffer  from  a  2k  ambiguity,  but  provide  feedback  for  fine  phase  control  and  surface 
figure  control.  The  edge  sensors  are  not  susceptible  to  the  27t  ambiguities  and  provide 
feedback  for  coarse  alignment  and  phasing  of  the  segments.  The  coarse  phasing  control 
approach  uses  both  the  edge  sensors  and  the  phase  diversity  sensor  data  to  remove  edge 
sensor  bias  terms.  The  SMT  edge  sensors  do  not  meet  design  requirements  and  currently 
prohibit  segment  phasing  using  the  existing  control  architecture. 

More  recent  segmented  telescope  research  introduces  wavefront  sensing  and 
control  techniques  including  a  laser  truss  metrology  system  and  dispersed  fringe  sensing 
[46],  [55],  [56].  NASA  JPL  developed  laser  truss  metrology  techniques  to  perform 
segment  phasing  in  the  ATLAST  and  OpTIIX  programs.  Figure  18  shows  the  laser  truss 
metrology  concept.  A  laser  truss  system  measures  the  full  optical  state  of  the  telescope  at 
high  bandwidth  and  offers  advantages  over  edge  sensors  used  in  the  SMT.  The  system 
uses  laser  light  to  measure  the  position  of  the  primary,  secondary  and  tertiary  mirrors  in 
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relation  to  each  other.  According  to  NASA  JPL  the  laser  truss  system  provides  better 
observably  of  segment  rigid  body  errors,  can  operate  when  there  is  no  guide  star, 
provides  measurements  that  support  alignment  across  the  full  field  of  view,  and  can 
estimate  all  the  required  degrees  of  freedom  [55]. 


Figure  18.  Space  telescope  laser  truss  metrology  concept,  from  [55] 


3.  SMT  State  Space  Model  Background 

Analytical  modeling  of  segmented  telescopes  is  necessary  to  design  the  wavefront 
sensing  and  control  system.  The  control  design  in  this  dissertation  is  based  on  analytical 
models  developed  by  Yingling  [54],  Horn  and  Bampton  [57].  They  developed  integrated 
optical-mechanical  SMT  models  resulting  in  a  state  space  plant  model  for  the  SMT.  The 
state  space  model  represents  the  structural  and  optical  response  of  the  SMT.  The  state 
space  model  inputs  include  all  the  actuator  inputs  and  seismic  disturbances.  The  model 
outputs  the  optical  path  difference  (OPD)  of  light  at  the  focal  plane  array  pupil  and  the 
displacement  of  the  gap  sensors.  The  OPD  output  is  in  the  form  of  a  282  x  282  matrix. 
The  sensor  outputs  are  converted  from  the  OPD  output.  Table  2  shows  the  model  inputs 
and  sensor  outputs.  Figure  19  shows  the  state  space  model  layout. 
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Table  1.  SMT  state  space  model  inputs  and  outputs 


Control  Inputs 

Source 

Number 

per  segment 

Face  Sheet  Actuators 

936 

156 

Fine  Control  Actuators 

18 

3 

Coarse  Control  Actuators 

36 

6 

Fast  Steering  Mirror 

4 

Seismic 

3 

Total 

997 

Outputs 

Source 

Number 

per  segment 

Shack-Hartmann  Sensor 

732 

122 

Phase  Diversity  Sensor 

18 

3 

Gap  Sensor 

18 

3 

Jitter  Sensor 

2 

Modal  States 

166 

Total 

936 

FSA  displacement 
(m)  156x6 


FCA  displacement 
(m)  5x6 


CCA  displacement 
(m)  6x6 

FSM  Reacted 
Torque  (N  m)  2 


Seismic  (m/sA2)  3 


FSM  Command 
(rad)  2 


Figure  19.  SMT  Opto-mechanical  state  space  model  flow  diagram,  from  [57] 
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4.  SMT  Control  Approach 

The  SMT  wavefront  sensing  and  control  system  is  required  to  move  the  segments 
from  their  post  deployment  un-phased  positions  to  their  operational  diffraction  limited 
position.  The  initialization  process  begins  after  telescope  deployment  and  improves  the 
wavefront  error  from  several  millimeters  to  below  the  diffraction  limit.  The  three- 
wavefront  sensing  and  control  modes  are  coarse  segment  alignment,  coarse  phasing,  and 
fine  phasing.  The  coarse  segment  alignment  mode  identifies  the  segment  on  the  camera 
focal  plane  and  moves  the  segment  to  center  of  the  field  reducing  the  alignment  error 
from  millimeters  to  micrometers.  The  coarse  phasing  mode  co-aligns  the  segments 
resulting  in  a  wavefront  error  of  less  than  a  wavelength.  The  fine  phasing  mode  uses  the 
phase  diversity  sensor  and  FCAs  to  adjust  segment  piston,  tip/tilt,  and  uses  the  Shack 
Hartmann  wavefront  sensor  to  reduce  the  wavefront  error  of  the  segment  using  the  FSAs. 
Table  2  shows  the  developer’s  SMT  control  loops,  the  associated  inputs  and  outputs,  and 
the  phasing  mode  they  support.  The  SMD  program  did  not  demonstrate  a  real-time  high- 
bandwidth  closed  loop  control  system  [58]. 


Table  2.  SMT  control  loops 


Control  Mode 

Controller 

Sample  Rate 
(Hz) 

Inputs 

Outputs 

Coarse  Alignment 

Bipod  Offload 

0.1 

Gap  Sensor 

CCAs 

Coarse  Phasing 

Segment  Rigid  Body 

1000 

Gap  Sensor,  Full 
Aperture  Phase 
Diversity 

FCAs 

Fine  Phasing 

Segment  Shape 

200 

Shack-Hartmann 

FSAs 

Phase  Diversity 

1 

Phase  Diversity 
Sensor 

FCAs 

The  initialization  process  brings  the  segments  into  phase  and  reduces  the  primary 
mirror  surface  wavefront  error.  The  telescope  points  toward  a  coherent  light  source 
reference  like  a  star.  The  initialization  process  also  establishes  the  relative  set  point 
conditions  for  the  edge  sensors.  The  maintenance  control  takes  place  during  mission 
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operations  and  uses  the  CCAs,  FCAs,  edge  sensors,  and  phase  diversity  sensor  to 
maintain  the  segment  phasing  using  a  closed  loop  controller.  The  edge  sensors  provide 
high  bandwidth  relative  positions  of  the  segments  while  the  phase  diversity  sensor 
provide  low  bandwidth  absolute  reference  to  correct  for  edge  sensor  bias  and  drift  using  a 
Kalman  filter.  The  FSAs  are  designed  to  hold  the  segment  shape  for  extended  periods  of 
time  allowing  the  segment  shape  to  be  updated  periodically  during  calibration 
opportunities. 

The  SMT  control  system  also  compensates  for  jitter  using  a  fast  steering  mirror 
and  a  jitter  beacon.  The  fast  steering  mirror  is  also  the  field  steering  mirror  (FSM).  The 
FSM  is  placed  in  a  common  path  with  the  imaging  sensor  and  compensates  for  wavefront 
tip  and  tilt.  The  FSM  control  loop  is  operates  at  1kHz. 

E.  SMT  RESEARCH  AT  NAVAL  POSTGRADUATE  SCHOOL 

Segmented  mirror  telescope  research  at  the  Naval  Postgraduate  School  has 
included  both  analytical  and  experimental  work  by  students,  staff  researchers,  and 
external  organizations.  The  development  of  a  segmented  mirror  space  telescope  is  a 
multi-disciplinary  effort.  Past  areas  of  research  include  controls,  structures,  vibrational 
analysis,  segmented  telescope  design,  and  wavefront  reconstruction.  Key  works  that 
inform  this  dissertation  are  surveyed  below. 

Past  SMT  control  research  was  conducted  analytically  with  some  experiments 
using  test  beds  other  than  the  SMT.  Burtz  [59]  applied  a  multi-input  multi-output 
(MIMO)  FL  controller  to  a  flexible  segmented  space  telescope  dynamic  model  that 
included  model  reduction  techniques  to  reduce  computer  simulation  complexity.  The 
control  technique  was  experimentally  tested  using  an  adaptive  optics  test  bed.  The  work 
by  Burtz  was  continued  by  Looysen  [60]  who  combined  Hx  control  techniques  with  low 
pass  filters  to  improve  plant  stability. 

Several  researchers  have  addressed  the  computational  burden  of  a  large  MIMO 
control  problem  and  complex  structural  model.  Nagashima  and  Agrawal  [61]  reduced 
the  model  dimensions  by  decoupling  the  dynamic  and  static  portions  of  the  telescope 

plant  using  a  sensitivity  decoupling  method.  They  applied  an  controller  to  the  new 
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plant  and  showed  improved  performance  compared  to  the  uncoupled  plant.  Axtell  [62] 
developed  a  wavelet  reconstruction  technique  to  improve  the  computational  efficiency  for 
systems  with  many  sensor  outputs  like  the  SMT.  Dras,  Jennings,  and  Cobb’s  [63] 
developed  a  reduced  order  model  for  the  SMT.  They  experimentally  validated  the  model 
using  SMT  vibration  data  enabling  faster  model  tuning  and  analysis. 

In  addition  to  controls  and  wavefront  reconstruction  researchers  developed  opto¬ 
mechanical  design  techniques  and  performed  detailed  structure  and  vibration  analysis. 
Yingling  [54]  developed  an  integrated  opto-mechanical  design  method  resulting  in  a 
model  that  was  experimentally  validated  using  the  SMT.  Yingling’s  model  is  the  basis 
for  the  models  developed  in  this  research.  Yingling  also  added  tuned  mass  dampers  to 
the  SMT  to  reduce  structural  disturbances  allowing  improved  experimental  controls 
research  in  this  dissertation. 

Most  recently  Watson  and  Bagnasco’s  [64],  [65]  added  pupil  relay  optics  and  a 
deformable  mirror  the  SMT  experimental  setup.  This  work  provided  a  starting  point  for 
the  experimental  woofer-tweeter  research  in  this  dissertation. 

This  dissertation  builds  on  previous  research  in  segmented  mirror  telescopes, 
adaptive  optics,  wavefront  sensing,  and  control.  This  dissertation  responds  to  the  residual 
surface  error  of  lightweight  active  mirrors  by  offering  control  techniques  to  compensate 
for  the  residual  surface  errors  using  an  additional  deformable  mirror  in  a  woofer-tweeter 
configuration.  The  control  approaches  are  specific  to  a  space  telescope  application,  and 
the  simulation  and  experimental  results  are  demonstrated  using  an  AHM. 
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III.  SPACE  TELESCOPE  WOOFER-TWEETER  MODEL 


In  preparation  for  the  development  of  space  telescope  woofer-tweeter  control 
techniques  this  chapter  presents  the  segmented  telescope,  deformable  mirror,  and  woofer- 
tweeter  analytical  models.  This  requires  understanding  the  SMT  woofer-tweeter  design 
assumptions  and  an  understanding  of  the  SMT  design.  The  SMT  design  drives  several 
trades  when  determining  the  size  and  placement  of  the  deformable  mirror  in  the  optical 
path.  The  chapter  provides  a  design  recommendation  based  on  the  assumptions  and 
design  drives  which  informs  the  development  of  the  woofer-tweeter  model.  The  chapter 
begins  by  presenting  design  assumptions  and  design  drivers  for  the  optical  configuration. 
Next  the  chapter  describes  the  placement  of  the  deformable  mirror.  Then  the  chapter 
describes  analytical  models  of  the  surface  error,  the  SMT,  the  deformable  mirror,  and  the 
woofer-tweeter  system. 

A.  DESIGN  ASSUMPTIONS 

Three  key  assumptions  are  made  when  designing  the  SMT  woofer-tweeter 
system.  The  first  assumption  is  that  wavefront  disturbances  are  low  frequency  allowing 
the  two  mirrors  to  operate  at  the  same  bandwidth.  The  second  assumption  is  that  the 
optical  layout  only  requires  one  wavefront  sensing  system.  The  third  assumption  is  that 
the  SMT  segment  edges  are  already  phased  prior  to  woofer-tweeter  fine  phasing.  These 
assumptions  are  discussed  further  below. 

1.  Disturbance  Bandwidth 

The  purpose  of  the  proposed  space  telescope  woofer-tweeter  system  is  to  correct 
primary  mirror  static  residual  errors,  and  low  bandwidth  disturbances.  O’meara,  Swigert, 
and  Brown  noted  that  surface  figure  control  of  an  orbiting  telescope  is  different  from  a 
ground-based  telescope  based  on  the  frequency  of  the  error  sources  [17].  This  allows  the 
two  mirrors  to  operate  at  the  same  bandwidth  unlike  a  woofer-tweeter  atmospheric 
compensation  system  where  the  mirrors  correct  disturbances  at  different  temporal 
frequencies  requiring  the  woofer  and  tweeter  to  operate  at  different  bandwidths  [66]— 
[68], 
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2.  Wavefront  Sensing  Approach 

In  order  to  simplify  design  and  controls  the  wavefront  sensing  system  measures 
the  combined  effects  of  the  primary  segmented  mirror  and  the  tweeter  deformable  mirror. 
The  wavefront  sensing  system  assumes  that  the  incoming  light  is  coherent  and  that  the 
only  source  of  wavefront  error  is  from  the  primary  mirror  and  deformable  mirror.  This 
single  output  approach  allows  a  relationship  between  the  woofer-tweeter  inputs  and  the 
measured  output. 

3.  Segment  Phasing  Assumption 

The  proposed  woofer-tweeter  system  improves  the  telescope  fine  phasing  by 
improving  the  surface  figure  error  of  individual  primary  mirror  segments.  This  requires 
that  the  individual  segments  already  be  phased  relative  to  each  other.  In  order  to  preserve 
actuator  stroke  and  prevent  actuator  saturation  the  woofer-tweeter  control  system  does 
not  compensate  for  piston,  tip,  and  tilt.  Any  residual  segment  piston,  tip  or  tilt  error  is 
removed  using  the  segment  FCAs. 

B.  DESIGN  DRIVERS 

The  optical  design  of  the  SMT  drives  the  proposed  woofer-tweeter  design.  This 
section  discusses  two  key  design  considerations  when  adding  a  deformable  mirror  to  the 
optical  telescope.  The  first  is  the  fact  that  the  SMT  has  a  steerable  FOV  and  the  second 
consideration  is  the  impact  of  field  angle  magnification  based  on  the  size  of  the  exit  pupil 
and  deformable  mirror  size. 

1.  SMT  Field  Steering 

The  SMT  is  a  steerable  field  of  view  system  where  the  telescope  has  a  wide  FOR 
and  the  sensor  has  a  smaller  FOV.  The  telescope  uses  the  FSM  located  at  the  exit  pupil  to 
steer  the  FOV  within  the  telescope  FOR.  In  order  for  the  telescope  to  perform  well  at  all 
FOV  locations  the  wavefront  error  due  to  the  primary  mirror  must  be  corrected  at  all  field 
angles.  This  may  require  the  deformable  mirror  to  update  the  correction  based  on  the 
telescope  FOV  setting. 
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The  logical  place  to  make  this  correction  is  at  the  exit  pupil  because  all  the  field 
angles  are  present  and  it  is  conjugate  to  the  primary  mirror.  However,  the  SMT  exit  pupil 
location  is  already  occupied  by  the  FSM.  A  design  trade  is  discussed  in  a  later  section  to 
determine  the  best  place  to  locate  the  deformable  mirror. 

2.  Field  Angle  Magnification 

The  next  consideration  is  field  angle  magnification.  The  current  SMT  exit  pupil  is 
de-magnified  to  a  diameter  is  150mm,  or  20x.  This  has  the  effect  of  magnifying  the  field 
angles  due  to  the  Lagrange  invariant.  The  Lagrange  invariant  relates  the  ray  height  and 
angle  for  all  rays  as  they  propagate  through  an  optical  system,  Equation  (3.1).  The 
variable  DPM  is  the  diameter  of  the  primary  mirror  and  the  variable  6pm  is  the  field  angle 
at  the  primary  mirror.  The  variable  Ddm  is  the  diameter  of  the  deformable  mirror  and  the 
variable  6dm  is  the  field  angle  at  the  deformable  mirror.  If  the  pupil  plane  is  de-magnified 
the  field  angles  magnify  and  the  surface  error  correction  becomes  more  difficult  at  the 
exit  pupil.  Therefore  the  field  angle  magnification  must  be  taken  into  account  when 
selecting  the  size  of  the  deformable  mirror  and  size  of  the  exit  pupil.  Equation  (3.2) 
shows  the  relationship  between  the  field  angle  at  the  exit  pupil,  and  the  field  angle  at  the 
primary  mirror  in  terms  of  angular  magnification  where  Ma  is  the  magnification  factor. 

DPMePM=DDMeDM  (3-1) 

9DM=MaGPM  (3.2) 

McComas  [69]— [71]  studied  the  field  angle  magnification  effect  and  proposed  a 
conceptual  configurable  adaptive  optics  system  that  could  correct  for  on-orbit  aberrations 
and  field  of  view  wavefront  errors.  He  reviewed  the  narrow  field  of  view  theory  and  the 
limitations  of  narrow  field  of  view  theory  using  a  wide  FOV  system.  Equation  (3.3) 
describes  the  residual  wavefront  error  after  the  deformable  mirror  for  a  narrow  FOV 
system.  The  variable  mR(x,y,t)  is  the  residual  wavefront  error,  a>(x,y;t)  is  the  wavefront 
error  from  the  primary  mirror  incident  on  the  deformable  mirror  at  angle  90 ,  and 
dm(x,y;t)  is  the  deformable  mirror  correction.  This  assumes  the  wavefront  entering  the 
telescope  is  perfect,  only  the  primary  mirror  contributes  error,  and  the  primary  and 
deformable  mirror  are  conjugates. 
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< oR(x,y\t)  =  (o(x,y\t) — —cos(00)dm(x,y,t )  (3.3) 

McComas  reformulated  the  above  equation  in  terms  of  optical  path  length,  </>( x,y) , 
and  includes  both  field  angle  at  the  primary  mirror  and  angular  magnification  due  to  the 
Lagrange  invariant.  The  OPD  is  shown  in  equation  (3.4)  where  dPM(x,y )  is  the  surface 
height  of  the  primary  mirror,  dDM  (x,y)  is  the  surface  height  of  the  deformable  mirror,  and 
0PM  is  the  field  angle  at  the  primary  mirror.  The  result  of  equation  (3.4)  is  that  a 
deformable  mirror  may  not  be  able  to  correct  all  the  field  angles  because  as  the  term 
Ma0PM  approaches  n/2  the  deformable  mirror  cannot  compensate. 

(/)R  —  2dPM  (x,y)cos(dPM )  —  2dDM  (x,_y)cos(A/a@PM )  (3.4) 

The  SMT  has  a  FOR  of  0.5  degrees  and  an  angular  magnification  of  20.  Figure 
20  shows  that  the  relative  path  length  error  at  varying  field  angles  assuming  a  primary 
mirror  displacement  of  one  wave.  The  path  length  error  is  relatively  small  because  the 
deformable  mirror  is  the  same  size  as  the  exit  pupil  and  both  the  FOR  and  magnification 
are  small. 

However,  if  a  smaller  diameter  deformable  mirror  is  used  the  exit  pupil  requires 
de-magnification.  For  example,  a  commercially  available  25mm  diameter  MEMs 
deformable  mirror  requires  us  to  de-magnify  exit  pupil  diameter  by  an  additional  factor 
of  six  resulting  in  a  field  angle  magnification  of  120  at  the  deformable  mirror.  Figure  21 
shows  the  increased  angular  magnification  substantially  increases  the  off-axis  path  length 
error.  The  deformable  mirror  can  still  compensate  for  the  primary  mirror  at  off  axis  FOV 
locations,  but  requires  more  actuator  stroke. 

The  two  angular  magnification  examples  show  how  the  Lagrange  invariant  is  a 
design  driver  when  choosing  a  deformable  mirror.  In  both  examples  the  term  MadPM  does 
not  approach  n  12  over  the  FOR  and  the  deformable  mirror  can  compensate  for  field  angle 
magnification  at  the  edge  of  the  FOR.  However  the  mirror  stroke  must  be  sized  for  the 
anticipated  primary  mirror  surface  error  and  the  field  angle  magnification  at  off  axis  FOV 
locations. 
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The  Lagrange  invariant  impacts  the  deformable  mirror  wavefront  control  system. 
In  the  20x  magnification  example  the  optical  path  length  contribution  due  to  the 
Lagrange  invariant  can  be  ignored  because  the  error  contribution  is  minimal.  However, 
the  deformable  mirror  needs  to  update  correction  to  account  for  the  Lagrange  invariant  at 
off  axis  FOV  locations  in  the  120x  magnification  case.  This  is  accomplished  by 
including  the  field  angle  in  the  deformable  mirror  controller.  Figure  20  shows  the 
relative  path  length  error  for  different  field  angles  given  a  field  angle  magnification  of  20. 
The  contours  represent  the  additional  wavefront  error  at  the  deformable  mirror  for  a 
single  wave  of  error  at  the  primary  mirror.  Figure  21  shows  the  relative  path  length  error 
for  different  field  angles  given  a  field  angle  magnification  of  120.  The  yellow  box 
represents  the  field  of  view  and  the  red  box  represents  the  field  of  regard. 

Field  angle  magnification  at  the  DM  also  complicates  the  design  of  the  optical 
system  after  the  deformable  mirror.  The  optical  components  between  the  exit  pupil  and 
image  plane  must  handle  the  magnified  field  angles.  Optical  designs  become  more 
complex  and  require  additional  components  to  remove  large  field  angle  aberrations. 
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Figure  20.  Relative  path  length  error  (waves)  for  20x  angular  magnification 
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Figure  21.  Relative  path  length  error  (waves)  for  120x  angular  magnification 
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C.  PLACEMENT  OF  DEFORMABLE  MIRROR 

This  section  presents  two  possible  configurations  for  the  addition  of  a  deformable 
mirror  to  the  SMT.  Option  one  replaces  the  FSM  with  a  deformable  mirror  on  a  gimbal. 
Option  two  adds  an  optical  relay  and  uses  a  smaller  diameter  deformable  mirror.  The 
first  option  is  selected  for  this  research  because  the  effects  of  field  angle  magnification 
are  relatively  small.  Both  configurations  are  presented  below  to  provide  context  for  the 
selected  option,  and  because  the  second  option  can  be  implemented  in  future  work  using 
a  commercial  deformable  mirror. 

The  natural  placement  of  the  deformable  mirror  in  the  SMT  optical  design  is  at 
the  exit  pupil  since  the  exit  pupil  is  conjugate  with  the  entrance  pupil.  However,  the  FSM 
currently  resides  at  the  exit  pupil  requiring  the  optical  design  to  be  modified.  This 
position  requires  the  deformable  mirror  size  to  match  the  exit  pupil  diameter,  150  mm. 
The  FSM  functionality  can  be  maintained  by  placing  the  DM  on  a  gimbal  at  the  exit  pupil 
location.  Figure  16  shows  the  optical  layout  of  the  SMT  and  Figure  22  shows  the 
proposed  DM  location. 
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Figure  22.  SMT  optical  layout  replacing  the  FSM  with  a  deformable  mirror  on 

a  gimbal,  from  [72] 


The  second  option  relays  and  reduces  the  size  of  the  exit  pupil  to  accommodate  a 
small  deformable  mirror.  This  option  reduces  the  changes  to  the  SMT  optical  layout  and 
adds  a  pupil  relay  and  deformable  mirror  after  the  FSM  and  after  the  last  fold  mirror,  M6. 
This  approach  can  match  the  exit  pupil  size  to  a  small  commercially  available  deformable 
mirror  and  allow  use  of  the  existing  wavefront  sensing  system  and  science  camera. 
Although  the  increased  exit  pupil  de-magnification  increases  the  field  angles  at  the 
deformable  mirror  the  proposed  option  can  demonstrate  a  proof  of  concept  using  the 
SMT.  Figure  23  shows  the  proposed  optical  layout  and  Figure  24  shows  the  details  of  the 
optical  relay.  The  biggest  challenge  to  this  option  is  incorporating  the  optical  relay  into 
the  back  optics  of  the  SMT. 
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Figure  23.  SMT  with  deformable  mirror  optical  relay  and  pupil  size 
reduced,  from  [72] 
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Figure  24.  Layout  of  SMT  pupil  optical  relay  with  deformable 
mirror,  from  [72] 
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D.  RESIDUAL  WAVEFRONT  ERROR  MODELING 


This  section  presents  the  residual  error  of  a  single  SMT  AHM  segment  used 
throughout  this  research.  To  demonstrate  the  breadth  of  application  the  residual  error  of 
the  CFRP  mirror  mentioned  in  the  background  is  also  presented.  These  errors  are 
representative  of  the  residual  errors  found  in  large  lightweight  primary  optics  and  include 
print-through,  fixture  error,  and  low  order  aberrations.  The  residual  errors  identified  in 
these  mirrors  are  used  to  create  a  generalized  residual  error  model  that  can  be  used  to 
simulate  hypothetical  primary  mirror  residual  surface  errors. 

The  residual  wavefront  error  is  represented  as  the  optical  path  difference  (OPD) 
from  a  flat  wavefront.  Assuming  the  light  entering  the  telescope  is  coherent,  the  OPD 
corresponds  to  the  phase  shift  of  the  wavefront  as  it  passes  through  the  imperfect  optics 
of  the  telescope.  The  telescope  residual  wavefront  error  is  measured  as  a  two- 
dimensional  OPD  at  the  exit  pupil. 

1.  SMT  AHM  Residual  Error 

The  SMT  AHM  segment  provides  a  reference  for  expected  residual  surface  error. 
Figure  25  shows  the  measured  OPD  of  SMT  segment  3.  The  RMS  wavefront  error  is 
approximately  0.32  waves  at  632.8  nm  after  applying  a  closed  loop  controller  to  reduce 
the  wavefront  error.  The  SMT  residual  error  contains  high  spatial  frequency  print- 
through  error  caused  by  the  segment  ribs  and  low  spatial  frequency  error  caused  by  fine 
control  actuator  attachment  points  and  mirror  defects  [16].  A  two-dimensional  Fourier 
transform  of  the  SMT  residual  wavefront  error  identifies  the  prominent  spatial 
frequencies.  Figure  25  shows  the  Fourier  transform  log  plot  and  the  dominant  spatial 
frequencies  are  below  25  cycles/aperture. 
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Figure  25.  SMT  segment  residual  OPD  (left),  two-dimensional  FFT  log  plot  of 
SMT  residual  surface  error  (right),  from  [16] 

2.  Carbon  Fiber  Reinforced  Polymer  Mirror  Residual  Error 

The  1 -meter  CFRP  mirror  developed  by  CMA  provides  a  reference  of  expected 
surface  error  for  a  large  CFRP  mirror.  The  mirror  has  troubled  areas  on  the  surface  due 
to  manufacturing  defects.  For  example,  the  center  core  of  the  substrate  failed  when  the 
mirror  was  removed  from  the  mandrel.  The  mirror  is  characterized  using  a  Zygo 
interferometer.  Figure  26  shows  the  measured  OPD  with  piston,  tip,  and  tilt  removed. 
The  interferometer  measurements  have  dropouts  in  the  center  and  in  the  top  right 
quadrant  due  to  manufacturing  defects.  The  hexagonal  back  structure  is  also  apparent  as 
print  through. 


CFRP  Mirror 
RMS  WFE:  0.69 


CFRP  Mirror  CFRP  Mirror 

105  Zcrnikc  Terms  Zcrnike  +  Residual 

RMS  WFE:  0.64  RMS  WFE:  0.68 


Figure  26.  CFRP  measured  (left),  CFRP  Zemike  fitted  (center),  CFRP  Zemike 
fitted  plus  residual  wavefront  error  from  measurement  (right). 
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The  CFRP  mirror  is  modeled  by  fitting  the  measured  data  to  Zemike  polynomials. 
This  gives  a  good  approximation  of  the  mirror  surface  and  fills  in  the  discontinuities  due 
to  the  interferometer  dropouts.  The  difference  between  the  fitted  Zemike  surface  and  the 
measured  surface  is  found  and  the  difference  is  added  to  the  fitted  surface.  This  provides 
a  continuous  model  of  the  CFRP  surface  error  to  develop  deformable  mirror  models  to  fix 
the  residual  surface  error. 

3.  Generalized  Residual  Error  Model 

This  section  presents  a  model  developed  to  simulate  the  residual  error  of  a 
hypothetical  large  aperture  lightweight  mirror.  The  model  outputs  a  2D  OPD  that  has 
characteristics  similar  to  current  replicated  primary  mirror  aberrations  including  low 
order  aberrations,  print  through,  and  fixture  error.  The  print  through  geometry  is 
determined  by  the  rib  structure  of  the  mirror  substrate.  The  model  can  include  hexagonal 
ribs  or  rectangular  ribs,  and  can  vary  the  number  of  ribs  in  the  support  structure  to 
generate  print-through  in  the  OPD.  Figure  27  shows  an  example  of  a  simulated  OPD  for  a 
hexagonal  segment  with  hexagonal  ribs,  three  fixture  points,  and  astigmatism.  Figure  28 
shows  a  simulated  OPD  for  a  circular  aperture  with  rectangular  ribs.  This  model  can  be 
used  to  predict  the  number  of  tweeter  actuators  and  actuator  spacing  to  remove  the 
surface  errors. 
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Figure  27.  Simulated  OPD  of  hexagonal  mirror  with  hexagonal  print  through 
error,  fixture  error,  and  astigmatism  (left),  two-dimensional  FFT  log  plot  of 
simulated  residual  surface  error  (right) 


Figure  28.  Simulated  OPD  of  circular  mirror  with  rectangular  print  through 
error,  fixture  error,  and  astigmatism  (left),  two-dimensional  FFT  log  plot  of 
simulated  residual  surface  error  (right) 


E.  SMT  MODEL  DEVELOPMENT 

This  section  describes  the  SMT  model  development  and  includes  background  on 
the  original  SMT  state  space  model.  The  SMT  model  is  based  on  the  2nd  order  MIMO 
state  space  model  developed  by  the  Lockheed  Martin  and  updated  by  Yingling  [54],  This 
section  describes  the  addition  of  actuator  and  sensor  modeling  as  well  as  the  development 
of  a  state  feedback  observer  used  to  estimate  the  states  from  the  sensor  outputs. 
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1.  State  Space  Model 

The  SMT  optical  telescope  assembly  includes  the  control  system  as  part  of  the 
optical  prescription.  The  state  space  model  developed  by  Horn,  Bampton  [57],  and  later 
updated  by  Yingling  [54]  is  a  starting  point  for  the  research  in  this  dissertation.  The  state- 
space  model  represents  the  structural  and  optical  response  of  the  telescope  and  includes 
the  flexible  dynamics  of  the  six  segments  and  the  payload  support  structure.  The  models 
primary  outputs  are  a  2-dimensional  282  x  282  OPD  at  the  exit  pupil  plane  and  the  edge 
sensor  displacement  of  the  segments  relative  to  the  SMT  structure.  The  model  secondary 
outputs  are  the  wavefront  sensor  outputs.  Sensor  outputs  matrices  convert  the  OPD  to 
Shack-Hartmann  slope  measurements,  and  phase  diversity  piston  and  slopes 
measurement.  Table  1  presented  earlier  in  this  chapter  summarizes  the  model  inputs  and 
outputs. 

The  state  space  model  inputs  and  outputs  reference  either  a  local  or  global 
coordinate  system.  The  global  coordinate  system  is  with  respect  to  the  center  of 
telescope  aperture.  The  local  coordinate  system  is  centered  at  each  segment  with  the  x- 
axis  pointing  radially  away  from  the  center  of  the  telescope.  Figure  29  and  Figure  30 
show  the  SMT  global  and  local  coordinate  systems. 
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Figure  29.  SMT  global  (red)  and  local  coordinate  systems,  the  z-axis  points 
toward  the  observer,  from  [57] 


Figure  30.  SMT  global  (red)  and  local  (black)  coordinates  viewed  from  the 

side,  from  [57] 


Horn,  Bampton,  and  Yingling  used  finite  element  models  (FEMs)  of  the  segment 
and  optical  telescope  assembly  (OTA)  to  create  a  segment  and  an  OTA  state  space 
models.  To  reduce  model  computational  complexity  they  extracted  eigenmodes  within 
the  frequency  band  of  interest  from  the  FEMs  and  truncated  the  high  frequency  modes. 
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They  identified  the  truncated  high  frequency  modes  to  the  segment  FSAs  response.  To 
capture  the  FSA  modes  they  augmented  the  eignemodes  with  additional  residual  modes. 
They  then  built  a  second  order  state  space  model  from  the  FEMs  represented  by  A,  B,  C, 
and  D  matrices,  Eq  (3.5). 
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A  B 
C  D 


(3.5) 


Generalized  modal  coordinates  represent  the  states,  Equation  (3.6).  Equation 
(3.7)  shows  the  structure  of  the  A  matrix  where  An  is  a  zero  matrix,  A n  is  an  identity 
matrix,  A21  and  A22  are  diagonal  matrices,  equations  (3.10)  and  (3.11),  of  modal 
frequencies  where  (0;  is  the  frequency  and  g  is  the  damping  ratio  of  the  ith  mode. 
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Equation  (3.12)  shows  the  structure  of  the  B  matrix.  Bi  is  a  zero  matrix  and  B2 
represents  the  mode  shapes  at  the  actuator  locations.  Equation  (3.13)  shows  the  C  matrix 
structure  where  Ci  represents  the  modes  shapes  at  the  sensor  locations  and  C2  is  a  zero 
matrix.  The  D  matrix  is  nominally  a  zero  matrix  because  the  states  only  change  due  to 
structural  interactions  [54], 


[5] 


B, 

B2 


[C]  =  [  C,  C2  ] 


(3.12) 

(3.13) 
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Representing  the  high  frequency  FSA  modes  statistically  further  reduced 
complexity  of  the  state  space  model.  This  was  accomplished  by  representing  the  residual 
modes  associated  with  the  segment  FSAs  as  static  influence  functions.  The  influence 
functions  were  removed  from  the  dynamic  model  and  represented  using  the  D  matrix. 
The  B  matrix  then  represents  only  normal  mode  sensitivity  to  specific  actuator  input  and 
the  C  matrix  represents  the  normal  modes  shapes.  The  dynamic  state  space  model  from 
equation  (3.5)  is  rewritten  into  two  state  space  equations  below. 

{x(,)}=WW0}+[B]{»(>)}  (3.14) 

M>)}  =  [C]M<)}+[o]{«(<)}  (3.15) 

A  static  influence  matrix  model  from  the  dynamic  model  is  found  by  assuming  a 
steady  state  response,  {i(f)}  =  0  .  Substituting  equation  (3.16)  into  equation  (3.17)  directly 

relates  the  inputs  to  the  measured  output.  The  matrix  T  is  the  analytical  influence  matrix 
from  the  static  model.  The  static  model  is  used  in  cases  where  the  telescope  structural 
dynamics  are  slow  compared  to  the  FSA  response.  An  example  is  a  case  of  thermal 
disturbances  where  the  dynamics  are  slow  allowing  and  structural  dynamics  are  neglected 
while  controlling  the  face  sheet  actuators.  This  quasi-static  approach  simplifies  the 
control  of  the  segment  face  sheets  as  long  as  the  bandwidth  is  low  frequency;  /  <k  i Hz 


[54]. 

0  =  [A]{x(i)}  +  [b]{m(i)}  =>  {jc(r)|  =  -[A]-1  [b]{m(i)}  (3.16) 

{>’(?)L^_^  =  [cKJC(f)}  +  [Z)]{M(?)}  =  -[c][4r1[fi]{M(0}  +  [Z)]{«(0}  (3-17) 

r^Hcpn*]  (3-18) 

{y(0}  =  r{«(0}  (3.19) 


2.  Actuator  Models 

MathWorks  Simulink  is  used  to  perform  computer  simulation  of  the  state  space 
model  and  quasi-static  state  space  model.  The  state  space  opto-mechanical  model  models 
the  structural  dynamics  of  the  SMT.  Time  delays  and  noise  are  added  to  the  CCA,  FCA, 
FSA,  and  FSM  inputs  of  the  opto-mechanical  model  to  simulate  the  actuators.  The  FSM 
torque  inputs  to  the  optical  mechanical  model  are  found  by  using  an  additional  2nd  order 
model  and  the  jitter  sensor  inputs.  The  state  space  model  is  updated  to  include  the  plant 

noise  where  {v(t)}  is  the  plant  noise  and  [G]  is  the  noise  gain  matrix. 
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{i(r)}-[A]{x(f)}  +  [5]{w(f)}  +  [G]{v(f)} 


(3.20) 


3.  Sensor  Models 

The  jitter  sensor  model  uses  the  OPD  output  matrix  and  a  jitter  matrix  to  compute 
the  jitter  outputs.  The  outputs  are  the  x  and  y  slopes  measured  in  radians  of  the  OPD  in 
global  coordinates.  These  outputs  combine  with  the  FSM  reaction  torques  in  the  x  and  z 
moments  and  the  FSM  angle  commands  in  the  x  and  z  rotations  both  in  global 
coordinates  to  control  the  FSM. 


The  SMT  Shack-Hartmann  wavefront  sensor  model  converts  the  OPD  output  to 
Shack-Hartmann  wavefront  sensor  outputs  by  multiplying  the  OPD  output  matrix  by  the 
wavefront  sensor  matrix,  WFS.  This  matrix  multiply  approximates  the  gradient  of  the 
wavefront  measured  by  the  Shack-Hartmann  wavefront  sensor.  Each  segment  has  61 
lenses  that  produce  61  slopes  measurements  in  the  x  direction  and  61  slope  measurements 
in  the  y  direction.  The  WFS  matrix  multiply  approximates  the  slopes  using  central 
differencing  based  on  the  location  of  the  lens  on  the  OPD  in  both  the  x  and  y  direction. 
Equations  (3.21)  and  (3.22)  find  the  slopes  from  the  OPD  given  the  displacement  values 
at  pixels  i  and  j.  The  pixels  i  and  j  correspond  to  the  Shack-Hartmann  lens  locations 
superimposed  on  the  OPD,  Figure  31. 
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(3.21) 

(3.22) 


The  Shack-Hartmann  wavefront  sensor  provides  a  large  capture  range  of  the 
segment  surface  figure,  but  cannot  determine  segment  piston  error.  The  lenslet 
developed  for  the  SMT  only  has  61  lenslets  per  a  segment.  There  are  fewer  lenslets  than 
actuators  and  high  order  surface  aberrations  cannot  be  measured.  Figure  31  shows  the 
Shack-Hartmann  lens  locations  overlaid  on  a  SMT  segment  with  the  face  sheet  actuator 
locations  represented  by  red  dots.  This  sensor  can  initially  correct  the  surface  figure  and 
the  phase  diversity  wavefront  sensor  supports  fine  surface  figure  control. 
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Figure  31.  Shack-Hartmann  wavefront  sensor  lens  locations  on  single  SMT 

segment,  from  [57] 

The  SMT  phase  diversity  sensor  model  computes  the  local  segment  tip,  tilt,  and 
piston  from  the  OPD  output  using  a  matrix  multiply.  This  approach  simplifies  the  model 
by  not  simulating  the  phase  retrieval  process.  However,  the  phase  diversity  sensor  can 
retrieve  the  entire  phase  of  the  wavefront  and  not  just  tip,  tilt,  and  piston.  The  model  is 
updated  to  include  a  phase  retrieval  algorithm.  To  generate  a  point-spread  function  at  the 
image  plane  the  OPD  output  is  used.  Defocus  is  added  to  the  image  to  use  the 
Gerchberg-Saxton  routine  to  retrieve  the  phase.  Figure  32  shows  a  phase  retrieval 
example.  The  top  row  is  generated  from  the  state  space  model  output  OPD,  the  second 
row  is  the  estimated  point  spread  function  and  pupil  estimate  using  a  phase  retrieval 
algorithm,  and  the  third  row  is  the  difference  between  the  estimate  and  the  model  output. 
The  sensor  models  output  both  the  sensor  data  and  include  wavefront  reconstruction.  The 
model  outputs  the  reconstructed  wavefront  as  an  OPD  image  and  modal  coefficients. 
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Figure  32.  Example  of  phase  retrieval  from  OPD  using 
Gerchberg-Saxton  algorithm 


The  sensors  are  further  modeled  by  adding  noise  and  a  time  delay  to  the  Simulink 
model.  Gaussian  readout  noise  is  added  to  the  OPD  output  to  simulate  sensor  noise  in  the 
Shack-Hartmann  wavefront  sensor.  Similarly  Gaussian  readout  noise  is  added  to  the  PSF 
images  in  the  phase  diversity  sensor  simulation.  To  simulate  the  sensor  bandwidth  time 
delays  are  added  to  the  sensor  output.  The  output  state  equation  is  updated  to  include 
measurement  noise  shown  below  where  {w(t)}  is  the  measurement  noise  and  [H]  is  the 
measurement  noise  gain. 

{>■(')}  =  [C  ]{*(>)} +[D]{«(<)} +[«]{»'(<)}  (3-23) 

4.  Observer 

A  full  state  observer  estimates  the  state  variables  based  on  the  measurements  of 
the  sensor  model  outputs  and  the  control  variables.  The  observer  in  this  model  uses 
modal  coefficients  sensor  outputs  to  represent  the  reconstructed  OPD.  The  state  variables 
of  the  state  space  model  are  generalized  modal  coordinates  and  the  output  is  the  OPD. 
The  observation  error  is  defined  as  the  difference  between  the  measured  output  and  the 

estimated  output.  The  observer  model  is  defined  in  equation  (3.24)  where  x  is  the 
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estimated  state  and  cx+du  is  the  estimated  output.  Collecting  the  states,  controls,  and 
outputs  gives  equation  (3.25). 

{i}  =  [A]{j}  +  [B]{«}  +  [s:,]({,}-[C]{i}-[Z)]W)  (3.24) 

{«H|4H^][C]){.v}+([b]-[k.][o]){«}+[s.]W  (3.25) 

The  matrix  Ke  is  the  observer  gain  matrix.  This  gain  matrix  weights  the 
correction  term  involving  the  difference  between  the  measured  output  and  the  estimated 
output.  The  error  vector  is  defined  by  (3.26)  and  the  observer  error  equation  is  given  by 
(3.27).  The  matrix  Ke  is  selected  such  that  [A-KcC]  is  a  stable  matrix  and  the  error  vector 
converges  to  zero.  This  is  accomplished  by  finding  the  steady  state  gain,  Ke,  using  the 
MATLAB  kalman.m  function.  The  gain  matrix  Ke  provides  the  optimal  solution  to  the 
steady  state  observer  estimation  problem  given  the  state  equations  (3.20)  and  (3.23)  with 
known  inputs  and  Gaussian  plant  noise  and  measurement  noise.  The  gain  matrix  Ke  is 
used  to  find  the  state  estimate,  {*},  that  minimizes  the  steady  state  error  covariance 
shown  in  equation  (3.28). 


(3.26) 

M=(M-[*r.][c])M 

(3.27) 

*  =  UmE^[x- x}{x- x}T  j 

(3.28) 

F.  DEFORMABLE  MIRROR  SIZING  AND  MODELING 

This  section  determines  the  tweeter  deformable  mirror  actuator  count  and  spacing 
required  to  compensate  for  the  primary  mirror  residual  error.  The  actuator  count  and 
spacing  is  used  to  model  tweeter  deformable  mirror  using  an  influence  matrix.  First  the 
methodology  used  for  sizing  the  deformable  mirror  is  presented.  Next  deformable  mirror 
models  are  presented  to  correct  the  residual  error  of  the  SMT  AHM  segment  and  the  1- 
meter  CFRP  mirror.  Finally  a  model  of  the  Boston  Micromachine  Corporation  (BMC) 
Multi-DM  used  in  the  experimental  setup  presented  in  Chapter  VI  is  described. 

1.  Sizing  the  Deformable  Mirror 

The  deformable  mirror  actuator  count  and  spacing  depends  on  the  residual  error 
being  corrected.  The  desire  is  to  reduce  the  residual  RMS  wavefront  error  to  A./15  waves 
RMS.  Scholl  and  Lawrence  [73]  investigated  the  feasibility  of  using  a  deformable  mirror 
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to  correct  quasi-static  changes  of  on-orbit  primary  mirrors  and  showed  that  to  achieve  a 
Strehl  ratio  of  0.95  the  actuator  spacing  corresponded  to  approximately  half  the  Nyquist 
frequency.  The  first  step  to  sizing  the  deformable  mirror  is  to  determine  the  dominant 
spatial  frequencies  using  the  2D  FFT  residual  error  plots.  This  spatial  frequency 
information  determines  the  number  of  actuators  and  actuator  spacing  of  the  tweeter  in 
order  to  spatially  filter  the  dominant  frequencies. 

Next  the  deformable  mirror  influence  function  is  modeled.  Two  common 
approaches  to  modeling  a  continuous  surface  normal  deformable  mirror  actuator 
influence  functions  is  to  use  a  cubic  function  or  Gaussian  function  [48].  This  research 
uses  a  Gaussian  function  to  model  the  actuator  influence  to  present  a  generalized 
deformable  mirror  instead  of  a  specific  mirror  technology  like  bimorph  piezo-electric  or 
MEMs.  The  Gaussian  function  was  chosen  for  simplicity  but  does  not  model  deformable 
mirror  edge  constraints  as  noted  by  [74],  Equation  (3.29)  represents  a  two-dimensional 
Gaussian  influence  function  where  rs  is  the  actuator  spacing,  a  is  the  actuator  coupling 
coefficient,  b  is  the  Gaussian  index,  and  S  is  the  stroke  of  the  mirror  deflection.  This 
function  was  modified  to  generate  a  3D  response  that  represents  the  2D  OPD  output  of  a 
single  actuator. 

'  /  /|  |  \ y N 

IFGaussian{r)  =  S  eXP  ln(«)  ~ “  (3-29) 

X  X  ^  ''  J 

The  actuators  were  modeled  over  a  square  grid  with  the  actuator  spacing  matching 
the  dominant  spatial  frequency  found  using  the  2D  FFT.  The  actuator  spacing  was 
assumed  to  be  constant  over  the  mirror.  To  represent  the  hexagonal  SMT  aperture  shape 
or  circular  CFRP  mirror  shape  an  aperture  mask  was  applied  to  the  square  actuator  grid 
and  the  unused  actuators  were  removed.  The  influence  function  stroke  was  set  to  2.5 
times  the  RMS  wavefront  error.  Figure  33  shows  an  example  of  a  rectangular  actuator 
grid  with  a  SMT  segment  hexagonal  mask,  and  Figure  34  shows  an  example  of  a 
rectangular  actuator  grid  with  the  SMT  pupil  plane  mask.  The  deformable  mirror 
influence  matrix  model  is  built  from  the  individual  influence  function  responses. 
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Figure  33.  Deformable  mirror  model  with  333  actuators  in  rectangular  grid 
(left),  model  influence  function  for  actuator  100  (right) 

Actuator  100  Influence  Funetion 


Figure  34.  Deformable  mirror  model  with  1778  actuators  in  rectangular  grid 
and  SMT  pupil  plan  mask  (left),  model  influence  function  for  actuator  100  (right) 

2.  Deformable  Mirror  Model  for  SMT  AHM 

A  2D  Fourier  transform  of  the  SMT  residual  wavefront  error  identified  the 

prominent  spatial  frequencies.  Figure  25  shows  the  Fourier  transform  log  plot  and  the 

dominant  spatial  frequencies  are  below  25  cycles/aperture.  The  corresponding  actuator 

spacing  results  in  a  25  x  25  actuator  grid  for  the  tweeter.  The  actuator  spacing,  response 

function  from  equation  (3.29),  and  the  SMT  hexagonal  mask  results  in  an  influence 
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0.2, 


matrix  model  with  a  total  of  543  active  actuators.  The  actuator  stroke  is  set  to  2.5  times 
the  RMS  wavefront  error. 

The  tweeter  deformable  mirror  model  is  included  in  a  simulated  control  loop  to 
correct  the  residual  SMT  AHM  residual  error  model.  The  control  loop  is  a  global 
iterative  constrained  minimization  controller  and  presented  in  Chapter  IV.  The  closed 
loop  controller  can  use  either  the  measured  OPD  or  SVD  modal  coefficients  as  feedback. 
The  SVD  modes  represent  the  spatial  modes  of  the  woofer- tweeter  mirror  system.  Figure 
35  shows  the  error  history  using  OPD  feedback  and  SVD  modal  feedback  including  the 
first  25  SVD  modes,  50  SVD  modes,  75  SVD  modes,  and  100  SVD  modes.  The  543 
channel  tweeter  achieves  a  RMS  wavefront  error  below  A/ 15  and  reduces  the  peak-to- 
valley  wavefront  error  below  1  wave.  The  simulation  also  shows  that  including  higher 
order  SVD  modes  in  the  feedback  improves  the  wavefront.  The  two-dimensional  FFT 
log  plot,  Figure  36,  shows  that  the  deformable  mirror  model  successfully  removed  the 
dominant  spatial  frequencies. 


Simulated  Iterative  CLS  Control 
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Figure  35.  Simulated  543-channel  tweeter  with  global  constrained  least  squares 
controller,  RMS  wavefront  error  history  (left),  and  peak-to-valley  wavefront  error 

history  (right),  from  [16] 
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Figure  36.  2D  FFT  log  plot  of  simulated  dual  deformable  mirror  residual  surface 

3.  Deformable  Mirror  Model  for  CFRP  Mirror 

The  same  deformable  mirror  sizing  approach  used  for  the  SMT  AHM  mirror  is 
used  for  the  CFRP  mirror.  A  2D  Fourier  transform  of  the  SMT  residual  wavefront  error 
identified  the  prominent  spatial  frequencies  contributing  to  the  residual  error.  Figure  37 
shows  the  Fourier  transform  log  plot  and  the  dominant  spatial  frequencies  are  below  75 
cycles/aperture.  The  corresponding  actuator  spacing  results  in  a  75  x  75  actuator  grid  for 
the  tweeter  [16].  This  actuator  count  is  very  large  and  greater  than  the  actuator  count  of 
commercially  available  deformable  mirrors. 

Instead  of  correcting  with  an  unrealistic  deformable  mirror  model  the  deformable 
mirror  sizing  process  was  incorporated  into  a  simulation  where  the  number  of  actuators 
was  considered  a  variable.  This  allowed  multiple  deformable  mirror  models  to  be  tested 
against  the  CFRP  mirror  residual  wavefront.  This  simulation  used  equal  actuator  spacing 
and  the  response  function  from  equation  (3.29).  Figure  38  shows  a  plot  of  predicted 
residual  wavefront  error  versus  deformable  mirror  actuator  count  for  the  CFRP  mirror, 
and  Figure  39  shows  the  wavefront  error  after  using  a  437  actuator  deformable  mirror. 
These  results  show  significant  amount  of  wavefront  error,  but  a  significant  amount  of 
improvement  is  seen  using  as  little  as  a  32-actuator  deformable  mirror  as  shown  in  Figure 
38.  The  addition  of  more  actuators  continues  to  improve  the  wavefront,  but  the  high 
spatial  frequency  aberrations  require  an  impractical  number  of  actuators  to  further  reduce 
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the  wavefront  error.  A  better  approach  is  to  further  improve  the  CFRP  mirror  high  spatial 
frequency  surface  errors  and  allow  the  deformable  mirror  to  correct  the  lower  spatial 
frequencies. 


CFRP  Segment  Residual 
RMS  WFE  =  0.68465 


2D  FFT  of  CFRP  Segment 
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Figure  37.  CFRP  mirror  residual  error  (left)  and  two-dimensional  FFT  log  plot 
of  1-m  CFRP  mirror  residual  surface  error  (right) 
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RMS  WFE  vs.  Tweeter  Actuators 


Figure  38.  CFRP  mirror  rms  wavefront  error  versus  number  of  deformable 

mirror  actuators 


Figure  39.  CFRP  mirror  residual  error  after  correction  using  deformable  mirror 
with  437  actuators  (left),  and  2D  FFT  log  plot  of  corrected  CFRP  mirror  showing 

low  frequencies  removed  (right) 
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The  majority  of  the  CFRP  mirror  residual  error  is  attributed  to  the  discontinuities 
in  the  OPD  caused  by  the  measurement  dropouts,  and  the  substrate  hexagonal  print- 
through.  Future  versions  of  CFRP  mirrors  may  fix  these  problems  and  the  mirror  may  be 
closely  modeled  using  a  Zemike  approximation.  Figure  40  and  Figure  41  show  the 
Zemike  approximated  OPD  and  2-dimensional  FFT  before  and  after  correction  using  176 
actuators.  Figure  42  shows  that  the  simulation  achieves  A/20  RMS  wavefront  error  using 
a  deformable  mirror  with  as  few  as  150  actuators. 


Figure  40.  CFRP  mirror  Zemike  approximation  residual  error  (left)  and  2D 
FFT  log  plot  of  1-m  CFRP  mirror  Zemike  approximation  residual  surface  error 

(right) 


Figure  41.  CFRP  mirror  Zemike  approximation  residual  error  after  correction 
(left),  and  2D  FFT  log  plot  of  corrected  CFRP  mirror  Zemike  approximation 
showing  low  frequencies  removed  (right) 
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RMS  WFE  vs.  Tweeter  Actuators 


Figure  42.  CFRP  mirror  RMS  wavefront  error  approximated  with  105  Zemike 
modes  versus  number  of  deformable  mirror  actuators 

4.  Boston  Micromachine  Corporation  (BMC)  Multi-DM  Model 

In  addition  to  the  analytical  deformable  mirror  models  a  Boston  Micromachine 
Corporation  (BMC)  Multi-DM  with  140  actuators,  a  3.3mm  by  3.3mm  square  aperture, 
and  1.5pm  of  stroke  was  modeled.  This  deformable  mirror  is  used  in  the  experimental 
test  bed  described  in  Chapter  VI.  The  mirror  model  is  an  experimentally  derived 
influence  matrix.  The  BMC  deformable  mirror  actuators  were  mapped  to  the  SMT 
segment  and  the  unused  actuators  were  masked  and  removed  them  from  the  influence 
matrix.  Based  on  the  location  of  the  deformable  mirror  in  the  experimental  test  bed,  61 
of  the  140  actuators  had  influence  over  the  SMT  segment  wavefront. 
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G.  SMT  WOOFER-TWEETER  MODEL 

Two  woofer-tweeter  models  were  developed.  The  SMT  woofer-tweeter  model 
represents  the  entire  SMT  system  while  the  SMT  segment  woofer-tweeter  model 
represents  a  single  segment  woofer-tweeter  model.  Both  models  place  the  deformable 
mirror  at  the  exit  pupil  and  assume  segment  phasing  is  complete.  The  models  are 
described  in  more  detail  in  the  following  sections. 

1.  SMT  Woofer-Tweeter  Model 

The  SMT  and  deformable  mirror  system  model  combines  the  SMT  quasi-static 
state  space  model  and  the  deformable  mirror  model.  The  model  output  is  the  sum  of  the 
two  mirror  OPDs  and  the  input  is  vector  of  actuator  inputs  for  the  SMT  FSAs.  The  quasi¬ 
static  approach  takes  advantage  of  the  fact  the  structural  dynamics  are  slow  compared  to 
the  response  of  the  FSAs  and  the  deformable  mirror.  The  response  bandwidth  of  the 
SMT  segment  and  the  deformable  mirror  are  assumed  to  be  the  same.  The  assumption  is 
valid  because  the  wavefront  sensor  measures  the  combined  contribution  of  both  mirrors 
and  is  slower  than  the  mirrors  response  bandwidth. 

The  SMT  model  uses  the  sensor  outputs  to  reconstruct  the  wavefront  of  the 
individual  segments  and  recombines  the  segments  to  form  the  entire  aperture.  In  a 
similar  fashion  portions  of  the  deformable  mirror  are  mapped  to  a  particular  SMT 
segment  to  individually  correct  the  residual  error  of  the  segments.  The  model  allows  for 
the  SMT  FSAs  and  deformable  mirror  actuators  to  be  commanded  either  in  serial  or 
parallel  depending  on  the  controller. 

Zero  mean  Gaussian  plant  noise  and  sensor  noise  is  introduced  into  the  model. 
Plant  noise  is  modeled  as  a  disturbance  on  the  actuator  inputs  to  the  quasi-static  model. 
The  Shack-Hartmann  sensor  noise  is  introduced  at  the  model  OPD  output,  and  the  phase 
diversity  noise  is  introduced  after  converting  the  OPD  pupil  image  to  the  image  plane. 
The  manufacturer’s  actuator  and  sensor  variances  are  used  to  model  the  noise. 
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2.  SMT  Segment  Woofer-Tweeter  Model 

The  woofer-tweeter  deformable  mirror  model  represents  the  experimental  woofer- 
tweeter  system  presented  in  Chapter  VI.  The  model  reduces  the  complexity  of  the  larger 
system  model  by  only  modeling  one  of  the  six  SMT  AHM  segments.  This  is  suitable 
because  the  system  level  fine  phase  controller  corrects  the  surface  of  individual  segments. 
The  single  SMT  segment  and  BMC  deformable  mirror  are  modeled  using  experimental 
influence  matrices.  The  experimental  influence  matrices  built  using  equation  (2.6)  where 
the  columns  represent  the  influence  function.  The  influence  matrices  represent  the 
physical  response  of  the  system.  The  influence  functions  are  represented  as  either  the 
measured  OPD,  or  SVD  modal  coefficients.  The  influence  matrix  composed  of  OPD 
influence  functions  result  in  an  over-determined  system. 

The  experimental  influence  matrices  are  examined  to  determine  how  well  the 
influence  matrix  can  correct  the  mirrors.  The  SMT  AHM  and  BMC  deformable  mirror 
influence  matrix  normalized  singular  values  are  shown  in  Figure  43  and  can  identify 
poorly  observable  or  controllable  modes.  Small  singular  values  can  cause  large  gains 
when  the  influence  matrix  is  inverted  and  included  in  a  feedback  controller.  To  prevent 
the  model  from  causing  large  feedback  gains,  modes  represented  by  normalized  singular 
values  of  0.1  or  less  are  removed.  These  modes  represent  less  than  10%  of  the  peak 
magnitude  and  are  not  well  corrected  by  the  mirrors.  These  modes  also  represent  high 
spatial  frequencies.  For  the  SMT  AHM  this  includes  SVD  modes  greater  than  12  and  for 
the  BMC  deformable  mirror  it  includes  modes  greater  than  40.  Depending  on  the 
simulation  the  modal  feedback  is  removed  above  these  modes  or  the  influence  matrix  is 
truncated  to  prevent  the  small  singular  values  from  generating  large  gains  in  the  control 
matrix.  Zero  mean  Gaussian  plant  noise  is  used  in  the  system  level  model. 
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SMT  Influence  Matrix  BMC  Influence  Matrix 


Figure  43.  SMT  AHM  and  BMC  deformable  mirror  influence  matrix 

singular  values 

The  sensor  is  modeled  as  the  laser  interferometer  used  in  the  experimental  setup. 
The  laser  interferometer  measures  the  wavefront  error  at  the  pupil  plane.  The  laser 
interferometer  is  modeled  as  the  OPD  output  from  the  SMT  and  deformable  mirror  model 
with  zero  mean  Gaussian  sensor  noise.  The  laser  interferometer  output  standard 
deviation  is  very  small  at  0.002  waves  (632.8  nm)  RMS. 

H.  CONCLUSION 

This  chapter  presented  the  models  used  to  simulate  the  SMT,  the  deformable 
mirror,  and  woofer-tweeter  system.  The  system  level  SMT  dynamic  model  is  used  to 
simulate  coarse  alignment  and  coarse  phasing  control,  while  the  fine  phasing  control 
design  is  based  on  the  quasi-static  model.  The  residual  error  model  is  used  to  identify 
dominant  residual  surface  error  in  order  to  select  the  number  and  spacing  of  the  tweeter 
actuators.  The  analytical  models  developed  in  this  chapter  are  the  basis  for  the  control 
design  in  Chapter  IV  and  the  simulation  performed  in  Chapter  V. 
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IV.  SPACE  TELESCOPE  WAVEFRONT  CONTROL 


This  chapter  develops  the  control  techniques  for  a  space  telescope  with  a 
segmented  active  primary  mirror  and  an  additional  deformable  mirror  in  the  optical 
assembly.  First  the  chapter  presents  SMT  incremental  wavefront  control  techniques 
including  coarse  alignment,  coarse  phasing,  and  fine  phasing  techniques.  Next  the 
chapter  presents  techniques  to  improve  fine  phasing  residual  error.  This  includes  using 
an  additional  deformable  mirror  in  a  woofer-tweeter  configuration  in  order  to  improve  the 
fine  phasing  residual  wavefront  error,  and  updating  the  SMT  AHM  segment  model. 
Existing  closed  loop  woofer-tweeter  control  techniques  are  presented  and  a  woofer- 
tweeter  gradient  control  technique  is  introduced.  Open-loop  woofer-tweeter  control 
techniques  are  also  introduced.  Finally  an  adaptive  influence  matrix  technique  is 
presented  as  a  means  to  improve  fine  phasing  residual  wavefront  error  by  improving  the 
SMT  AHM  segment  model. 

A.  SMT  WAVEFRONT  CONTROL 

In  preparation  for  the  woofer-tweeter  fine  phasing  control  this  section  reviews  the 
wavefront  sensing  and  control  architecture  of  the  SMT.  State  feedback  control 
techniques  are  discussed  as  a  means  to  incrementally  phase  the  telescope  segments  and 
improve  the  wavefront  error.  The  incremental  control  approach  includes  segment 
stacking,  coarse  alignment,  coarse  phasing,  and  fine  phasing. 

1.  SMT  Wavefront  Control  Approach 

The  SMT  model  developed  in  Chapter  III  presents  a  traditional  MIMO  control 
problem.  The  traditional  approach  to  a  MIMO  control  problem  uses  modem  control 
techniques  like  a  linear  quadratic  regulator  (LQR)  or  linear  quadratic  Gaussian  (LQG) 
controller.  Correia  et  al.  [75]  proposed  an  LQG  control  approach  for  future  large  aperture 
ground-based  astronomical  telescopes  where  the  dynamics  of  the  mirrors  cannot  be 
ignored.  A  similar  state-feedback  control  approach  is  applied  to  the  SMT  using  multiple 
control  loops  to  perform  and  maintain  coarse  alignment,  coarse  phasing,  and  fine 
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phasing.  This  control  approach  uses  the  dynamic  and  quasi-static  models  developed  in 
the  previous  chapter. 

The  SMT  control  approach  is  an  incremental  approach  that  first  aligns  the 
segments,  then  phases  the  segments,  and  then  improves  the  wavefront  error  caused  by 
surface  imperfections  on  the  primary  mirror  segments.  Figure  44  shows  the  incremental 
alignment  and  wavefront  control  approach  for  the  SMT.  Segment  alignment  begins  after 
telescope  deployment. 

After  telescope  deployment  the  segments  are  set  in  place  using  actuator  settings 
determined  during  ground  testing.  The  initial  telescope  wavefront  errors  are  expected  to 
be  on  the  order  of  millimeters.  The  telescope  begins  the  incremental  segment  phasing 
routines  by  pointing  the  telescope  toward  a  star  to  serve  as  a  coherent  light  source.  The 
wavefront  error  of  the  SMT  incrementally  improves  until  the  total  wavefront  error  of  the 
telescope  is  less  than  A/15  rms  for  0.6  |am  light. 

State  feedback  control  is  used  for  the  coarse  alignment,  coarse  phasing,  and  the 
FCA  fine  phasing  control.  These  controllers  are  simulated  with  the  SMT  dynamic  model. 
The  FSA  control  is  simulated  using  the  SMT  quasi-static  model.  After  initial  alignment 
and  phasing  the  coarse  phasing  and  fine  phasing  controllers  provide  maintenance  control 
to  maintain  the  optical  quality  of  the  telescope  and  prevent  error  growth. 


70 


Figure  44.  SMT  incremental  alignment  and  wavefront  control  approach 


2.  Segment  Stacking 

The  first  alignment  step  identifies  the  SMT  segments  on  the  camera  focal  plane. 
If  a  segment  sub-image  is  not  on  the  focal  plane  the  segment  is  steered  using  the  CCAs 
until  the  sub-image  appears  on  the  focal  plane.  Identifying  the  sub-images  may  require 
deliberate  misalignment  in  order  to  identify  them  individually  at  the  focal  plane.  To 
associate  the  sub-images  to  a  segment  they  are  individually  moved  and  new  images  are 
taken  and  compared  to  the  previous  images.  Figure  45  shows  the  initial  unaligned  sub¬ 
images  and  Figure  46  shows  a  composite  history  of  segment  identification.  The  centroid 
location  of  the  segment  sub-image  on  the  focal  plane  provides  feedback  to  the  stacking 
control  law.  The  stacking  control  law  actuates  the  CCAs  to  drive  the  six  sub-images  to 
the  center  and  reduces  the  alignment  error  from  millimeters  to  10s  of  micrometers. 
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Figure  45.  Unaligned  SMT  segments  OPD  at  exit  pupil  (left),  and  unstacked 
segment  images  at  image  plane  (right) 


Stacked  Segments 


Figure  46.  SMT  sub-image  segment  identification  (left),  and  stacked  segment 

sub-images  (right) 

3.  Coarse  Alignment 

Coarse  alignment  coarsely  aligns  the  segments  using  the  edge  sensors  and  the 
CCAs.  The  SMT  development  program  showed  that  the  edge  sensors  originally  used  on 
the  program  were  not  accurate  enough  to  phase  the  segment  edges.  NASA  JPL  developed 
a  laser  truss  metrology  system  to  aid  in  the  phasing  of  segmented  telescopes.  The  laser 
truss  metrology  system  could  provide  continuous  feedback  to  the  segment  hexapod 
actuators  to  capture  and  maintain  phasing. 


Segment  Identification 
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The  SMT  model  uses  the  edge  sensors  to  build  a  LQG  state  feedback  controller  to 
drive  the  edge  sensor  error  to  zero.  The  edge  sensors  output  is  a  differential  displacement 
in  the  global  z  direction  between  the  segment  and  the  center  core  structure.  The  model 
introduces  plant  noise  at  the  actuator  and  edge  sensor  measurement.  The  state  space 
model  from  Chapter  III  is  modified  to  only  include  edge  sensor  measurements  and  CCA 
commands.  The  state  equations  are  updated  to  include  zero  mean  Gaussian  plant  noise, 
w,  and  sensor  noise,  v.  The  matrices  G  and  H  are  the  respective  covariance  matrices.  We 
use  a  Kalman  estimator  to  predict  the  state  where  L  is  the  Kalman  gain,  equation  (4.3). 


{i}  =  [A]{x}  +  [fl]{«}  +  [G]{W} 

(4.1) 

{y}  =  [C]{x}  +  [D]{u}  +  [H]{W}  +  M 

(4.2) 

{i}  =  [A]{x}  +  [5]{M}  +  [L](W-[C]{x}-[D]{M}) 

(4.3) 

The  state  feedback  control  law  (4.4)  minimizes  the  quadratic  cost  function  (4.5) 
subject  to  the  system  dynamics,  (3.14).  The  Q  matrix  weights  the  error  while  the  R  matrix 
weights  the  amount  of  control  expended.  The  Q  and  the  R  matrices  are  both  positive 
definite.  To  find  the  optimal  gain  K  the  linear  control  law  is  substituted  into  the  dynamics 
and  cost  function  and  then  the  optimization  problem  is  solved  [76].  The  matrix  K  must 
ensure  that  the  eigenvalues  of  [A-BK]  have  real  and  negative  parts  to  ensure  the  system  is 
stable. 


{u}  =  -[K]{x} 

(4.4) 

o 

(4.5) 

{. x}  =  [A-BK]{x } 

(4.6) 

After  coarse  alignment  a  Shack-Hartmann  wavefront  sensor  measurement  is  taken 
and  an  initial  surface  correction  is  applied  to  the  segment  surface  using  the  adaptive 
optics  control  approach  like  (2.10).  The  Shack-Hartmann  sensor  has  a  large  capture 
range  and  does  not  measure  segment  piston.  Tip  and  tilt  are  removed  from  the 
measurement  to  prevent  the  saturation  of  the  FSA  actuators.  The  initial  segment 
correction  improves  the  focal  plane  measurements  and  further  reduces  the  wavefront 
error. 
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4.  Coarse  Phasing 

The  coarse  phasing  control  loop  further  improves  phasing  using  the  FCAs.  Both 
the  edge  sensors  and  phase  diversity  wavefront  sensor  provide  feedback.  The  phase 
diversity  sensor  removes  bias  and  drift  from  the  edge  sensor.  However,  the  segments 
must  be  phased  within  the  capture  range  of  the  phase  diversity  sensor.  A  state  feedback 
controller  and  Kalman  estimator  are  used  in  a  similar  manner  as  the  coarse  alignment 
controller  where  the  state  space  model  is  modified  to  only  input  FCA  commands  and 
output  edge  sensor  measurements. 

The  SMT  system  depends  on  the  edge  sensors  and  coarse  controllers  to  align  the 
segments  within  the  capture  range  of  the  phase  diversity  sensor  which  is  less  than  1  wave 
peak- to- valley  [11].  In  this  simulation  the  edge  sensor  models  are  used.  However  on 
more  recent  segmented  optical  telescope  programs  like  ATLAST  and  OPTIIX  the 
preferred  sensor  technique  is  DFS.  Once  the  segments  are  aligned  within  the  capture 
range  of  the  phase  diversity  sensor,  the  phase  diversity  measurements  provide  more 
accurate  tip,  tilt,  and  piston  measurements  (<  6  nm). 

5.  Fine  Phasing 

The  fine  phasing  control  loops  include  the  FSAs  and  the  FCAs.  The  FSAs 
primarily  control  the  surface  of  the  segment.  The  FCA  actuators  are  updated  to  offload 
piston,  tip,  and  tilt  from  the  FSAs.  The  fine  phasing  control  loops  use  the  Shack- 
Hartmann  wavefront  sensor  and  phase  diversity  wavefront  sensor.  The  Shack-Hartmann 
wavefront  sensor  measures  the  low  spatial  frequency  surface  errors  at  high  rates,  while 
the  phase  diversity  wavefront  sensor  also  measures  high  spatial  frequency  surface  error 
and  piston  error  at  low  rates. 

The  wavefront  is  incrementally  improved  by  first  reducing  the  piston,  tip,  and  tilt 
errors  using  the  FCAs  and  phase  diversity  sensor.  Next  the  lower  spatial  frequency 
surface  errors  are  improved  using  the  FSAs  and  Shack-Hartmann  wavefront  sensors.  The 
segment  wavefront  errors  are  further  improved  using  the  FSAs  and  the  phase  diversity 
sensor  with  the  Gerchberg-Saxton  phase  retrieval  algorithm.  This  phase  retrieval 
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algorithm  captures  the  segment  surface  errors  that  are  unobserved  by  the  Shack- 
Hartmann  wavefront  sensor. 

The  FCA  controller  uses  the  dynamic  SMT  model  and  the  FSA  controller  uses  the 
quasi-static  model.  To  capture  realistic  residual  surface  error  of  the  SMT  segments  the 
model  uses  measured  SMT  segment  data  for  the  segment  face  sheet  initial  conditions. 
This  data  is  taken  at  the  SMT  biased  position.  The  FCA  fine  phase  controller  uses  the 
segment  piston,  tip,  and  tilt  measurements  from  the  phase  retrieval  process  in  a  state 
feedback  controller.  The  FSA  controllers  use  integral  control  using  an  influence  matrix 
calibrated  to  the  specific  segment. 

Another  FSA  fine  phasing  approach  is  to  apply  open  loop  control  voltages 
determined  from  a  constrained  minimization.  This  approach  assumes  that  the  influence 
matrix  model  accurately  predicts  the  surface  figure  of  the  SMT  segment  and  that  the 
SMT  segment  maintains  the  surface  figure  over  an  extended  period  of  time.  This 
approach  has  the  advantage  of  not  requiring  continuous  feedback  and  preventing  FSA 
saturation.  Solving  a  constraint  based  minimization  may  not  be  practical  onboard  the 
spacecraft.  However,  solving  the  minimization  on  the  ground  and  periodically  updating 
the  FSA  control  signals  is  feasible. 

B.  FINE  PHASING  IMPROVEMENTS 

This  section  introduces  three  approaches  to  improve  the  SMT  fine  phasing 
wavefront  error.  The  first  approach  uses  an  additional  deformable  mirror  in  a  closed  loop 
woofer-tweeter  configuration  to  correct  residual  wavefront  error.  The  second  approach 
uses  an  additional  deformable  mirror  in  an  open  loop  optimal  woofer-tweeter 
configuration  to  correct  residual  wavefront  error.  The  third  approach  addresses  variation 
between  the  SMT  AHM  segment  model  and  the  physical  system.  An  adaptive  influence 
matrix  technique  is  introduced  to  improve  the  performance  of  the  SMT  AHM  model  and 
improve  the  fine  phasing  residual  wavefront  error. 
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1.  Woofer-Tweeter  Closed-Loop  Control 

This  section  presents  the  closed  control  of  the  additional  deformable  mirror  in  a 
woofer-tweeter  configuration  with  an  SMT  segment  as  part  of  the  fine  phasing  control 
loop.  Existing  closed  loop  woofer-tweeter  control  techniques  are  described  with 
consideration  to  the  SMT  application,  and  a  new  woofer-tweeter  gradient  control 
technique  is  introduced.  The  woofer-tweeter  mirrors  are  controlled  simultaneously  using 
the  feedback  from  the  Shack-Hartmann  wavefront  sensor  and  phase  diversity  phase 
retrieval.  The  woofer-tweeter  control  techniques  are  developed  and  modeled  using  the 
quasi-static  SMT  segment  woofer-tweeter  model  developed  in  Chapter  III.  The  woofer- 
tweeter  control  techniques  presented  are  as  follows:  serial  control,  global  control,  offload 
control,  modal  control,  and  optimal  control. 

a.  Serial  Control 

The  serial  woofer-tweeter  controller  uses  two  sequential  iterative  constrained 
minimization  controllers  to  correct  the  woofer  and  then  correct  the  tweeter.  This 
technique  does  not  simultaneously  control  the  woofer  and  tweeter.  First  the  woofer 
correction  is  applied  until  a  steady  state  wavefront  error  is  reached  and  then  the  tweeter 
correction  is  applied. 

b.  Global  Control 

A  simple  technique  to  controlling  the  woofer-tweeter  system  is  to  build  a  global 
influence  matrix  and  control  the  system  as  a  single  device  as  suggested  in  [38].  The 
global  influence  matrix  concatenates  the  woofer  influence  matrix  (m  x  nw)  and  tweeter 
influence  matrices  (m  x  nt)  into  one  matrix.  This  is  possible  because  the  SMT  segment 
influence  matrix  and  deformable  mirror  influence  matrix  have  the  same  number  of  rows 
since  the  influence  matrices  are  measured  from  the  same  wavefront  sensing  system.  The 
result  is  a  single  influence  matrix  for  the  woofer-tweeter  system,  equation  (4.7). 

[rG]  =  [  rSM7.  rTweeter  ]  (4.7) 
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(1)  Integral  Control 

A  standard  adaptive  optics  integral  control  technique  can  use  the  global  influence 
matrix  to  form  a  control  matrix  by  taking  the  pseudo  inverse,  equation  (4.8).  Hampton  et 
al.  [38]  consider  this  method  impractical  because  the  ratio  of  actuator  motion  to  control 
signal  for  the  woofer  can  be  significantly  higher  than  that  of  the  tweeter  causing  woofer 
measurement  noise  to  be  much  greater  than  tweeter  noise.  However,  the  SMT  woofer- 
tweeter  system  combines  two  mirrors  with  similar  stroke  and  control  voltages,  and 
bandwidth  reducing  the  likelihood  of  this  concern.  The  active  SMT  segments  and  the 
Figure  47  shows  a  block  diagram  of  this  control  approach  where  r  represents  the  global 
influence  matrix. 


K+i}  =  ^[rc]th}+K} 


(4.8) 


error 
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Figure  47.  Global  woofer-tweeter  control  block  diagram 


(2)  Global  Iterative  Constrained  Minimization 

Integral  feedback  control  can  saturate  the  segment  actuators  and  second 
deformable  mirror  actuators.  To  prevent  actuator  saturation  an  iterative  constrained 
optimization  approach  minimizes  the  wavefront  and  ensures  the  control  signals  remains 
within  the  limits  of  the  actuators.  Equation  (4.9)  represents  the  woofer-tweeter  system  as 
a  single  step  linear  model.  The  output  is  in  terms  of  wavefront  measured  in  optical  path 
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difference,  {(/)},  and  {(/>o\  is  wavefront  that  represents  the  initial  condition.  A  constrained 
minimization  approach  uses  the  quadratic  cost  function  in  equation  (4.10)  subject  to  the 
lower  bounds  and  upper  bounds  of  available  control.  This  technique  is  a  least  squares 
problem  with  inequality  constraints.  An  open  loop  or  closed  loop  controller  can  apply 
the  resulting  control  signal.  The  simulation  and  experiment  use  iterative  feedback  control 
since  there  is  uncertainty  between  the  model  and  the  plant.  Equation  (4.11)  updates  the 
control  signal  iteratively  using  the  solution  to  the  constrained  minimization.  This 
requires  the  constrained  minimization  to  be  solved  at  each  step  with  new  lower  and  upper 
bounds  represented  by  uc,k .  The  gain  ji  is  chosen  between  0  and  1  [16]. 


M=[ro]M+{0o} 

argmin  J{u)  =  ^  ([rG ]{«}  +  {0O })r ([TG ]{w}  +  {^0 }) 
subject  to:  lb  <  [u]  <  ub 

K+i}=a'K.J+K} 


(4.9) 

(4.10) 

(4.11) 


(3)  Gradient  Control 


Another  global  woofer-tweeter  control  technique  uses  a  local  gradient  approach. 
A  similar  approach  was  shown  by  Zhu  et  al.  [77]  to  control  a  single  micro-machined 
deformable  mirror  using  Zemike  coefficients.  This  work  uses  a  woofer-tweeter  global 
influence  matrix,  rc ,  represented  by  singular  value  decomposition  (SVD)  mode  shape 

coefficients.  The  global  influence  matrix  is  transformed  to  a  modal  influence  matrix,  f  G , 
using  the  singular  value  decomposition. 

[rG]  =  [t/pM  (4.12) 

[fG]  =  ([rG][v])t[rG]  (4.13) 

SVD  modal  coefficients  represent  the  wavefront  using  equation  (4.14).  The 
variable  ac  is  the  wavefront  correction  and  aaber  is  the  wavefront  aberration  caused  by  the 
imperfect  mirror  surface. 

{a}  =  {ac}  +  {aaber} 

k}=[fc]{«} 

M=[fG]M+k^} 
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(4.14) 

(4.15) 

(4.16) 


Equation  (4.17)  defines  a  cost  function  as  the  square  of  the  wavefront  error  and 
equation  (4.18)  defines  the  wavefront  error.  The  variable  y  is  the  measured  wavefront  and 
the  variable  yo  is  the  desired  flat  wavefront. 


(4.17) 

{y-v}=[rG][v]M  =  [fG]M 

(4.18) 

•/=([fo]M)r([fG]W) 

(4.19) 

([k]([fo]{«}+Kk)f([k]([fG]{«}+{«^})) 

(4.20) 

Equation  (4.21)  takes  the  partial  derivative  of  the  cost  function  with  respect  to  the 
control  variable  u,  and  equation  (4.22)  simplifies  the  previous  equation.  Next  the  partial 
derivative  of  the  cost  function  is  taken  with  respect  to  the  control  variable  u,  equation 
(4.21),  and  reduces  to  equation  (4.22).  The  gradient  of  the  cost  is  incorporated  into  a 
deepest  decent  algorithm  to  iteratively  reduce  the  variance  of  the  wavefront.  The 
gradient  control  update  equation  is  shown  in  equation  (4.23).  The  steepest  decent 
approach  includes  the  positive  small  scalar  gain  a  that  is  greater  than  zero.  There  are 
many  techniques  to  choose  the  value  of  a  to  produce  a  minimum  for  the  function  J. 
However,  a  fixed  step  is  the  simplest  and  produces  reliable  results  shown  later  in  Chapter 
VI.  For  this  controller  the  gain  a  is  set  between  zero  and  the  reciprocal  of  the  largest 
eigenvalue  of  [fG]r[fG]  which  is  also  the  largest  squared  singular  value. 


^  =  (2[fG]T[fc]{M}  +  2[fcJ{aai)er})[fG]r[fc]  (4.21) 

|^  =  2[fJ([fJ[fG]M)  (4.22) 

K+1 }  =  k  } -  2«  [f  c  J  ([f  c  J  [f  G  ]{at })  (4.23) 

0  <a<f  (4.24) 

A 


A  benefit  of  this  control  technique  is  that  the  square  term,  [fG]r[fG] ,  in  equation 

(4.23)  is  an  SVD  modal  weighting  matrix.  The  feedback  vector  of  SVD  modal 
coefficients  {a}  is  weighted  by  the  square  of  the  singular  values  of  fG,  shown  in 
equation  (4.25)  and  equation  (4.26)  since  U  is  a  unitary  matrix.  This  is  beneficial 
because  the  weighting  matrix  is  specific  to  the  performance  of  the  deformable  mirrors. 
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Additionally  the  weighting  matrix  can  be  adjusted  to  remove  specific  modes  or  weight 
them  differently. 

[fcHrJv'Ht'p]  (4.25) 

rr„r  rr-.i=[-]'['-']'['-'][-]=[-]'[x]  <4-26> 

c.  Offload  Control 

Another  woofer-tweeter  control  technique  considers  the  woofer  and  tweeter 
separately  and  implements  an  offload  control  technique  [48],  [67],  [78].  This  approach 
deliberately  offloads  stroke  and  spatial  correction  based  on  the  correlation  between  the 
mirrors.  The  singular  values  of  the  two  mirrors  are  used  to  identify  the  mirror  to  offload. 
For  example,  Figure  48  shows  the  zonal  and  modal  normalized  singular  value  plot  of  the 
SMT  segment  and  BMC-140  influence  matrices.  In  this  case  the  BMC- 140  deformable 
mirror  influence  matrix  singular  values  indicate  the  mirror  can  better  recreate  more 
spatial  modes  than  the  SMT  segment.  Therefore  the  BMC-140  is  considered  the  tweeter 
and  the  low  spatial  modes  are  offloaded  from  the  tweeter  onto  the  woofer.  This 
technique  spatially  separates  the  control  between  the  woofer  and  tweeter  using  an  offload 
matrix,  O.  The  offload  matrix  is  the  projection  matrix  of  the  tweeter  onto  the  woofer. 

The  offload  matrix  is  found  by  taking  the  pseudo  inverse  of  the  woofer  influence  matrix, 
Tw,  and  multiplying  it  by  the  tweeter  influence  matrix,  Tt.  Each  column  of  the  offload 
matrix  represents  one  influence  function  of  the  tweeter  in  woofer  space.  When  used  with 
an  integral  controller  the  control  update  is  first  determined  for  the  tweeter.  The  tweeter 
control  signal  updates  the  tweeter  and  updates  the  woofer  by  projecting  the  tweeter 
control  signal  into  the  woofer  control  space  by  multiplying  the  offload  matrix,  O,  by  the 
tweeter  control  signal,  equation  (4.28)  where  {u,}  is  the  tweeter  control  vector  and {«,„}  is 

the  woofer  control  vector.  Figure  49  shows  the  offload  matrix  included  in  an  integral 
feedback  controller. 

[0]  =  [r*nrr]  (4.27) 

k}  =  [o]M  (4-28) 
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Figure  48.  Influence  matrix  singular  values,  OPD  influence  matrix  (left), 
Zemike  coefficient  influence  function  (right) 
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Figure  49.  Offload  woofer-tweeter  control  block  diagram 


d.  Modal  Control 

The  modal  woofer-tweeter  control  techniques  separate  the  wavefront  errors  based 
on  spatial  modes  like  Zemike  modes  or  SVD  modes  and  assign  the  modes  to  either  the 
woofer  or  tweeter.  This  section  presents  several  modal  control  techniques  including 
modal  coefficient  control,  confinement  correction,  and  distributed  modal  control. 
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(1)  Modal  Coefficient  Control 

A  simple  spatial  mode  control  technique  is  modal  coefficient  control  (MCC). 
This  technique  represents  the  wavefront  error  in  terms  of  modal  coefficients  and  splits  the 
error  signal  between  the  woofer  and  tweeter  by  modal  coefficients.  Zemike  polynomials 
typically  represent  the  wavefront  error  because  it  is  a  common  basis  for  both  mirrors. 
Modal  controllers  use  modal  influence  matrix  models  to  represent  the  deformable 
mirrors.  The  SVD  modal  approach  is  used  with  SMT  system  because  Zemike  modes  are 
not  orthogonal  over  the  hexagonal  shape  of  the  segments.  Figure  50  shows  the  modal 
control  block  diagram. 

The  wavefront  error  represented  by  </> ,  is  broken  into  two  separate  wavefronts,  0, 
and  </>2 .  The  wavefront  is  represented  using  Zemike  polynomials  where  an  is  the 
coefficient,  and  Zn  is  the  Zemike  term.  In  this  example  0,  represents  the  low  order  modes 
corrected  by  the  woofer  and  02  represents  the  high  order  modes. 

0  =  0i+02  (4.29) 

0,  =i>„z  (4-3°) 

n= 1 

02=i>„z„  (4-31) 

n=k+ 1 

The  controller  sends  the  low  order  modes  to  the  woofer  and  high  order  modes  to 
the  tweeter.  Two  separate  modal  woofer  and  tweeter  influence  matrices  are  used  to 
command  the  woofer  and  tweeter  mirrors  using  an  integral  modal  feedback  controller. 
This  control  technique  spatially  separates  the  control  between  the  two  mirrors.  A 
downside  to  this  technique  is  that  the  mirrors  cannot  share  stroke  to  correct  the  same 
mode.  Therefore  the  deformable  mirrors  need  to  be  independently  sized  to  correct  the 
anticipated  errors. 
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Figure  50.  Modal  control  block  diagram 


(2)  Confinement  Correction  Algorithm 

Another  modal  control  technique  is  the  confinement  correction  algorithm  from 
[79].  This  technique  skips  the  modal  estimation  of  the  wavefront  error  by  directly 
splitting  the  gradient  sensor  error  between  the  woofer  and  tweeter  using  a  direction 
gradient  method.  This  is  accomplished  by  building  a  tweeter  influence  matrix  that 
concatenates  a  gradient  response  influence  matrix  and  modal  influence  matrix.  The 
resulting  tweeter  command  voltages  are  used  to  compute  the  high  order  gradient  error 
corrected  by  the  tweeter  and  compute  the  low  order  gradient  error  to  be  corrected  by  the 
woofer.  This  technique  also  uses  integral  feedback  control. 

The  matrix  H  is  used  to  convert  the  OPD  outputs  to  slopes  using  the  61  lenslet 
Shack-Hartmann  wavefront  sensor  geometry.  This  results  in  122  outputs,  61  slopes  in 
the  x  and  y-axis.  The  influence  matrices  are  transformed  to  sensor  space  by  multiplying 
by  H.  The  gradient  and  modal  tweeter  reconstruction  matrix  used  in  the  confinement 
correction  algorithm  is  found  by  concatenating  the  sensor  space  influence  matrix,  rTSlopes , 

and  the  modal  influence  matrix,  f  T .  Figure  51  shows  the  block  diagram  of  the 
confinement  correction  algorithm.  The  confinement  correction  controller  measures  the 
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residual  error  between  the  error  signal  and  tweeter  correction  and  sends  the  residual  error 
to  the  woofer. 


[s}  =  [H][OPD] 

(4.32) 

r„es]=[tf][r„,] 

(4.33) 

rr_5toJ)es]  =  [tf][rr] 

(4.34) 

r 

-i  T  _  Slopes 

t_cca\=  * 

l  T 

(4.35) 
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Figure  51.  Confinement  correction  algorithm  controller  block  diagram 


(3)  Distributed  Modal  Control 

Another  technique  is  to  remove  modes  from  the  influence  matrices  using 
projection  matrices.  Several  techniques  to  separate  woofer-tweeter  modes  using 
distributed  modal  control  (DMC)  have  been  proposed  in  [39],  [66],  [79].  This  section 
presents  a  technique  by  presented  by  Conan  et  al.  [66]. 

Specific  modes  like  tip,  tilt,  and  piston  can  be  removed  from  the  influence 
matrices  to  reduce  the  likelihood  of  saturation  due  to  these  modes.  This  can  also 
decouple  the  woofer  and  tweeter  modes  by  removing  the  influence  functions  of  one 
mirror  from  the  other.  The  resulting  influence  matrices  are  inverted  to  form  two 
reconstruction  matrices  used  in  an  integral  feedback  controller. 
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Removing  the  woofer  modes  from  the  tweeter  prevents  the  tweeter  from 
reproducing  woofer  modes.  A  projection  matrix  of  the  tweeter  in  woofer  space,  Pt,  is 
found  using  equation  (4.36)  where  the  tweeter  is  projected  into  woofer  space.  The 
tweeter  DM  modal  basis  is  then  found  using  equation  (4.37).  The  woofer  modes  are 
removed  from  the  tweeter  using  equation  (4.38)  resulting  in  a  new  influence  matrix  rr  WR . 

In  a  similar  manner  one  can  remove  the  tweeter  modes  from  the  woofer,  however  as  seen 
in  the  example  singular  value  plot,  Figure  48,  the  woofer’s  influence  on  higher  order 
modes  is  negligible. 


[pr]  =  [rlv]t[rr] 

(4.36) 

[Br]  =  [/]  — [Prf[Pr] 

(4.37) 

]  =  [rr  ][#r  ] 

(4.38) 

Removing  the  woofer  modes  from  the  tweeter  helps  prevent  tweeter  saturation.  In 
most  AO  applications  the  tweeter  operates  at  a  higher  bandwidth  than  the  woofer  causing 
the  tweeter  to  apply  the  initial  correction  and  then  allowing  the  woofer  to  slowly  remove 
the  lower  order  modes.  This  causes  the  tweeter  to  initially  saturate.  However,  this 
method  does  not  allow  the  woofer  and  tweeter  to  share  dynamic  range  for  modes  covered 
by  both  mirrors.  This  is  undesirable  for  AO  systems  where  the  combined  stroke  of  the 
woofer  and  tweeter  is  required. 

2.  Woofer-Tweeter  Open-Loop  Control 

The  global  iterative  constrained  minimization  previously  discussed  presents  the 
idea  of  applying  an  optimal  control  approach  to  the  SMT  and  deformable  mirror  control 
system.  However,  this  technique  assumes  the  quasi-static  model  of  the  SMT  segments 
and  deformable  mirror  is  an  accurate  representation  of  the  physical  hardware.  The 
constrained  quadratic  cost  function  was  already  presented  and  used  in  the  iterative 
constrained  least  squares  controller  to  prevent  actuator  saturation  while  minimizing  the 
RMS  wavefront  error.  This  section  presents  other  cost  functions. 

RMS  wavefront  error  is  a  good  characterization  of  optical  quality.  Peak-to-valley 
wavefront  error  is  also  important  especially  when  trying  to  reduce  wavefront 
discontinuities  between  segment  edges  or  reducing  the  wavefront  error  to  within  the 
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capture  range  of  a  wavefront  sensor.  Peak-to-valley  wavefront  error  is  introduced  into 
the  cost  function  below  where  </>  is  the  wavefront  represented  as  a  vector  from  equation 

(4.9),  and  ci  and  C2  are  weighting  coefficients. 


argmin  j(u)  =  ^(c,  ({(f)}T  {(/)})  +  c2  (max({0})  -  min({^}))) 
subject  to:  lb  <{;/}<  ub 


(4.39) 


The  cost  function  can  be  reformulated  to  include  a  cost  for  peak-to-valley 
wavefront  error  as  well  as  other  control  constraints.  Previous  work  by  Redding  [80] 
investigated  adding  actuator  penalties  to  the  cost  function.  These  penalties  are  intended 
to  preserve  actuator  stroke  by  penalizing  the  deviation  of  the  new  control  signal  from  the 
previous  control  signal,  vA ,  and  the  magnitude  of  the  control  change. 


argmin  J (u)  =  ^-(c,  ({^>}r  {</>})  +  c2  (max(j>}) -  min({</>}))  +  c3 ({ vA  f  { v& })  +  c4  ({«}T  {«}))  ^  ^ 
subject  to:  lb  <  {«}  <  ub 

The  cost  function  above  is  re-written  in  terms  of  the  control  variable  u  using  the 
SMT  and  deformable  mirror  influence  matrices.  The  control  signal  is  found  at  the 
stationary  point  where  the  derivative  of  J  is  equal  to  zero  and  the  constraints  are  met. 
These  points  are  found  using  MATLAB’s  Isqlin.m  function  for  the  constrained  least 
squares  problem  described  previously,  and  fmincon.m  for  the  nonlinear  cost  functions. 


3.  Adaptive  Influence  Matrix 

The  SMT  active  mirrors  are  considered  linear  systems  and  use  traditional  adaptive 
optic  control  techniques  like  an  influence  matrix  to  relate  actuator  inputs  to  sensor 
outputs.  The  model  may  vary  from  the  physical  hardware,  or  there  may  be  regions  of 
operation  where  the  active  mirror  behavior  is  not  linear.  This  research  proposes  an 
approach  to  improve  correction  of  an  active  space  telescope  primary  mirror  by  iteratively 
updating  the  influence  matrix.  This  approach  may  be  advantageous  for  improving 
wavefront  error  on  orbit  without  having  to  recollect  influence  functions  saving  precious 
collection  time. 

The  adaptive  influence  matrix  technique  is  developed  using  a  modal  influence 
matrix.  The  influence  matrix  and  sensor  outputs  are  represented  as  SVD  modal 
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coefficients.  Equation  (4.41)  is  the  SVD  factorization  of  the  influence  matrix.  Equation 
(4.42)  makes  the  influence  matrix  orthogonal  and  equation  (4.43)  converts  the  influence 
matrix  to  SVD  modal  coefficients.  Equation  (4.44)  converts  the  measured  output  to 
modal  coefficients. 


[rHi/][z][vf 

(4.41) 

['■Hll'l 

(4.42) 

[f]=[f]’[r] 

(4.43) 

M=[r]’W 

(4.44) 

The  influence  matrix  update  begins  by  predicting  the  output  using  equation  (4.45) 
and  measuring  the  output.  The  prediction  and  measurement  error  is  computed  using 
equation  (4.46).  A  quadratic  cost  function  compares  the  measured  wavefront  output  and 
predicted  wavefront  output  error,  equation  (4.47).  The  gradient  of  the  cost  function  with 
respect  to  the  influence  matrix,  equation  (4.48),  is  used  to  update  the  influence  matrix. 
The  influence  matrix  is  using  updated  using  equation  (4.49)  where  the  gradient  is 
multiplied  by  a  small  gain  and  added  to  the  previous  influence  matrix.  The  gain  p  is 
empirically  chosen  to  be  small  in  order  to  slowly  change  the  influence  matrix.  This 
process  is  repeated  several  times  incrementally  reducing  the  cost  function  J.  The  control 
algorithm  uses  the  updated  influence  matrix. 


(4.45) 

kHvHv} 

(4.46) 

(4.47) 

(4.48) 

(4.49) 

The  influence  matrix  update  is  equivalent  to  multiplying  the  static  modal 
influence  matrix  by  a  weighting  matrix,  equation  (4.50).  The  modal  basis  of  the  feedback 
is  not  updated.  Equation  (4.51)  shows  the  updated  influence  matrix  included  in  the 
controller  update  equation  from  in  Chapter  II  where  the  error  is  represented  in  SVD 
modal  coefficients. 
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(4.50) 


K+i } = <?i  [f  j+l  ]f  (H }  -  K }) + K } 


(4.51) 


C.  CONCLUSIONS 

This  chapter  presented  the  current  segmented  mirror  telescope  wavefront  sensing 
and  control  approach  for  the  SMT  and  described  the  coarse  alignment,  coarse  phasing 
and  fine  phasing  techniques.  Existing  woofer-tweeter  control  techniques  were  discussed 
and  a  new  gradient  woofer-tweeter  control  technique  using  singular  value  weighting  was 
introduced.  A  new  influence  matrix  update  technique  was  also  presented  to  improve 
wavefront  correction  due  to  variations  between  the  model  and  physical  system.  The 
controller  developments  presented  in  this  chapter  along  with  the  models  presented  in 
Chapter  III  allow  computer  simulation  of  the  control  techniques  for  the  SMT  and  the 
SMT  segment  woofer-tweeter  system.  The  results  from  these  simulations  are  presented 
in  the  next  chapter. 
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V.  SIMULATION  RESULTS 


This  chapter  presents  the  simulation  results  of  the  models  and  controllers 
developed  in  the  earlier  chapters.  First  the  chapter  presents  a  wavefront  sensing  and 
control  simulation  results  for  the  entire  SMT  using  measured  SMT  segment  surface  error 
data  in  the  initial  conditions.  Next  segment  level  fine  phasing  is  conducted  because  it 
represents  the  experimental  setup  in  Chapter  VI.  The  segment  fine  phasing  controllers 
include  an  integral  controller,  iterative  minimization  controller,  and  gradient  controller. 
The  segment  level  control  simulation  results  show  that  the  wavefront  error  could  be 
improved  using  a  deformable  mirror.  The  chapter  then  presents  the  fine  phasing 
improvement  results  including  closed  loop  woofer-tweeter  control,  open  loop  woofer- 
tweeter  control,  and  the  adaptive  influence  matrix  model  update  technique.  These  results 
are  presented  to  assess  the  capability  of  the  controllers  prior  to  applying  them  to  the 
experimental  test  bed  in  the  next  chapter. 

A.  SMT  WAVEFRONT  CONTROL 

These  results  show  the  predicted  wavefront  error  at  each  control  step  from  coarse 
alignment  to  fine  phasing.  The  simulation  is  different  from  previous  SMT  control 
simulations  because  wavefront  sensing  and  control  had  been  simulated  from  deployment 
through  fine  phasing  using  measured  segment  surface  figure  error  data  in  the  initial 
conditions.  Including  the  actual  segment  surface  errors  demonstrates  the  wavefront 
sensing  and  control  challenges  of  the  current  SMT  system  and  reinforces  the  need  to 
improve  the  primary  mirror  surface  error.  For  example,  this  simulation  shows  that  the 
current  SMT  Shack-Hartmann  wavefront  sensor  does  not  have  enough  spatial  resolution 
to  reduce  the  segment  surface  error  within  the  capture  range  of  the  phase  diversity  sensor 
(<  1  wave  peak-to-valley).  The  simulations  also  predicted  that  even  with  a  better  sensor 
the  current  SMT  segments  could  not  reduce  the  peak-to-valley  error  to  perform  phase 
retrieval  using  the  phase  diversity  sensor.  These  results  support  the  need  to  either 
improve  the  existing  SMT  segments  or  add  an  additional  deformable  mirror  to  the  optical 
path. 
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The  SMT  wavefront  sensing  and  control  simulation  is  broken  into  the  three 
sequential  control  approaches:  coarse  alignment  control,  coarse  phasing  control,  and  fine 
phasing  control.  The  SMT  wavefront  sensing  and  control  simulation  begins  after 
segment  stacking.  The  resulting  RMS  wavefront  error  after  segment  stacking  is  on  the 
order  of  100s  of  |im.  Figure  52  shows  a  pre-stacking  OPD  with  a  RMS  wavefront  error  of 
49.8  waves  and  the  post  staking  OPD  with  a  RMS  wavefront  error  of  48.5  waves.  The 
stacking  process  drives  the  tip  and  tilt  of  the  segments  so  their  sub-images  overlap.  The 
coarse  alignment  phase  continues  to  align  the  segments  in  tip  and  tilt  as  well  as  removes 
the  piston  mismatch  between  the  segments. 


Unstacked  OPD 
RMS  WFE  =  49.8 
P-V  WFE  =  152 


Stacked  OPD 
RMS  WFE  =  48.5 
P-V  WFE  =  119 


Figure  52.  Alignment  and  segment  stacking  OPD  before  (left)  and  after  (right) 


1.  Coarse  Alignment 

The  stacked  sub-image  in  Figure  53  shows  a  large  amount  of  blurring  at  the  focal 
plane.  This  is  due  to  the  segment  piston  error  and  surface  imperfection  of  the  individual 
segments.  An  initial  surface  correction  is  applied  to  the  segments  using  the  Shack- 
Hartmann  wavefront  sensor  and  the  segment  FSAs.  The  Shack  Hartmann  wavefront 
sensor  does  not  observe  the  piston  mode  allowing  us  to  correct  for  the  SMT  segment 
surface  without  correcting  the  piston.  Figure  54  shows  the  resulting  OPD  after  surface 
correction  and  the  image  at  the  focal  plane  with  much  less  scattering  shown  in  the 
previous  figure.  Figure  55  shows  the  OPD  with  piston  removed  before  and  after  surface 
figuring  using  the  Shack-Hartmann  wavefront  sensor. 
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Segment  Identification  Stacked  Segments 


Figure  53.  Segment  identification  (left)  and  stacked  segment  sub-images  (right) 


Stacked  OPD  with  Surface  Correction 

RMS  WFE  =  48.5  Stacked  Segments  with  Surface  Correction 


Figure  54.  Initial  segment  surface  correction  using  the  Shack-Hartmann 
wavefront  sensor  and  segment  FSAs 
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SMT  OPD  without  piston  SMT  OPD  without  piston  after  surface  figuring 


RMS  WFE  =  0.423  RMS  WFE  =  036 


Figure  55.  OPD  with  piston  removed  before  surface  figuring  (left)  and  after 

surface  figuring  (right) 

The  next  step  uses  the  CCAs  with  feedback  from  the  edge  sensors.  The  dynamic 
model  and  LQG  controller  developed  in  Chapter  III  and  IV  are  used  to  drive  the  edge 
sensor  error  toward  zero.  The  Q  matrix  is  set  as  a  diagonal  matrix  of  the  dynamic 
model’s  A  matrix  singular  values.  This  places  greater  emphasis  on  the  higher  order  mode 
shapes  when  correcting  the  SMT  surface.  The  R  matrix  is  set  equal  to  the  identity 
matrix.  Figure  56  shows  the  resulting  OPD  after  coarse  alignment  along  with  a  time 
history  of  the  edge  sensor  displacements. 
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SMT  OPD  -  Coarse  Alignment  Edge  Sensor  Displacement 

RMS  WFE:  0.803  x  KT4  LQG  Control 

P-V  WFE:  5.52 


.51 - - - * - * - 

0  5  10  15  20 

seconds 


Figure  56.  OPD  after  coarse  alignment  using  CCAs  (left),  edge  sensor 
displacement  history  during  LQG  control  (right) 

2.  Coarse  Phasing 

The  coarse  phasing  step  begins  after  the  coarse  alignment  and  assumes  the  CCA 
coarse  alignment  continues  to  track  the  steady  state  edge  sensor  displacements.  The 
coarse  phasing  uses  the  FCA  with  feedback  from  the  edge  sensors  to  further  align  the 
segments.  Figure  57  shows  that  the  coarse  phasing  continues  to  improve  both  the  RMS 
wavefront  error  and  peak-to-valley  wavefront  error  due  to  the  fine  control  steps  of  the 
FCAs. 
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SMT  OPD  -  Coarse  Phasing 
RMS  WFE:  0J59 
P-V  WFE:  3.49 


Edge  Sensor  Displacement 
x  io'6  LQG  Control 


0  2  4  6  8  10 

seconds 


Figure  57.  OPD  after  coarse  phasing  use  FCAs  (left),  edge  sensor  displacement 
history  during  LQG  control  (right) 


3.  Fine  Phasing 

The  fine  phasing  routine  uses  both  the  Shack-Hartmann  wavefront  sensor  and 
phase  diversity  sensor  to  measure  the  OPD.  This  simulation  uses  the  SMT  segment 
quasi-static  model  and  the  control  inputs  are  found  using  a  constrained  least  squares 
routine.  The  peak-to-valley  wavefront  error  is  high  due  to  surface  imperfections  of  the 
AHM  segments.  This  does  not  prevent  phasing  the  segment  edges.  Figure  59  shows  the 
piston,  tip,  and  tilt  of  the  wavefront  prior  to  fine  phasing  and  the  edges  are  phased  to 
within  A/4.  However,  the  surface  peak-to-valley  wavefront  error  exceeds  1  wave  and  is 
not  within  the  capture  range  of  the  phase  retrieval  algorithm.  In  this  simulation  we  ignore 
this  problem  by  feeding  back  the  OPD  output  of  quasi-static  dynamic  model  to  predict 
the  performance  of  the  segments. 

There  are  two  problems  causing  the  high  peak-to-valley  wavefront  error.  The 
first  is  the  capability  of  the  SMT  AHM  mirror  and  the  second  is  the  Shack-Hartmann 
mask  spatial  resolution.  The  Shack-Hartmann  spatial  resolution  is  not  sufficient  to 
estimate  the  segment  high  spatial  frequency  residual  errors.  However,  when  using  the 
OPD  from  the  quasi-static  dynamic  model  as  feedback,  the  peak-to-valley  wavefront 
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error  still  exceeds  1  wave.  This  shows  that  the  current  system  segments  do  not  have  the 
actuator  density  to  correct  the  segment  residual  errors. 


Shack-Hartmann  Correction 
RMS  WFE  =  0.356 
P-V  WFE  =  3.03 


Phase  Retrieval  Correction 
RMS  WFE  =  0.0917 
P-V  WFE  =  1.5 


Figure  58.  OPD  after  surface  correction  using  Shack-Hartmann  wavefront 
sensor  (left),  OPD  after  surface  correction  using  phase  retrieval 
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0.1 

0.05 

? 
c 
00 
fN 

0  ^ 

Yl 

> 

53 

s 

-0.05 

-0.1 

Figure  59.  SMT  OPD  piston,  tip,  tilt  only  prior  to  fine  phasing 

B.  SMT  SEGMENT  CONTROL  SIMULATION 

The  SMT  test  bed  can  only  measure  a  single  AHM  segment  in  the  experimental 

setup.  Therefore  this  section  presents  control  simulations  of  a  single  segment  using 

integral  control,  iterative  constrained  minimization  control,  and  gradient  control.  These 

simulations  are  used  as  a  baseline  to  compare  the  performance  of  the  single  segment 
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woofer-tweeter  simulations.  The  results  below  show  the  constrained  minimization  and 
gradient  control  technique  outperforms  the  integral  controller.  The  controllers  are 
simulated  using  100  Monte-Carlo  runs  where  the  noise  is  varied.  The  update  gain  for  the 
integral  control  and  iterative  minimization  was  0.2. 

Figure  60  shows  the  mean  error  history  with  the  standard  deviation  represented  as 
error  bars.  The  integral  controller  and  iterative  constrained  minimization  technique  have 
similar  results  with  the  integral  controller  experiencing  integral  windup  toward  the  end  of 
the  simulation  and  the  constrained  minimization  outperforming  the  integral  controller. 
The  gradient  control  technique  shows  improved  performance  over  the  other  two 
controllers.  Despite  the  improved  performance  of  the  gradient  control  technique  the 
overall  surface  error  is  still  about  A/5  waves  RMS  requiring  more  surface  error 
improvement  supporting  the  addition  of  a  deformable  mirror  to  the  optical  telescope 
assembly. 


RMS  WFE  Historv 


Peak-to-Vallev  WFE  Historv 


Figure  60.  Simulated  SMT  segment  error  history  using  gradient  control, 
integral  control,  and  iterative  constrained  minimization 


This  was  the  first  time  the  gradient  control  technique  was  simulated  with  a  SMT 
AHM  segment.  The  baseline  simulation  results  show  that  the  gradient  control  technique 
improved  the  steady  state  SMT  AHM  mean  RMS  surface  error  by  38%  and  mean  peak- 
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to-valley  surface  error  by  32%  compared  to  the  SMT  iterative  minimization  controller. 
These  simulation  results  are  significant  because  current  AHM  control  does  not  use 
gradient  control  techniques. 


C.  FINE  PHASING  IMPROVEMENT  SIMULATION 

This  section  presents  three  simulated  approaches  to  improve  the  SMT  fine 
phasing  wavefront  error.  These  simulations  assume  that  coarse  alignment,  and  coarse 
phase  are  complete.  The  first  approach  simulates  an  additional  deformable  mirror  in  a 
closed  loop  woofer-tweeter  configuration.  The  second  approach  simulates  open  loop 
optimal  control  techniques  to  correct  residual  surface  error  with  the  woofer-tweeter 
system.  The  third  simulates  the  adaptive  influence  matrix  technique  to  improve  the 
performance  of  the  SMT  AHM  model  and  improve  the  fine  phasing  residual  wavefront 
error. 


1.  Woofer-Tweeter  Closed-Loop  Control 

This  section  presents  the  closed  loop  simulated  feedback  woofer-tweeter 
controllers  from  Chapter  IV.  The  results  show  that  the  global  woofer-tweeter  controller 
and  the  gradient  controller  outperform  the  traditional  modal  and  offload  woofer-tweeter 
control  techniques  for  the  woofer-tweeter  configuration.  The  simulation  results  also 
show  the  global  woofer-tweeter  controllers  significantly  improve  wavefront  error 
compared  to  a  serial  sequential  woofer-tweeter  constrained  minimization  controller. 

The  woofer-tweeter  system  simulated  is  a  single  SMT  segment  configured  as  the 
woofer  and  the  BMC  140  actuator  deformable  mirror  as  the  tweeter.  Only  61  of  the  140 
BMC  deformable  mirror  actuators  are  used  as  described  in  Chapter  III.  The  simulations 
used  a  static  bias  as  the  initial  condition  of  the  SMT  segment.  The  influence  matrix 
models  from  Chapter  III  are  used  to  model  the  woofer-tweeter  system.  Saturation  limits 
are  placed  on  the  SMT  FSAs  and  the  BMC  deformable  mirror  actuators  limiting  the 
available  actuator  stroke.  Plant  noise  and  sensor  noise  are  introduced  by  adding  zero 
mean  Gaussian  noise  to  the  actuator  inputs  and  OPD  output. 


97 


The  woofer-tweeter  controllers  are  compared  against  the  SMT  controller 
simulations  presented  earlier,  a  serial  woofer-tweeter  controller,  and  a  tweeter  only 
controller.  The  serial  woofer-tweeter  controller  applies  iterative  constrained 
minimization  control  to  the  woofer,  and  then  applies  iterative  constrained  minimization 
control  to  the  tweeter.  The  tweeter  only  controller  is  similar  to  the  results  presented  in 
[64], 

The  controllers  are  simulated  using  100  Monte  Carlo  runs  and  the  mean  steady 
state  RMS  and  peak-to-valley  wavefront  error  are  presented.  Appendix  A  contains  error 
history  plots  showing  the  mean  error  with  standard  deviation  error  bars  for  the  different 
controller  simulations.  Figure  61  shows  the  different  controller  steady  state  mean  RMS 
wavefront  error,  and  Figure  62  shows  the  steady  state  mean  peak-to-valley  wavefront 
error. 


Figure  61.  Simulation  of  feedback  controllers  for  woofer-tweeter  system,  RMS 
steady  state  mean  wavefront  error  the  initial  error  is  shown  in  red,  single  mirror 
control  shown  in  green,  and  woofer-tweeter  control  shown  in  blue 
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Figure  62.  Simulation  of  feedback  controllers  for  woofer- tweeter  system,  peak- 
to-valley  steady  state  mean  wavefront  error  the  initial  error  is  shown  in  red,  single 
mirror  control  shown  in  green,  and  woofer-tweeter  control  shown  in  blue 


a.  Serial  Simulated  Control 

The  serial  woofer-tweeter  controller  uses  two  sequential  iterative  constrained  least 
squares  control  loops  to  correct  the  woofer  and  then  correct  the  tweeter.  Figure  63  shows 
the  simulated  serial  woofer-tweeter  controller  error  history  where  the  SMT  segment 
applied  correction  and  then  the  deformable  mirror  applied  correction.  The  final 
wavefront  error  plot  above  shows  that  that  the  serial  woofer-tweeter  controller  improves 
the  RMS  wavefront  error  more  than  the  SMT  integral  and  constrained  minimization 
control  techniques.  However  the  serial  woofer-tweeter  controller  does  not  perform  as 
well  as  the  SMT  gradient  control  or  the  tweeter  only  control  when  comparing  RMS 
wavefront  error. 

The  under  performance  of  the  serial  controller  compared  to  tweeter  controller 
indicates  that  the  initial  condition  of  the  controller  can  impact  performance.  The  tweeter 
only  control  and  the  tweeter  portion  of  the  serial  control  are  the  same,  but  the  initial 
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conditions  are  different.  The  wavefront  error  after  applying  the  woofer  controller  has 
high  spatial  frequency  errors  that  are  not  present  at  the  SMT  biased  position.  The  BMC- 
140  tweeter  cannot  correct  these  errors.  These  results  show  that  a  better  control  strategy 
can  be  applied  to  reduce  the  overall  wavefront  error. 


Serial  CLS  Simulated  Control 
RMS  WFE  =  0.245 


Serial  CLS  Simulated  Control 
P-V  WFE  =  1.79 


Figure  63.  Simulated  woofer-tweeter  serial  control 


b.  Global  Simulated  Control 

The  global  influence  matrix  control  techniques  included  the  global  integral 
controller,  global  iterative  constrained  minimization  controller  and  the  woofer-tweeter 
gradient  controllers.  All  three  global  influence  matrix  controllers  showed  significant 
improvement  in  wavefront  error  compared  to  a  serial  woofer-tweeter  constrained 
minimization  control  approach.  The  woofer-tweeter  gradient  controller  was  the  only 
global  controller  to  outperform  the  SMT  gradient  controller  in  both  residual  RMS  and 
peak- to- valley  wavefront  error.  Figure  64  compares  the  error  history  of  the  simulated 
global  influence  matrix  control  results. 
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RMS  WFE  History 


Peak-to-Vallev  WFE  History' 


Figure  64.  Simulated  woofer-tweeter  error  history  using  gradient  control, 
integral  control,  and  iterative  constrained  minimization 


c.  Offload  Simulated  Control 

The  offload  controllers  did  improve  the  RMS  wavefront  error  compared  to  the 
serial  woofer-tweeter  controller,  but  did  not  perform  as  well  as  the  global  woofer-tweeter 
controllers.  The  offload  control  technique  projected  the  tweeter  commands  into  the 
woofer  commands.  The  result  was  a  controller  that  could  not  take  full  advantage  of  the 
woofer  capabilities  limiting  the  total  stroke  of  the  system  resulting  in  poor  performance 
in  this  simulation. 

d.  Modal  Simulated  Control 

The  modal  coefficient  control  (MCC)  technique  did  not  reduce  the  RMS 
wavefront  error  below  the  baseline  serial  woofer-tweeter  controller.  The  MCC  controller 
used  Zemike  modes  as  a  common  basis  and  spatially  separated  the  correction  based  on 
Zemike  modes.  This  approach  did  not  allow  the  two  mirrors  to  work  together  to  correct 
the  same  mode  shapes  reducing  the  overall  correction  capability  available  from  the  two 
mirrors. 

The  distributed  modal  control  (DMC)  technique  was  simulated  for  two  scenarios. 

The  first  simulation  removed  the  first  7  woofer  SVD  modes  from  the  tweeter,  and  the 

second  scenario  removed  the  first  1 1  woofer  SVD  modes  from  the  tweeter.  These  two 
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simulations  performed  better  than  the  offload  and  MCC  techniques  because  the  controller 
used  the  full  capability  of  the  woofer.  Removing  more  woofer  modes  from  the  tweeter 
improved  the  performance  by  making  more  tweeter  stroke  available  to  remove  higher 
spatial  frequency  errors. 

The  coefficient  confinement  algorithm  (CCA)  was  not  simulated  because  the 
simulation  used  the  Shack-Hartmann  wavefront  sensor  as  feedback  to  measure  the 
wavefront  slopes.  Figure  65  shows  that  the  controller  drives  the  Shack-Hartmann 
wavefront  slopes  toward  zero.  However,  the  resulting  wavefront  obtained  from  the 
quasi-static  model  shows  significant  residual  wavefront  error,  Figure  66.  This  is  because 
the  spatial  resolution  of  the  Shack-Hartmann  wavefront  sensor  was  not  high  enough  to 
accurately  estimate  the  SMT  segment  and  deformable  mirror  response  between  the 
lenslet  locations. 


102 


CCA  Final  OPD 
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P-V  WFE  =  3.821 
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Figure  66.  Final  wavefront  error  (waves)  after  apply  CCA  to  the  SMT  woofer- 

tweeter  simulation 

2.  Woofer-Tweeter  Open-Loop  Control  Simulation 

This  section  presents  the  simulation  results  from  the  open  loop  control  approaches 
applied  to  a  single  SMT  AHM  segment  and  the  woofer-tweeter  system.  Three  open  loop 
minimization  controllers  were  simulated  with  cost  functions  shown  below.  These  are  the 
cost  functions  presented  in  Chapter  IV  with  the  exception  of  (5.3)  where  the  penalty  on 
the  change  in  the  control  signal  was  removed  to  simplify  the  minimization.  The  variable 
4>  is  the  wavefront  represented  as  a  vector  from  equation  (4.9).  Then  minimization  for  the 
SMT  uses  the  SMT  segment  influence  matrix,  and  the  woofer-tweeter  minimization  uses 
the  global  influence  matrix. 


argminy(M)  =  W) 


(5.1) 


subject  to:  lb  <  {  w  }  <  ub 


(5.2) 


subject  to:  lb  <{;/}<  ub 


(5.3) 


subject  to:  lb  <{;/}<  ub 
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The  solution  to  the  first  cost  function  (5.1)  was  used  as  the  initial  guess  to  the 
minimization  routine  for  the  other  two  cost  functions.  The  simulations  applied  equal 
weights  to  the  wavefront  errors  and  a  smaller  weight  to  control  penalty  so  the  cost  is  not 
dominated  by  control.  However,  these  weights  can  change  to  achieve  different  results. 
Figure  67  and  Figure  68  show  the  simulated  open  loop  results  for  both  the  SMT 
correction  and  global  woofer-tweeter  correction.  The  results  for  cost  function  (5.1) 
closely  matched  the  SMT  and  woofer-tweeter  iterative  least  squares  results  as  expected. 
The  results  using  cost  function  (5.2)  also  matched  the  weighting  where  the  RMS 
wavefront  error  was  worse  than  using  cost  function  (5.1),  but  the  peak-to-valley 
wavefront  error  was  improved.  Cost  function  (5.3)  results  were  slightly  worse  than  the 
previous  cost  function  because  control  was  penalized. 

Including  the  peak-to-valley  penalty  in  the  cost  function  was  an  improvement 
beyond  previous  work  [80],  and  relevant  to  the  SMT  since  the  current  peak-to-valley 
wavefront  error  is  outside  the  capture  range  of  the  phase  diversity  sensor.  However,  the 
simulation  predicted  that  the  experimental  woofer-tweeter  system  does  not  have  enough 
actuators  to  reduce  the  peak-to-valley  below  1  wave. 
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Figure  67.  Open  loop  cost  function  minimization  control  RMS  wavefront  error 
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Figure  68.  Open  loop  cost  function  minimization  peak-to-valley  wavefront 

error 
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3.  Adaptive  Influence  Matrix 

The  adaptive  influence  matrix  can  improve  performance  when  the  SMT  AHM 
model  varies  from  the  physical  hardware.  To  simulate  the  adaptive  influence  matrix  two 
SMT  AHM  experimental  influence  matrices  were  used.  The  influence  matrices  were 
measured  using  the  approach  described  in  equations  (2.5)  and  (2.6).  One  influence 
matrix  was  measured  at  voltages  of  +/-25  volts,  and  another  at  +/-20  volts.  If  the  mirror 
response  was  linear  and  there  was  no  uncertainty  in  the  system  the  influence  matrices 
would  be  the  same.  The  two  influence  matrices  are  compared  by  projecting  the  influence 
functions  of  one  influence  matrix  onto  the  influence  functions  of  the  other.  If  the 
influence  matrices  were  the  same  the  projection  matrix  should  be  the  identity  matrix. 
Figure  69  shows  the  projection  matrices  r^5r±20 . 


Projection  Matrix 


Figure  69.  Projection  matrix  of  +/-  20  volt  influence  matrix  on  +/-25  volt 

influence  matrix 

The  simulation  assumed  the  influence  matrix  used  to  control  the  mirror  is  r±25 , 
and  the  plant  was  r±20 .  This  represented  uncertainty  between  our  influence  matrix  and 
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our  plant.  The  simulation  ran  using  the  SMT  modal  integral  controller  and  includes  the 
updated  influence  matrix.  Figure  70  shows  the  error  history  without  updating  the 
influence  matrix,  with  updating  the  influence  matrix,  and  if  the  influence  matrix  matches 
the  plant.  The  simulated  results  show  that  updating  the  influence  matrix  improved  the 
RMS  wavefront  error,  but  slightly  increased  the  peak-to-valley  wavefront  error.  A 
challenge  was  choosing  the  influence  matrix  update  gain  and  number  of  iterations  to  not 
change  the  influence  matrix  too  much  that  the  reconstruction  gain  matrix  (  rt  )  saturated 

the  control  voltages. 


RMS  WFE  P-V  WFE 

Influence  Matrix  =  r±25  Influence  Matrix  =  T±2s 

Plant  =  r±2o  Plant  =  T±2o 


Figure  70.  Error  history  of  SMT  AHM  modal  integral  control  simulation  with 

adaptive  influence  matrix 


D.  CONCLUSIONS 

The  SMT  wavefront  control  simulation  and  SMT  segment  wavefront  control 
simulation  show  that  the  current  SMT  system  cannot  meet  the  fine  phasing  performance 
required  of  an  optical  telescope.  This  supports  the  addition  of  a  deformable  mirror  in  a 
woofer-tweeter  configuration  to  improve  the  residual  wavefront  error.  The  following 
woofer-tweeter  control  techniques  all  improved  the  RMS  wavefront  error  more  than  the 
current  SMT  AHM  constrained  minimization  controller:  serial  control,  global  control, 
offload  control,  and  modal  control. 


107 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


108 


VI.  EXPERIMENTAL  RESULTS 


This  chapter  presents  experimental  results  from  the  SMT  woofer-tweeter  test  bed. 
The  chapter  begins  by  describing  the  experimental  test  bed  at  the  Naval  Postgraduate 
School.  The  influence  functions  are  developed  experimentally  and  the  SMT  AHM  is 
experimentally  characterized  using  a  three-point  linearity  test  and  by  comparing  influence 
functions  over  the  actuator  range.  Next  the  SMT  segment  control  results  are  presented 
including  an  integral  controller,  iterative  minimization  controller,  and  gradient  controller. 
The  SMT  segment  experimental  results  are  used  evaluate  the  woofer-tweeter  results.  The 
chapter  then  presents  the  fine  phasing  improvement  experimental  results  including  closed 
loop  woofer-tweeter  control,  open  loop  woofer-tweeter  control,  and  the  adaptive 
influence  matrix  model  update. 

A.  EXPERIMENTAL  SETUP 

The  SMT  test  bed  does  not  have  a  large  diameter  coherent  light  source  to  serve  as 
a  reference.  This  prevents  the  use  of  the  telescope  wavefront  sensors.  To  overcome  this 
the  SMT  test  bed  uses  a  center  of  curvature  test  to  measure  the  SMT  primary  mirror 
surface  error  shown  in  Figure  71.  The  center  of  curvature  test  consists  of  a  4D 
Technologies  Corp.  PhaseCam  4020  laser  interferometer  with  an  f/3  diverging  lens,  a  null 
corrector,  the  SMT  primary  mirror,  and  a  hexapod  to  align  the  interferometer  and  null 
corrector  along  the  SMT  principle  axis  [16].  This  experimental  setup  does  not  use  the 
SMT  optics  after  the  primary  mirror,  but  the  primary  mirror  is  the  main  source  of 
wavefront  error  in  the  SMT  system. 

Adding  a  deformable  mirror  to  the  center  of  curvature  test  bed  tests  the  woofer- 
tweeter  configuration.  This  is  achieved  by  adding  pupil  relay  optics  to  the  center  of 
curvature  test  in  a  double  pass  configuration.  Figure  72  shows  the  center  of  curvature  test 
bed  layout  with  the  deformable  mirror.  The  telescope  pupil  plane  is  relayed  from  the  f/3 
objective  to  the  deformable  mirror  and  then  relayed  to  the  interferometer.  The  entire 
SMT  pupil  plane  does  not  fit  on  the  deformable  mirror  surface.  Therefore,  only  one 
segment  is  viewed  and  controlled  at  a  time. 
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Figure  71.  Center  of  curvature  test  bed,  from  [64] 
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Figure  72.  SMT  center  of  curvature  experimental  layout,  from  [16] 

The  laser  interferometer  is  tolerant  of  small  structural  vibrations  and  minor  air 
turbulence  in  the  test  bed  setup.  The  interferometer  measures  the  wavefront  error  of  the 
light  reflected  by  the  primary  mirror  in  a  double  pass  configuration  and  includes  error 
contributions  of  all  the  components  in  the  optical  path.  The  interferometer  generates  a 
wavefront  measurement  averaging  10  measurements  at  approximately  1  Hz.  The  sensor 
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bandwidth  is  sufficient  for  the  closed  loop  control  to  correct  static  surface  errors  of  the 
primary  mirror  and  low  bandwidth  disturbances  [16]. 

The  second  deformable  mirror  is  a  Boston  Micromachines  Corporation  (BMC) 
continuous  surface  140-actuator  MEMs  deformable  mirror.  The  mirror  has  a  3.3mm  by 
3.3mm  square  aperture,  and  a  maximum  actuator  stroke  of  1 ,8fim.  The  woofer-tweeter 
configuration  uses  the  top  SMT  segment  in  this  experiment.  The  BMC  deformable  mirror 
position  is  adjusted  to  place  the  pupil  plane  image  of  the  top  SMT  segment  on  the 
deformable  mirror.  The  BMC  deformable  mirror  actuators  are  experimentally  mapped  to 
the  SMT  segment  with  61  actuators  having  measured  influence  over  on  the  SMT  segment 
wavefront.  Figure  73  shows  the  BMC  deformable  mirror  actuator  locations  with  respect 
to  the  SMT  segment  actuator  locations  at  the  pupil  plane  [16]. 

Actuator  Locations 


Figure  73.  Deformable  mirror  actuator  locations  projected  onto  SMT  segment 
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The  SMT  AHM  and  BMC  deformable  mirror  actuators  have  non-linear  regions. 
The  SMT  AHM  actuator  actuators  are  electrostrictive  lead  magnesium  niobate  (PMN) 
with  a  range  between  0-100V  and  have  a  quadratic  non-linear  response.  However  the 
actuator  response  is  near  linear  between  the  operating  ranges  of  20V-80V.  The  BMC 
deformable  mirror  actuator  voltage  range  is  from  0V-202V.  The  relationship  between 
voltage  and  displacement  for  the  BMC  MEMs  mirror  is  quadratic.  Watson  showed  a 
linear  relationship  between  voltage  and  displacement  can  be  approximated  between  20V 
and  180V  for  this  particular  mirror  [65].  Limiters  were  applied  to  prevent  commanding 
the  mirrors  outside  these  ranges. 

A  closed  loop  control  system  is  built  by  connecting  the  laser  interferometer,  SMT 
AHM  segment,  and  deformable  mirror  to  a  control  computer.  The  laser  interferometer 
measurements  are  used  as  feedback.  The  control  computer  uses  MATLAB  to  retrieve  the 
interferometer  measurements  and  commanded  the  SMT  segments  and  the  deformable 
mirror  actuators.  Wavefront  residual  piston,  tip,  and  tilt  were  removed  from  the 
wavefront  measurements  after  aligning  the  optical  setup.  This  allowed  testing  the  woofer- 
tweeter  fine  phasing  techniques  without  actively  compensating  for  piston,  tip,  and  tilt 
with  the  SMT  FSAs  and  deformable  mirror  actuators. 

B.  EXPERIMENTAL  SYSTEM  CHARACTERIZATION 

Deformable  mirrors  are  modeled  as  static  linear  devices  as  described  in  the 
previous  sections.  However  data  from  the  SMT  test  bed  shows  that  the  AHM  mirror 
response  does  not  perfectly  match  the  linear  model.  This  section  describes  the 
experimental  influence  matrices  and  shows  actuator  coupling.  A  three-point  linearity  test 
is  presented  to  compare  positive  and  negative  influence  functions  and  shows  that  the 
mirror  response  is  not  linear.  Finally  the  influence  functions  are  measured  and  compared 
at  different  voltages  for  a  single  actuator  showing  different  responses.  These  variations 
and  non-linear  responses  support  the  use  of  feedback  control  to  command  the  SMT  AHM 
and  the  woofer-tweeter  system. 


112 


1.  Experimental  Influence  Matrices 

The  influence  matrices  were  experimentally  derived  for  the  SMT  segment  and 
BMC- 140  mirror  using  the  PhaseCam  4020  interferometer  to  measure  the  mirror  OPD 
after  poking  individual  actuators  from  a  biased  position  in  both  a  positive  and  negative 
direction.  The  woofer  influence  function  was  measured  by  poking  the  actuators  +/-  30V 
from  the  50V  bias.  The  tweeter  influence  function  was  measured  by  poking  the  actuators 
+/-80V  from  the  101V  bias.  Figure  74  shows  the  influence  function  coupling  for  the 
woofer  influence  matrix,  Tw,  and  the  tweeter  influence  matrix,  rT.  The  SMT  actuator 
influence  is  highly  coupled  while  the  BMC  Multi-DM  actuator  influence  is  coupled  to 
local  actuators.  This  is  due  to  the  surface  parallel  actuators  of  the  SMT  and  the  surface 
perpendicular  actuators  of  the  BMC  Multi-DM. 
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Figure  74.  Woofer  actuator  space  (left)  and  tweeter  actuator  space  (right) 


2.  SMT  Three  Point  Linearity  Test 

The  active  mirror  under  study  is  a  1 -meter  hexagonal  SiC  AHM  mirror  with  156 
face  sheet  actuators.  The  actuators  are  epoxied  to  the  structural  ribs  of  the  SiC  substrate 
in  a  surface  parallel  configuration  at  the  50V  bias  position.  Activating  a  single  actuator 
produces  a  local  and  global  response  on  the  mirror  surface. 


113 


The  linear  response  of  the  AHM  actuators  is  investigated  using  a  three-point 
linearity  test.  This  test  identifies  actuators  that  have  significantly  different  negative  and 
positive  responses  by  comparing  the  negative  and  positive  influence  functions  measured 
from  the  50V  bias  with  the  actuators  at  the  extremes  of  their  linear  range.  Equation  (6.1) 
and  equation  (6.2)  shows  the  positive  and  negative  influence  function  respectively  and 
the  superscript  represents  the  control  voltage  during  the  sensor  measurement.  The 
maximum  absolute  value  of  the  influence  function  represents  the  maximum  deflection  of 
the  surface  caused  by  the  actuator.  To  compare  the  deflections  between  the  positive  and 
negative  influence  functions  the  ratio  between  the  largest  and  smallest  deflection  are 
taken.  Figure  75  shows  the  ratios  for  the  156  face  sheet  actuators  and  ratios  greater  than  1 
indicate  a  non-linear  response.  Measurement  errors,  actuator  bias  position  errors,  and 
anomalies  in  the  actuator  to  substrate  bond  can  cause  the  positive  and  negative  influence 
functions  to  be  different.  The  results  show  that  25  actuators  have  a  response  one  standard 
deviation  outside  the  sampled  156  actuators. 
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Figure  75.  Three  point  linearity  test 

The  experimentally  derived  influence  functions  average  the  positive  and  negative 
response.  The  linear  model  does  not  capture  what  actuators  generate  more  mirror 
response  in  the  positive  direction  or  the  negative  direction.  Therefore  it  is  expected  that 
the  physical  system  behave  differently  than  the  linear  model. 

3.  Influence  Function  Comparison 

Figure  76  shows  influence  functions  of  actuator  156  taken  at  varying  voltages 
over  the  linear  range  of  the  actuators.  The  influence  functions  should  be  identical  for  a 
linear  system.  However,  the  normalized  mirror  responses  are  different  when  taken  at 
different  voltages  from  the  bias.  The  influence  function  measured  at  +/-5V  and  +/-10V 
has  a  different  global  response  than  the  influence  functions  measured  at  higher  voltages. 
The  influence  functions  measured  at  +/- 25  and  +/-30  volts  are  almost  identical.  The 
mirror  center  imperfections  are  also  noticeable  in  for  the  influence  functions  taken  at 

lower  voltages.  This  could  be  due  to  the  age  of  the  SMT  segments  and  changes  in  the 
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surface  coatings,  substrate  properties,  and  the  actuator  epoxy.  The  results  show  that  the 
measured  mirror  surface  response  is  not  linear  over  the  voltage  range.  These  results 
further  support  the  use  of  feedback  control  to  compensate  for  in  accuracies  in  the  linear 
model. 


Actuator  156  +/-5V  Actuator  156+/-  10V  Actuator  156+/-  15V 


Figure  76.  Influence  function  comparison  for  actuator  156 


C.  SMT  SEGMENT  CONTROL  EXPERIMENTAL  RESULTS 

This  section  presents  experimental  control  results  of  a  single  SMT  AHM  segment 
using  integral  control,  iterative  constrained  minimization  control,  and  gradient  control. 
These  experimental  results  are  used  to  compare  the  performance  of  the  experimental 
woofer-tweeter  results.  The  results  below  show  the  constrained  minimization  and 
gradient  control  approach  outperforms  the  integral  controller.  The  SMT  segment  was 
initialized  to  the  biased  position  and  the  update  gain  for  the  integral  control  and  iterative 
minimization  was  0.2. 

The  SMT  AHM  closed  loop  experiments  produced  similar  performance  to  the 
simulation.  The  constrained  minimization  control  and  gradient  control  outperform  the 
integral  control  with  the  gradient  controller  improving  the  RMS  wavefront  error  slightly 
more  than  the  iterative  constrained  minimization,  Figure  77.  The  final  wavefront  error  of 
the  constrained  minimization  control  and  the  gradient  control  is  similar  because  the 
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constrained  minimization  and  gradient  approach  cost  functions  are  both  minimizing  the 
quadratic  wavefront  error,  but  using  different  approaches  and  constraints. 

The  improved  RMS  wavefront  error  performance  using  the  gradient  controller  is 
significant  because  the  current  AHM  control  techniques  update  the  control  signal  by 
minimizing  ||r«-y||2  with  or  without  constraints.  However,  the  gradient  approach  finds  a 
local  minimum  through  a  steepest  descent  algorithm.  Figure  78  compares  the  final 
control  voltages  of  the  three  controllers  and  the  gradient  approach  has  the  fewest 
actuators  approaching  the  saturation  limit  (+/-  30  volts).  Additionally  the  gradient 
approach  has  the  lowest  quadratic  control  cost,  Figure  79. 


Figure  77.  Experimental  SMT  segment  error  history  using  gradient  control, 
integral  control,  and  iterative  constrained  minimization 
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Figure  78.  Final  control  comparison  between  SMT  iterative  constrained  least 
squares  and  gradient  control 
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Figure  79.  SMT  AHM  quadratic  final  control  cost 
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D.  EXPERIMENTAL  FINE  PHASING  IMPROVEMENT 


This  section  presents  experimental  results  based  on  the  fine  phasing  improvement 
approaches  introduced  in  Chapter  IV  and  simulated  in  Chapter  V.  This  section  first 
presents  the  experimental  results  from  using  the  woofer-tweeter  system  with  the  closed 
loop  woofer-tweeter  control  techniques  recommended  in  Chapter  V.  Next  the  section 
presents  the  experimental  results  of  the  open  loop  woofer-tweeter  optimal  control 
techniques.  Finally  the  experimental  results  comparing  a  controller  using  an  updated 
influence  matrix  are  compared  to  a  controller  using  a  standard  influence  matrix. 

1.  Woofer-Tweeter  Closed-Loop  Control 

This  section  presents  the  experimental  closed  loop  woofer-tweeter  controller 
results.  The  results  show  that  the  global  woofer-tweeter  control  approaches  and  the 
gradient  control  approaches  outperform  the  traditional  modal  and  offload  woofer-tweeter 
control  approaches  for  the  experimental  woofer-tweeter  configuration.  The  experimental 
results  also  show  the  global  woofer-tweeter  controllers  significantly  improve  wavefront 
error  compared  to  a  serial  sequential  woofer-tweeter  constrained  minimization  controller. 
Figure  80  and  Figure  81  summarize  the  steady  state  RMS  and  peak-to-valley  wavefront 
error  for  the  feedback  control  approaches  respectively.  Appendix  B  contains 
experimental  error  history  plots  showing  the  RMS  and  peak-to-valley  wavefront  error  for 
the  different  controller  experiments. 

Setting  the  SMT  segment  3  to  the  biased  position  initialized  the  woofer-tweeter 
experiments.  The  starting  RMS  wavefront  error  was  approximately  0.57  waves. 
Similarly,  setting  the  BMC  deformable  mirror  to  the  pre-calibrated  flat  position 
initialized  the  experiment  and  resulted  in  0.02  waves  RMS  of  wavefront  error  over  the 
SMT  mask.  To  integral  controller  gain  was  set  to  0.2  and  number  of  iterations  was  set  to 
20.  The  experimental  results  are  slightly  different  than  the  predicted  results  due  to 
differences  between  the  model  and  the  experimental  hardware  performance. 
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Simulated  and  Experimental  Steady  State  RMS  Wavefront  Error 
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Figure  80.  Simulated  and  experimental  woofer-tweeter  steady  state  feedback 
control  RMS  wavefront  error  comparison 


Figure  81.  Simulated  and  experimental  woofer-tweeter  steady  state  feedback 
control  peak-to-valley  wavefront  error  comparison 
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The  steady  state  woofer-tweeter  experimental  results  are  compared  with  the 
simulated  results  from  Chapter  V.  The  experimental  results  generally  match  the 
simulated  results  validating  the  models.  This  section  presents  the  following  closed  loop 
woofer-tweeter  control  results;  serial  control,  global  control,  offload  control,  and  modal 
control. 


a.  Serial  Control 

The  woofer-tweeter  serial  controller  closely  matched  the  trend  seen  in  simulation. 
The  serial  controller  outperformed  the  SMT  integral  controller  steady  state  RMS 
wavefront  error.  The  serial  control  results  were  similar  to  the  SMT  iterative  constrained 
minimization  control  and  the  gradient  control  results  and  like  in  simulation  the  tweeter 
only  control  slight  produced  a  better  RMS  wavefront  error  result.  This  is  due  to  the  high 
spatial  frequency  residual  error  after  applying  the  SMT  AHM  correction  compared  to  the 
lower  order  errors  that  are  present  when  the  factory  bias  controls  are  applied.  In  this  case 
the  second  deformable  mirror  is  more  effective  at  removing  the  residual  error  from  the 
SMT  AHM  bias  position  than  from  the  SMT  AHM  corrected  position. 

b.  Global  Control 

The  global  woofer-tweeter  experimental  RMS  wavefront  error  results  are  similar 
to  the  simulated  global  controller  results.  The  woofer-tweeter  gradient  control  approach 
produced  the  best  RMS  wavefront  error  compared  to  the  other  woofer-tweeter 
controllers.  Figure  82  compares  the  global  control  error  histories.  The  woofer-tweeter 
gradient  control  showed  a  24%  RMS  wavefront  error  improvement  over  the  SMT 
segment  gradient  control  results.  As  with  the  SMT  AHM  correction  the  global  gradient 
approach  outperformed  the  global  iterative  constraint  based  minimization  approach.  The 
total  control  cost  was  also  slightly  less  using  the  gradient  approach  than  using  the  global 
integral  or  iterative  constraint  minimization  approach  (see  Figure  83). 

The  global  gradient  approach  offers  significantly  improved  performance  over  the 

integral  controller.  This  is  because  the  gradient  control  approach  is  more  numerically 

stable  and  does  not  require  the  multiplying  the  feedback  by  the  pseudo-inverse  of  the 

global  influence  matrix.  The  global  gradient  approach  also  weights  the  modal  feedback 
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by  the  square  of  the  influence  matrix  singular  values  placing  emphasis  on  the  modes  that 
can  be  well  corrected  by  the  mirrors. 

The  global  gradient  approach  showed  similar  performance  to  the  constrained 
minimization  control  because  the  constrained  minimization  and  gradient  approach  cost 
functions  are  both  minimizing  the  quadratic  wavefront  error,  but  using  different 
minimization  approaches  and  constraints.  The  global  gradient  approach  and  constrained 
based  minimization  approach  also  converged  in  a  similar  number  of  iterations,  however, 
the  global  gradient  control  approach  was  computationally  faster  because  it  did  not  require 
using  a  numerical  optimization  solver. 
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Figure  82.  Experimental  woofer-tweeter  error  history  using  gradient  control, 
integral  control,  and  iterative  constrained  minimization 
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v  m4  Quadratic  Control  Cost,  u1  u 
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Figure  83.  Woofer-tweeter  quadratic  final  control  cost 

c.  Offload  Control 

The  two  offload  controller  experimental  RMS  wavefront  error  results  were 
slightly  worse  than  the  global  control  results  like  with  the  simulated  results.  The 
constrained  minimization  offload  controller  produced  slightly  better  results  than  the 
integral  controller  as  expected.  The  offload  iterative  constrained  minimization  controller 
did  produce  better  RMS  wavefront  error  results  than  the  SMT  segment  control  and  the 
serial  woofer-tweeter  control.  Although  the  results  were  not  as  good  as  the  global  control 
results  in  this  experiment  this  technique  is  suitable.  This  approach  would  have  produced 
better  results  if  the  woofer  and  tweeter  shared  more  of  the  same  SVD  spatial  modes.  This 
control  approach  could  work  if  the  deformable  mirrors  used  in  the  woofer-tweeter  system 
were  selected  or  designed  to  operate  in  an  offload  woofer-tweeter  configuration. 

d.  Modal  Control 

This  section  presents  the  woofer-tweeter  distributed  modal  control  results.  Two 
experiments  were  conducted  to  match  the  simulated  distributed  modal  control.  The  first 
experiment  removed  the  first  7  woofer  SVD  modes  from  the  tweeter  and  the  second 
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experiment  removed  the  first  1 1  woofer  SVD  modes  from  the  tweeter.  The  experimental 
results  did  not  show  as  much  improvement  as  predicted  by  the  simulation.  Although  the 
full  capability  of  the  woofer  was  used  in  the  experiment  a  significant  amount  of 
capability  to  correct  low  order  modes  was  removed  from  the  tweeter.  This  is  evident 
because  the  RMS  wavefront  error  becomes  worse  as  more  woofer  modes  are  de-projected 
from  the  tweeter.  This  control  approach  could  work  if  the  woofer  and  tweeter  were 
designed  and  sized  appropriately  to  complement  their  capabilities. 

2.  Woofer-Tweeter  Open  Loop  Control 

This  section  presents  open  loop  control  experimental  results  the  woofer-tweeter 
system.  Experimental  results  for  the  SMT  AHM  segment  are  included  for  comparison. 
Three  open  loop  minimization  approaches  were  tested  using  the  cost  functions  from 
Chapter  III  repeated  here  in  equation  (6.3).  Table  3  shows  the  weighting  coefficients 
applied  for  the  three  experiments.  Figure  84  and  Figure  85  compare  the  simulated  and 
experimental  RMS  and  peak-to-valley  wavefront  error  after  applying  open  loop 
corrections. 

argminJ  (u)  =  ^-(c,  ({<p}T  {(/)})  +  c2  (max({fj-)-  min({fj-))  +  c3  ({vA}r  {vA})  +  c4  ({wf  {«})) 
u  L  '  (6.3) 

subject  to:  Vo  <  {«}  <  itb 


Table  3.  Cost  function  coefficients 


cl 

c2 

c3 

c4 

Cost  Function  1 

1 

0 

0 

0 

Cost  Function  2 

1 

2 

0 

0 

Cost  Function  3 

1 

1 

0 

5.00E-04 

The  open  loop  minimization  experimental  results  do  not  closely  match  the 
simulated  results.  This  is  due  to  the  fact  that  the  linear  influence  matrix  model  response 
varies  from  the  physical  hardware  response.  However,  the  open  loop  control  does  show 
improvement  in  RMS  and  peak-to-valley  wavefront  error  from  the  starting  wavefront 
error. 
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The  different  weighting  coefficient  results  generally  match  the  predicted  results. 
For  example,  the  experimental  results  show  that  adding  a  penalty  to  peak-to-valley 
wavefront  error  was  able  to  slightly  improve  the  peak-to-valley  wavefront  error.  The 
results  also  show  that  penalizing  the  control  increases  the  RMS  and  peak-to-valley 
wavefront  error  because  less  control  authority  is  available.  These  experimental  results 
shows  that  the  SMT  AHM  and  woofer-tweeter  system  perform  better  using  a  closed  loop 
feedback  control  system. 
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Figure  84.  Comparison  between  open  loop  simulated  and  experimental  control, 

RMS  wavefront  error 
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Figure  85.  Comparison  between  open  loop  simulated  and  experimental  control, 
peak-to-valley  wavefront  error 


3.  Adaptive  Influence  Matrix 

The  adaptive  influence  matrix  experimental  results  are  presented  in  this  section. 
The  results  show  a  small  improvement  in  RMS  and  peak-to-valley  wavefront  error  by 
applying  an  update  to  the  influence  matrix  prior  to  using  the  influence  matrix  in  feedback 
control.  Two  experiments  are  presented  using  two  different  experimentally  derived 
influence  matrices.  The  first  experiment  uses  an  influence  matrix  with  influence  functions 
collected  at  +/-15V.  The  second  experiment  uses  an  influence  matrix  with  influence 
functions  collected  at  +/-30V. 

The  adaptive  influence  matrix  is  tested  using  the  SMT  AHM  segment.  The 
influence  matrix  is  updated  using  the  approach  described  in  Chapter  IV.  After  the 
influence  matrix  is  update  it  is  used  with  an  integral  feedback  controller  to  correct  the 
SMT  AHM  surface.  The  error  history  results  are  compared  to  an  integral  feedback 
controller  without  the  influence  matrix  update. 
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Figure  86  and  Figure  87  show  the  experimental  results  for  the  two  experiments. 
Both  experiments  show  the  influence  matrix  update  slight  improves  the  RMS  and  peak- 
to-valley  wavefront  error.  The  results  also  show  the  variation  in  performance  based  on 
the  influence  matrix.  The  +/-15V  influence  matrix  does  not  produce  results  as  good  as 
the  +/-30V  influence  matrix.  An  adaptive  influence  function  may  prove  useful  to  update 
the  influence  function  to  better  represent  the  physical  system. 
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Figure  86.  SMT  AHM  wavefront  error  history  with  and  without  influence 
matrix  update  using  integral  control,  +/-15V  influence  matrix 
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Figure  87.  SMT  AHM  wavefront  error  history  with  and  without  influence 
matrix  update  using  integral  control,  +/-  30V  influence  matrix 
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E.  CONCLUSIONS 

The  experimental  results  validated  the  woofer-tweeter  models  from  Chapter  III 
and  Chapter  IV.  The  results  showed  that  a  space  telescope  woofer-tweeter  system  is 
viable  and  can  improve  the  wavefront  error  attributed  to  the  primary  mirror  segment. 
The  experimental  woofer-tweeter  performance  was  limited  by  the  size  of  the  additional 
deformable  mirror  and  available  actuators,  however  a  24%  RMS  wavefront  error  and 
16%  peak  to  valley  wavefront  error  improvement  was  shown  using  the  woofer-tweeter 
gradient  controller.  The  woofer-tweeter  global  integral  controller  and  global  iterative 
constrained  minimization  controller  also  showed  wavefront  improvement  compared  to 
the  SMT  segment  constrained  minimization  controller.  Updating  the  SMT  segment 
influence  matrix  using  the  adaptive  influence  matrix  technique  also  showed  wavefront 
error  improvement. 
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VII.  SUMMARY,  CONTRIBUTIONS,  AND  FUTURE  WORK 


This  chapter  presents  a  final  summary  of  the  work  presented  in  this  dissertation. 
The  significant  contributions  of  this  research  are  summarized.  Recommendations  for 
future  research  are  also  described. 

A.  SUMMARY  OF  WORK 

The  goal  of  this  dissertation  was  to  develop  woofer-tweeter  control  techniques  for 
use  in  a  space  telescope.  In  order  to  do  this,  a  space  telescope  woofer-tweeter  model  was 
developed  to  simulate  woofer-tweeter  control  techniques.  The  woofer-tweeter  model  was 
based  on  the  Naval  Postgraduate  School’s  3-meter  SMT.  This  required  an  understanding 
the  SMT  woofer-tweeter  design  assumptions  and  an  understanding  of  the  SMT  design. 
The  result  was  a  woofer-tweeter  design  that  placed  a  deformable  mirror  at  the  SMT  exit 
pupil. 

Two  analytical  models  were  developed  based  on  woofer-tweeter  design.  The  first 
analytical  model  included  a  dynamic  model  of  the  SMT  with  a  static  influence  matrix 
model  for  the  additional  deformable  mirror.  The  second  model  included  a  quasi-static 
SMT  model  with  a  static  influence  matrix  deformable  mirror  model.  The  quasi-static 
woofer-tweeter  model  was  reduced  to  represent  a  single  SMT  segment  in  order  to  match 
the  experimental  setup.  The  models  were  updated  to  include  the  measured  primary 
mirror  residual  surface  error.  The  measured  surface  error  produced  simulated  results 
representative  of  the  physical  system. 

The  SMT  wavefront  control  was  simulated  beginning  from  segment  stacking 
through  fine  phasing  with  the  SMT  AHM  residual  error  as  an  initial  condition.  The 
simulated  results  showed  that  the  current  SMT  system  could  not  meet  the  fine  phasing 
requirements  due  to  the  residual  primary  mirror  surface  error.  Two  problems  were 
identified.  First,  the  Shack-Hartmann  wavefront  sensor  did  not  have  enough  spatial 
resolution  to  measure  the  high  frequency  wavefront  error.  Second,  the  residual  surface 
error  of  the  AHM  segment  was  outside  the  capture  range  of  the  phase  retrieval  algorithm. 
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These  results  supported  the  use  of  an  additional  deformable  mirror  to  improve  the 
residual  wavefront  error. 

Fine  phasing  improvement  approaches  were  presented  including  closed  loop 
woofer-tweeter  control  and  open  loop  woofer-tweeter  control.  The  simulated  and 
experimental  research  showed  that  the  closed  loop  control  approaches  outperformed  the 
open  loop  control  approaches  because  the  woofer-tweeter  influence  matrix  models  varied 
from  the  physical  hardware.  The  global  woofer-tweeter  control  techniques  reduced  the 
residual  wavefront  error  more  than  the  serial,  modal,  and  offload  control  techniques.  The 
woofer-tweeter  gradient  controller  produced  the  best  results  experimentally  reducing  the 
RMS  wavefront  error  by  24%  compared  to  the  best  SMT  AHM  corrected  wavefront 
error. 

Another  fine  phasing  improvement  presented  was  the  SMT  influence  matrix 
model  update.  The  SMT  AHM  influence  matrix  was  iteratively  updated  using  the 
gradient  of  a  cost  function  that  compared  the  difference  between  the  predicted  mirror 
response  and  the  measured  mirror  response.  Simulation  and  experimental  results  showed 
the  influence  matrix  update  technique  could  slightly  improve  the  performance  of  the 
SMT  AHM.  This  technique  may  be  advantageous  updating  the  influence  matrix  without 
recollecting  influence  functions  saving  precious  collection  time. 

Existing  closed  loop  woofer-tweeter  control  techniques  were  presented  and  a 
woofer-tweeter  gradient  control  technique  was  introduced.  Open  loop  woofer-tweeter 
control  techniques  were  also  introduced.  Finally  an  adaptive  influence  matrix  technique 
was  presented  as  a  means  to  improve  fine  phasing  residual  wavefront  error  by  improving 
the  SMT  AHM  segment  model. 

The  simulated  and  experimental  research  presented  in  this  dissertation  showed 
that  a  space  telescope  woofer-tweeter  system  is  feasible  and  can  improve  the  wavefront 
error  attributed  to  the  primary  mirror  segment.  Woofer-tweeter  space  telescope  control 
techniques  were  presented  and  experimentally  validated  showing  significant 
improvement  in  residual  wavefront  error. 
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B.  CONTRIBUTIONS 

Adaptive  optics  have  been  proposed  in  the  past  to  compensate  for  space  telescope 
primary  mirrors.  However,  the  current  state  of  the  art  has  not  addressed  wavefront 
control  techniques  for  a  woofer-tweeter  space  telescope  configuration.  The  contributions 
of  this  research  are  as  follows: 

1.  Woofer-Tweeter  Space  Telescope  Modeling  and  Simulation 

A  woofer-tweeter  space  telescope  model  was  developed  for  the  SMT  test  bed. 
This  model  built  upon  the  SMT  state  space  model  but  was  updated  to  include  an 
additional  deformable  mirror,  experimentally  derived  influence  matrices,  actuator 
models,  sensor  models,  full  state  observer,  and  residual  errors  associated  with  the  active 
primary  segments.  These  updates  enabled  more  realistic  results  during  simulation. 

The  SMT  incremental  wavefront  sensing  and  control  system  was  also  simulated 
for  the  first  time  from  deployment  through  fine  phasing  using  the  measured  SMT  AHM 
surface  figure  error  data.  The  coarse  alignment  and  coarse  phasing  steps  were 
successfully  controlled  using  LQG  state  feedback  controllers.  The  fine  phasing  step  was 
accomplished  using  a  quasi-static  adaptive  optic  integral  control  technique.  This 
simulation  showed  that  the  current  SMT  Shack-Hartmann  wavefront  sensor  lacked  the 
necessary  spatial  resolution  to  successfully  correct  the  SMT  AHM  residual  error.  The 
simulation  also  showed  that  the  SMT  AHM  residual  error  was  outside  the  capture  range 
of  phase  diversity  phase  retrieval  algorithms. 

The  model  was  reduced  to  a  single  segment  quasi-static  model  to  represent  the 
SMT  woofer-tweeter  experimental  test  bed  and  was  experimentally  validated  using  the 
SMT  test  bed.  The  model  was  used  to  simulate  fine  phasing  control  techniques  and 
evaluate  performance  prior  to  testing  the  controllers  on  the  test  bed. 

2.  Woofer-Tweeter  Space  Telescope  Control  Techniques 

Existing  closed  loop  woofer-tweeter  control  techniques  were  applied  to  a  space 
telescope  application.  The  techniques  were  simulated  and  experimentally  test  and 
performance  was  evaluated.  A  new  woofer-tweeter  gradient  control  technique  was 
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developed  using  SVD  modal  feedback.  The  gradient  controller  includes  singular  value 
weighting  that  automatically  weights  the  modal  coefficient  feedback  and  is  specific  to  the 
deformable  mirror  capabilities.  The  result  showed  the  gradient  controller  outperformed 
the  other  woofer-tweeter  controllers  in  simulation  and  experimental  tests.  The  SVD 
modal  gradient  controller  was  also  applied  to  the  SMT  AHM  mirror  and  showed 
improved  performance  over  the  integral  feedback  control  and  iterative  constrained 
minimization. 


3.  Woofer-tweeter  Space  Telescope  Experimental  Results 

This  research  showed  that  a  space  telescope  woofer-tweeter  system  is  feasible  and 
can  improve  the  wavefront  error  attributed  to  the  primary  mirror  segment.  The 
experiments  were  conducted  using  large  aperture  active  mirrors  designed  for  space 
applications.  The  errors  present  in  the  experimental  system  were  representative  of 
current  large  aperture  lightweight  optics.  The  research  showed  that  even  with  an 
undersized  additional  deformable  mirror  an  improvement  of  24%  RMS  wavefront  error 
was  achieved. 

4.  Experimental  Characterization  of  Actuated  Hybrid  Mirror 

A  SMT  actuated  hybrid  mirror  was  experimentally  characterized.  The 
experimental  results  showed  a  non-linear  response  of  the  SMT  segment  over  the  varying 
face  sheet  actuator  range.  Both  the  actuator  response  and  face  sheet  response  were 
shown  to  vary  from  the  linear  model.  This  information  provides  better  understanding  of 
the  behavior  of  these  mirrors  and  can  be  used  to  develop  new  control  strategies  to 
improve  the  performance  of  the  mirrors. 

5.  Adaptive  Influence  Matrix  Update 

An  adaptive  influence  matrix  update  technique  was  developed  in  order  to  improve 
the  performance  of  the  SMT  AHM  influence  matrix.  Experimental  results  showed  that 
using  the  updated  influence  matrix  with  an  integral  controller  improved  the  RMS 
wavefront  error  compared  to  a  controller  with  a  static  influence  matrix.  This  technique 
may  be  advantageous  for  updating  a  telescopes  active  primary  mirror  influence  matrix  on 
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orbit,  as  it  does  not  require  the  collection  of  influence  functions  for  each  actuator  saving 
valuable  collection  time. 

C.  RECOMMENDATIONS 

This  dissertation  showed  the  feasibly  of  using  an  additional  deformable  mirror 
with  a  large  active  primary  mirror  to  improve  the  residual  wavefront  error.  However  the 
simulated  and  experimental  results  were  based  on  the  performance  of  existing  mirror 
capabilities.  The  woofer-tweeter  fine  surface  control  approach  can  work  better  if  it  is 
considered  from  the  beginning  of  the  design  phase.  This  would  allow  the  woofer  and 
tweeter  performance  to  be  properly  traded  to  determine  the  individual  mirror  surface 
requirements  and  correction  capabilities. 

Future  woofer-tweeter  experimental  work  should  include  testing  a  properly  sized 
woofer-tweeter  space  telescope  system.  This  could  include  the  replacement  of  the 
existing  deformable  mirror  with  a  deformable  mirror  with  a  higher  actuator  count. 
Simulations  in  this  research  showed  a  minimum  of  543  actuators  were  need  to  reduce  the 
SMT  AHM  segment  residual  wavefront  error  of  A/ 1 5  waves  RMS.  Another  option  would 
be  to  demonstrate  a  woofer-tweeter  system  with  a  newer  large  AHM  with  less  residual 
surface  error. 

Another  recommendation  is  to  update  the  SMT  test  bed  to  include  a  deformable 
mirror  in  the  rear  optics.  This  could  add  confidence  to  the  woofer-tweeter  approach,  by 
demonstrating  a  more  realistic  implementation.  There  are  challenges  with  including  the 
deformable  mirror  in  the  optical  path  including  field  steering,  field  angle  magnification, 
and  the  optical  design.  The  SMT  Shack-Hartmann  wavefront  sensor  would  also  need  to 
be  upgraded  to  improve  spatial  resolution. 
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APPENDIX  A.  SIMULATION  RESULTS 


This  appendix  includes  the  error  history  and  wavefront  surface  plots  for  both  the 
feedback  control  techniques  and  open  loop  control  techniques  simulated  in  Chapter  V. 


A.  SMT  AIIM  SEGMENT  SIMULATED  RESULTS 

The  SMT  AHM  segment  was  simulated  using  the  156  FSAs  to  compare  the 
results  to  the  woofer-tweeter  configuration.  The  simulations  included  zero  mean 
Gaussian  plant  and  sensor  noise. 

1.  SMT  Segment  3  Static  Error 
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waves  (632.8  nm) 


2. 


SMT  Integral  Control 


SMT  Segment  Integral  Control  SMT  Segment  Integral  Control 


Figure  89.  SMT  integral  modal  control  for  different  modal  feedback 


SMT  Integral  Control 
RMS  WFE  =  0.324,  P-V  WFE  =  2.55 


Figure  90.  SMT  segment  integral  control  surface  error 
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3. 


SMT  Constrained  Minimization  Control 


SMT  Constrained  Minimization 
RMS  WFE  =  .251  P-V  WFE  =  2.25 


Figure  91.  SMT  constrained  least  squares  control  residual  surface  error 


4.  SMT  Gradient  Control 

This  simulation  only  adds  zero  mean  Gaussian  noise  with  a  standard  deviation  of 
0.5  volts  to  the  mirror  input.  To  represent  the  DM  fitting  capability  of  the  SMT  AHM 
segment  we  feedback  the  first  12  SVD  modes. 
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SMT  Segment  Gradient  Control 


Figure  92.  SMT  AHM  segment  gradient  control  wavefront  error  history 


SMT  Segment  Gradient  Control 
RMS  WFE  =  .174  P-V  WFE  =  1.3 


Figure  93.  SMT  AHM  gradient  control  residual  surface  error 


5.  Tweeter  Control 

The  tweeter  only  simulation  applies  a  correction  to  the  static  SMT  segment 
wavefront  using  the  BMC- 140  deformable  mirror.  The  process  noise  is  zero  mean 
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Gaussian  white  noise  with  a  standard  deviation  of  0.25  volts.  The  noise  is  only  applied  to 
the  tweeter  mirror.  The  simulation  uses  an  iterative  constrained  least  squares  controller 
with  an  update  gain  of  0.2  and  an  actuator  saturation  limit  of  +/-  80%  of  the  mirror’s 
biased  maximum  voltage. 


BMC  Constrained  Minimization 


Figure  94.  Tweeter  iterative  constrained 


BMC  Constrained  Minimization 


minimization  wavefront  error  history 
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BMC  Constrained  Minimization 
RMS  WFE  =  .209  P-V  WFE  =  2.68 
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Figure  95.  Tweeter  iterative  constrained  minimization  control  residual  surface 

error 

B.  WOOFER-TWEETER  CLOSED-LOOP  SIMULATED  RESULTS 

This  section  includes  the  woofer-tweeter  simulated  results.  The  woofer  in  this 
case  is  the  SMT  AHM  segment  and  the  tweeter  is  the  BMC  deformable  mirror  with  61 
active  channels. 

1.  Serial  Constrained  Minimization  Control 

The  simulation  runs  100  times  and  the  mean  RMS  and  peak-to-valley  wavefront 
error  are  plotted.  The  two-step  correction  uses  an  iterative  constrained  least  squares 
routine  to  minimize  the  wavefront.  First  the  SMT  AHM  segment  applied  correction,  and 
then  the  61  actuators  of  BMC- 140  deformable  mirror  applied  correction.  An  update  gain 
of  0.2  is  used  for  both  controllers  and  the  saturation  limit  is  set  to  +/-  30  volts  for  the 
SMT  and  the  BMC  deformable  mirror  where  the  BMC  deformable  mirror  inputs  are 
normalized  to  match  the  SMT  input  voltages.  Process  noise  is  zero  mean  Gaussian  with  a 
standard  deviation  of  0.5  volts  on  the  SMT  AHM  and  0.25  volts  on  the  BMC- 140  input. 
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Serial  CLS  Simulated  Control 
RMS  WFE  =  0.245 


Serial  CLS  Simulated  Control 
P-V  WFE  =  1.79 


Figure  96.  Serial  iterative  constrained  least  squares  woofer-tweeter  controller 

wavefront  error  history 


Serial  CLS  Control 
RMS  WFE  =  .245  P-V  WFE  =  1.79 


0  0 


Figure  97.  Serial  iterative  constrained  least  squares  controller  residual  surface 

error 
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2. 


Global  Control 


The  global  integral  control  simulations  run  100  times.  Results  are  shown  for 
integral  control,  iterative  constrained  minimization  control  and  gradient  control.  The 
figures  in  this  section  show  the  mean  RMS  and  peak-to-valley  wavefront  error.  The  error 
bars  show  the  standard  deviation.  The  process  noise  is  zero  mean  Gaussian  white  noise 
with  a  standard  deviation  of  0.5  volts. 

a.  Global  Integral  Control 


Woofer-Tweeter  Integral  Control 


Woofer-Tweeter  Integral  Control 


Figure  98.  Woofer-tweeter  integral  control  wavefront  error  history 
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Woofer-Tweeter  Integral  Control 
RMS  WFE  =  .184  P-V  WFE  =  1.51 


0.6 


Figure  99.  Woofer-tweeter  integral  control  residual  surface  error 


b.  Global  Iterative  Constrained  Minimization  Control 


Global  Iterative  Constrained  Minimization  Global  Iterative  Constrained  Minimization 


Figure  100.  Woofer-tweeter  global  iterative  minimization  control  wavefront 

error  history 
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Global  Iterative  Constrained  Minimization 
RMS  WFE  =  0.174,  P-V  WFE  =  1.4 


Figure  101.  Woofer-tweeter  global  iterative  minimization  control  residual 

surface  error 


c.  Global  Gradient  Control 

The  first  100  SVD  modes  are  used  as  feedback. 


Woofer-Tweeter  Gradient  Control 


Woofer-Tweeter  Gradient  Control 


Figure  102.  Woofer-tweeter  gradient  control  wavefront  error  history 
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Woofer-Tweeter  Global  Gradient  Control 
RMS  WFE  =  .154  P-V  WFE  =  1.39 


Figure  103.  Woofer-tweeter  gradient  control  residual  surface  error 


3.  Offload  Control 

The  offload  control  approach  is  simulated  using  an  integral  controller  and  a 
constrained  minimization  controller.  The  control  simulation  runs  100  times.  The  figures 
show  the  mean  RMS  and  peak-to-valley  wavefront  error.  The  error  bars  show  the 
standard  deviation.  The  process  noise  is  zero  mean  Gaussian  white  noise  with  a  standard 
deviation  of  0.5  volts. 
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Offload  Integral  Control 


Woofer -Tweeter  Offload  Integral  Control 
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Figure  104.  Woofer-tweeter  offload  integral  control  wavefront  error  history 


Woofer -Tweeter  Offload  Integral  Control 
RMS  WFE  =  .218  P-V  WFE  =  2.53 
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Figure  105.  Woofer-tweeter  offload  integral  control  residual  surface  error 
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b.  Offload  Iterative  Constrained  Minimization  Control 


Woofer-Tweeter  Offload  CLS  Control  Woofer-Tweeter  Offload  CLS  Control 


Figure  106.  Woofer-tweeter  offload  iterative  constrained  least  squares 
wavefront  error  history 


Woofer-Tweeter  Offload  CLS  Control 
RMS  WFE  =  .195  P-V  WFE  =  2.01 
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Figure  107.  Woofer- tweeter  offload  constrained  least  squares  residual  surface 

error 
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4. 


Modal  Control 


This  section  presents  simulated  results  for  two  modal  controllers.  The  first  is  the 
modal  coefficient  controller  and  the  second  is  the  distributed  modal  controller.  The 
control  simulations  run  100  times.  The  distributed  modal  control  simulation  is  run  for 
three  separate  scenarios  described  in  this  section.  The  figures  show  the  mean  RMS  and 
peak-to- valley  wavefront  error.  The  error  bars  show  the  standard  deviation.  The  process 
noise  is  zero  mean  Gaussian  white  noise  with  a  standard  deviation  of  0.5  volts. 


a.  Modal  Coefficient  Control 

This  modal  coefficient  controller  uses  Zemike  coefficients  as  the  common  basis 
for  the  two  SMT  segment  and  deformable  mirror.  The  woofer  is  assigned  the  first  7 
modes  while  the  tweeter  is  assigned  the  Zemike  modes  8-24.  The  residual  wavefront 
includes  the  high  order  modes  above  24  as  well  as  the  deformable  mirror  fitting  error. 
The  fitting  error  is  the  error  associated  with  the  mirror  not  reproducing  the  Zemike  mode 
exactly. 


Woofer-Tweeter  Modal  Coefficient  Control  Woofer-Tweeter  Modal  Coefficient  Control 


Figure  108.  Woofer-tweeter  modal  coefficient  control  wavefront  error  history 
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Woofer-Tweeter  Modal  Coefficient  Control 
RMS  WFE  =  .263  P-V  WFE  =  2.01 


Figure  109.  Woofer- tweeter  modal  coefficient  control  residual  surface  error 
b.  Distributed  Modal  Control 

Distributed  modal  control  simulations  are  run  for  three  scenarios.  The  first 
scenario  completely  removes  the  influence  of  the  woofer  from  the  tweeter.  The  second 
scenario  removes  the  first  7  modes  of  the  woofer  from  the  tweeter,  and  the  third  scenario 
removes  the  first  1 1  modes  of  the  woofer  from  the  tweeter.  The  figures  in  this  section 
show  the  steady  state  surface  error  and  the  RMS  and  peak-to-valley  wavefront  error 
history. 
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Woofer-Tweeter  DMC  Control  Woofer-Tweeter  DMC  Control 


Figure  110.  Woofer-tweeter  distributed  modal  control  wavefront  error  history 
woofer  removed  from  tweeter 


Figure  111.  Woofer-tweeter  distributed  modal  control,  woofer  removed  from 
tweeter,  residual  surface  error 
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Woofer-Tweeter  DMC  Control  Woofer-Tweeter  DMC  Control 


Figure  1 12.  Woofer-tweeter  distributed  modal  control  wavefront  error  history, 
mode  1-7  removed  from  tweeter 


Figure  113.  Woofer-tweeter  distributed  modal  control,  mode  1-7  removed  from 

tweeter,  residual  surface  error 
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Woofer-Tweeter  DMC  Control 
RMS  WFE  History 
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Figure  114.  Woofer-tweeter  distributed  modal  control  wavefront  error  history, 
mode  1-11  removed  from  tweeter 


Woofer-Tweeter  DMC  Integral  Control 
RMS  WFE  =  .218  P-V  WFE  =  1.58 


Figure  115.  Woofer-tweeter  distributed  modal  control,  mode  1-11  removed 
from  tweeter,  residual  surface  error 
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APPENDIX  B.  EXPERIMENTAL  RESULTS 


Included  in  this  appendix  are  the  SMT  test  bed  woofer-tweeter  controller 
experimental  results.  The  section  includes  wavefront  error  history  plots  presented  as 
RMS  wavefront  error  and  peak-to-valley  wavefront  error  measured  in  waves  (632.8  nm). 
A  trend  line  is  fitted  to  the  experimental  data  in  the  error  plots.  The  open  loop  control 
results  are  presented  as  before  and  after  surface  plots. 

A.  SMT  AHM  SEGMENT  EXPERIMENTAL  RESULTS 

This  section  includes  experimental  results  for  correcting  the  SMT  AHM  segment. 
This  includes  integral  control,  iterative  constrained  minimization  control,  gradient 
control,  and  tweeter  control. 

1.  SMT  Integral  Control 


SMT  Modal  Integral  Control 
P-V  WFE  History 


10 

Iterations 


Figure  1 16.  SMT  AHM  modal  feedback  integral  control 


153 


waves  (632.8nm)  w  waves  (632.8nm) 


2, 


SMT  Iterative  Constrained  Minimization  Control 


SMT  Iterative  Constrained  Minimization  SMT  Iterative  Constrained  Minimization 


Figure  117.  SMT  AHM  iterative  constrained  minimization  control 


SMT  Gradient  Control 

SMT  Gradient  Control  SMT  Gradient  Control 


Figure  118.  SMT  AHM  gradient  control 
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Tweeter  Correction 


Tweeter  Iterative  Constrained  Minimization  Tweeter  Iterative  Constrained  Minimization 


Figure  119.  SMT  AHM  tweeter  correction 


SMT  Segment  3 
RMS  WFE  =  .5603 


SMT  Segment  3  -  Global  Gradient 
RMS  WFE  =.1756 


o  o 


Figure  120.  Comparison  of  static  SMT  wavefront  error  before  and  after 
correction  using  deformable  mirror  with  global  gradient  control 
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B.  WOOFER-TWEETER  CLOSED  LOOP  EXPERIMENTAL  RESULTS 

This  section  includes  woofer-tweeter  closed  loop  control  results  including  the  2- 
step  iterative  constrained  least  squares  controller,  global  controllers,  and  distributed 
modal  controller. 

1.  Serial  Control  Constrained  Minimization  Control 


Serial  Woofer  &  Tweeter  Control 


Serial  Woofer  &  Tweeter  Control 


Figure  121.  Serial  integral  control 


Serial  Woofer -Tweeter 
Constrained  Minimization  Control 
RMS  Wavefront  Error 


Serial  Woofer-Tweeter 
Constrained  Minimization  Control 


Figure  122.  Serial  constrained  minimization  control 
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2. 


a. 


Global  Control 
Global  Integral  Control 


Global  Integral  Control 


Global  Integral  Control 


Figure  123.  Global  integral  control 


b. 


Global  Iterative  Constrained  Minimization  Control 


Global  Iterative  Constrained  Minimization 
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Figure  124.  Global  iterative  constrained  minimization  control 
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c. 


Gradient  Control 


Woofer-Tweeter  Gradient  Control 


Woofer-Tweeter  Gradient  Control 


Figure  125.  Woofer-tweeter  gradient  control 


3.  Offload  Control 

Experimental  offload  control  techniques  are  shown  in  this  section  using  an 
integral  controller  and  a  constrained  minimization  controller. 

a.  Offload  Integral  Control 


Figure  126.  Offload  integral  control 
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b.  Offload  Iterative  Constrained  Minimization  Control 


Woofer-Tweeter  Constrained  Minimization 
Offload  Control 


Woofer-Tweeter  Constrained  Minimization 
Offload  Control 


Iteration 


Iteration 


Figure  127.  Offload  iterative  constrained  minimization 


4.  Distributed  Modal  Control 

This  section  presents  the  distributed  modal  control  experimental  results  for  two 
cases.  The  first  case  removes  the  first  7  SVD  woofer  modes  from  the  tweeter  and  the 
second  case  removes  the  first  1 1  SVD  woofer  modes  from  the  tweeter.  The  experiments 
were  run  for  different  modal  feedback  settings.  These  settings  removed  higher  frequency 
modes  from  the  feedback. 
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Woofer-Tweeter  Distributed  Modal  Control  Woofer-Tweeter  Distributed  Modal  Control 


Figure  128.  Distributed  modal  control,  SVD  woofer  modes  1-7  removed  from 

tweeter 
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Figure  129.  Distributed  modal  control,  SVD  woofer  modes  1-11  removed  from 

tweeter 


C.  OPEN  LOOP  EXPERIMENTAL  RESULTS 

This  section  includes  the  open  loop  experimental  results  using  the  SMT  AHM 
segment  and  the  woofer-tweeter  experimental  setup.  The  cost  function  (4.40)  from 
Chapter  IV  is  used  and  presented  again.  Three  minimizations  are  performed  with 


160 


different  weighting  coefficients,  shown  in  Table  3.  The  minimization  was  solved  using 
the  MATLAB  function  fmincon.m. 

arg min  J (u)  s*4(c,  ({ >}r  {</>})  +  c2  (ma x({0})  -  min({0}))  +  c3  ({vA  }T  {vA  })  +  c4  ({u}r  {«})) 
subject  to:  lb  <  [u]  <  ub 


Table  4.  Cost  function  coefficients 


cl 

c2 

c3 

c4 

Cost  Function  1 

1 

0 

0 

0 

Cost  Function  2 

1 

2 

0 

0 

Cost  Function  3 

1 

1 

0 

5.00E-04 

The  experiments  were  completed  10  times  and  the  final  wavefront  error  results 


were  averaged.  The  surface  plots  below  show  the  resulting  surface  error  for  a  single 
experiment. 


SMT  Initial  WFE 
RMS  WFE:  0.559  P-V  WFE:  3.52 


SMT  Open  Loop  Control,  Cost  1 
RMS  WFE:  0.559  P-V  WFE:  3.52 
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Figure  130.  SMT  open  loop  minimization,  cost  function  1 
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SMT  Initial  WFE 
RMS  WFE:  0.599  P-V  WFE:  4.24 


SMT  Open  Loop  Control,  Cost  2 
RMS  WFE:  0.402  P-V  WFE:  3.22 


Figure  131.  SMT  open  loop  minimization,  cost  function  2 
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SMT  Open  Loop  Control,  Cost  3 
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Figure  132.  SMT  open  loop  minimization,  cost  function  3 
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SMT  Initial  WFE 
RMS  WFE:  0.533  P-V  WFE:  3.44 


Woofer-Tweeter  Open  Loop  Control 
Cost  1,  RMS  WFE:  0.319  P-V  WFE:  2.99 
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Figure  133.  Woofer-tweeter  open  loop  minimization,  cost  function  1 


SMT  Initial  WFE 
RMS  WFE:  0.599  P-V  WFE:  4.24 


SMT  Open  Loop  Control,  Cost  2 
RMS  WFE:  0.402  P-V  WFE:  3.22 
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Figure  134.  Woofer-tweeter  open  loop  minimization,  cost  function  2 
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SMT  Initial  WFE  Woofer-Tweeter  Open  Loop  Control 

RMS  WFE:  0.601  P-V  WFE:  3.85  Cost  3,  RMS  WFE:  0.358  P-V  WFE:  2.9 


Figure  135.  Woofer-tweeter  open  loop  minimization,  cost  function  3 
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